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Abstract 

The mechanisms of spontaneous and induced radiation produced by relativistic 
particles passing through natural and photonic crystals are reviewed. The theory 
of Volume Free Electron Lasers based on spontaneous radiation in natural and 
photonic crystals is presented. 
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Introduction 



Current development of Free Electron Lasers (FELs) operating in various 
ranges of wavelengths - from microwave to optical and even X-ray - which is 
carried out in some research centers gives reason to return to the problem of 
the possibility to design the FEL on the basis of the interaction between rel- 
ativistic electron (positron) beams and either natural or artificial (photonic) 
crystals. In the X-ray range of wavelengths, operation of this type of crystal 
X-ray FELs is based on several spontaneous radiation mechanisms: parametric 
X-ray radiation |17f85fl23fl24fl39] . diffracted radiation from relativistic chan- 




neled particles |17|19|22|139|142|143|144|145|153] . generation of X-ray radia- 
tion by particles moving in a crystal undulator p^l7f23|139|155fl56|157|158] . 
It was shown that despite serious difficulties, the law of radiative instability 
of relativistic beams passing through the crystal, which was discovered in |19] , 
enables achieving the generation threshold for induced X-ray radiation even at 
the beam current density j ~ 10^ A/cm^, while the current density required 
for X-ray laser generation in crystals by means of channeling X-ray radiation 
or radiation in a crystal undulator is j ~ 10^^ ^ 10^'^ |139] . 

Below, a brief review of the theory of the crystal FEL (based on either a 
natural or a photonic crystal) is presented. The review emphasizes the general 
significance of the law of radiative instability discovered in [19], which finally 
led to formulating the idea of the possibility to develop a new type of FELs, 
called Volume Free Electron Lasers 125) . 
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1 Diffraction phenomena accompanying spontaneous and stimu- 
lated radiation from relativistic particles in crystals ( |139] ) 



The emission of photons by relativistic particles in media has been calling 
attention for a long time. The main reason why this phenomenon arises such 
interest is because a wide variety of tasks can be solved by using different 
radiation mechanisms, such as bremsstrahlung, transition and Cherenkov ra- 
diations and so on. The last decades have seen a growing interest in the study 
of radiation by relativistic particles in both natural and artificial (photonic) 
crystals. A number of new radiation mechanisms associated with a periodic 
structure of crystals have been considered theoretically and confirmed exper- 
imentally pil7fi:-i9fl4()l[TiT] . 

All characteristic properties of radiation are in this case determined by the 
periodic structure of a crystal. The medium can influence radiation processes 
under passing of relativistic charged particles through crystals in two ways. 
First of all, it is well known that the trajectory of a charged particle inci- 
dent on a crystal at a small angle relative to crystallographic planes or axes is 
formed by a series of grazing collisions with the atoms of a crystal. As a result, 
the particle moves in an averaged potential of crystallographic planes or axes. 
In this case we speak about the channeling phenomenon and can consider 
the motion of channeled particles as the motion inside a potential well - one- 
dimensional (for plane channeling) or two-dimensional (for axial channeling). 
According to quantum mechanics a particle moving inside a potential well has 
a discrete energy spectrum, which is the spectrum of the particle transverse 
motion. Consequently, such particles can be considered as one-dimensional 
or two-dimensional atoms (oscillators) characterized by a spectrum of bound 
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states (zones) of transverse energy e„, Sf. The number of bound states and 
their characteristics depend on the longitudinal particle energy. One can con- 
clude that many phenomena observed for ordinary atoms will manifest them- 
selves when channeled particles pass through crystals. It is obvious that thus 
excited atoms should emit photons with the energy equal to the difference be- 
tween atomic state energies and Sf. The frequency of transition = EnSf 
and depends, in a laboratory frame, on the total particle energy. By analogy 
with an ordinary moving oscillator, the frequency of the emitted photons is 
evaluated by the Doppler effect and determined by the following expression: 



where 6 is the radiation angle, /3 = u/c, u is the longitudinal particle velocity, 
n{u) is the refraction index of the photon with a frequency u in the medium. 

The relativistic oscillator can be formed not only by an unperturbed crystal 
channel, but also by an external ultrasonic or laser wave which propagates in 
the crystal, creating a bent crystal channel [T3|17fl23fl39fl55fl56fl57fl58] 




On the other hand, when the wavelength of the emitted photons is of the or- 
der of the interplanar spacing of atoms in the crystal, radiation diffraction can 
essentially modify the photon state |123] . In this case the radiation process 
is characterized by several indices of refraction ni{u) dependent on the direc- 
tion of the photon momentum, which, in turn, leads to the modification of all 
mechanisms of radiation formation by relativistic particles in the X-ray range 
of spectrum. For example, radiation at a large angle relative to the direction 
of particle motion becomes possible. As a result, the diffraction pattern char- 
acterizing a given crystal is formed. The analysis of dielectric properties of a 



CO = 




(1) 



1 — f3n{uj) cos 6 
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crystal under diffraction conditions shows that at least one of several indices 
of refraction ni{u), characterizing the crystal under this condition, becomes 
greater than unity within a frequency interval. As a consequence, the Vavilov- 
Cherenkov condition can be fulfilled. In this case spontaneous and induced 
X-ray radiation analogous to optical Cherenkov radiation appears as well as 
diffracted transition radiation [123]. Spontaneous quasi- Cherenkov radiation, 
currently referred to as parametric (quasi- Cherenkov) X-ray radiation (PXR), 
is thoroughly investigated both theoretically and experimentally (see [85]). 

Under diffraction conditions, the radiation of a relativistic oscillator also mod- 
ifies essentially. Now the periodic structure of a crystal affects both the par- 
ticle motion and the particle state. This leads to the formation of diffracted 
radiation of oscillator (DRO), which cannot be reduced to the sequence of 
two independent processes: radiation by oscillator and diffraction of radiated 
photons. In this case the process of photon emission and its diffraction are 
developing simultaneously and coherently and result in the radiation with 
new properties. Radiation produced by channeled particles is often called the 
diffracted channeling radiation (DCR). 

These types of X-radiation (DRO, DCR), also associated with a change in 
the indices of refraction under diffraction conditions, was considered in [19] 
and then in [TT|12f2U] . Currently, there is an increased interest in DCR (see 
[142fl43fl44|145|153j ). If, in the absence of diffraction, the X-ray spectrum 
of the oscillator is determined by the complex Doppler effect {n{uj) < 1), 
then, under diffraction, the index of refraction can become greater than unity 
and, consequently, the anomalous Doppler effect is possible. In this case the 
photon emitted by the oscillator is accompanied by the excitation of the os- 
cillator itself. This is one of the important features of diffracted radiation of 
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the oscillator (DRO). 



So, the modification of refractive properties in periodic media under diffraction 
leads to the appearance of two types of radiation with angular distribution 
forming a diffraction pattern determined by the parameters of the periodic 
medium. Depending on the excited reflex, the spectrum of PXR and DRO for 
crystals with lattice parameters of the order of A lies in the range of X-ray 
and even higher frequencies. 

It is well known that in an amorphous medium, the ordinary Cherenkov radia- 
tion can be considered as a specific case of radiation of the oscillator with zero 
eigenfrequency |T8] . Similarly, in a periodic medium, the frequency of PXR can 
be expressed by ([T]) in the specific case of Qnf = 0, i.e., 1 — Prii^uj) cosO = 0, 
where ni{u) is the refraction index of the crystal under diffraction conditions. 



2 Dispersion characteristics of parametric X-ray radiation (PXR) 
and diffracted radiation of oscillator (DRO) 

Under diffraction the index of refraction depends on the direction of particle 
motion and the frequency of radiated photon, therefore equation ([1]) deter- 
mining the radiated photon spectrum is the equation with several solutions 
[TT1IT91I201I2T1I22]. 



For example, let us consider the case of two-wave diffraction, when the diffrac- 
tion condition is fulfilled only for a reciprocal lattice vector r. It means that 
two strong waves with wave vectors k and kf = k + f are excited under diffrac- 
tion. For the simplicity of analysis of photon frequencies, let us represent ([T]) 
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in the form: 



u — Q 



nlu] 



nf 



Lj/3 COS 9' 



(2) 



In this case the index of refraction in a crystal under diffraction conditions is 
characterized by two dispersion branches rii 2. Using a well-known expression 
for rii 2, one can rewrite ^ as follows: 



1 ^^n/(l - S) 



l3cos9 ujb/3 cos 9 



(3) 



1 - 
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:i + (3,)-Af3,— T 
\9o\ 



\ 



(A - 1) - Af3, 



m 



9o 



for the diffracted radiation of the oscillator (DRO) and 



1 



(3 cos 9 



(4) 



1 - 



;i + /3i)-A/3i^T 
1^0 1 



l^ol 



for parametric quasi- Cher enkov radiation (PXR), generated by a particle pass- 
ing through a crystal at constant velocity. 



We have introduced the following notations: 

2{kT) + 



a 



is the deviation from the exact Bragg conditions, 



' da\ 



ot = au^g + I — I [uj -ujb) 



where 



A 



> 
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5 = oj — ojb/ojb^ = T^c^/2|r/3| is the Bragg frequency, corresponding to 
a = 0, /3i = kz/k;^ + Tz is the geometry factor of diffraction asymmetry, the 
2;-axis is chosen as a normal to the target surface directed inward the crystaL 
Let us assume that the particle with a mean velocity u moves along the z- 
axis, /3 = u/c, (7o; Qt are the coefficients in a series expansion in terms of 
the reciprocal lattice vectors of the crystal susceptibility. For simplicity, we 
shall assume that a crystal is center-symmetric and absorption is neglected, 
uj\ = Aireno/mc is the Langmuir frequency of the medium. 



Let us consider in detail the conditions of the existence of parametric X-ray 
radiation. Indeed, for the existence of this radiation it is sufficient that relation 
(jl]) be fulfilled at least for one of the two indices of refraction characterizing 
the crystal at a given frequency. On the left-hand side of (jl]), there is a term 
greater than unity (/3 < 1 and | cos 6*1 < 1). Consequently, for the existence of 
the solution of equation (jl]), the expression between the braces on the right- 
hand side of (jl]) should be less than zero. Far from Bragg conditions 6 — ?> oo, 
and we transit to a well-known case of amorphous medium., i.e., to the index of 
refraction n{u)) = 1 — bj\l'2ij^ < 1 for any frequencies within the X-ray range. 
As a result, in this range Cherenkov radiation is impossible. The analysis of 
(jl]) near Bragg condition |a| < {qqI shows that the expression between the 
braces is always positive for one dispersion branch corresponding to the sign 
(— ); consequently, the fulfilment of (jlj) is impossible. For the second branch 
corresponding to the sign (-|-), this expression can be negative at |a| > \go\. 
For example, in the case of 







gr 






9o 
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we can approximately write 



/3 cos 6 



1 



1 - 




and, obviously, the fulfilment of the Cherenkov condition is possible for Laue 
diffraction case > 0) at frequencies for which 



In Bragg diffraction case < 0), the Cherenkov condition can be fulfilled not 
only for one dispersion branch but even for two branches at the degeneration 
point. 

Comparison of (jlj) and ([3]) shows that the relation ([3]) can be satisfied for 
both dispersion branches at Qnf > 0, i.e., for radiation accompanied by the 
transition of the channeled particle to a lower energy level (e^ > Ef). This 
means that photons with two different frequencies are radiated at a given 
angle. In this situation we observe a normal complex Doppler effect. In this 
case, equation ([3]) is satisfied for radiation angles larger than the angle of 
parametric (quasi- Cherenkov) radiation. At the same time, the fulfilment of 
leads to the strong limitation of the particle energy, the radiation angle 
and the value of deviation from the exact Bragg condition a for the radiation 
accompanied by the oscillator excitation Q^f < (£„ < Sf). In this case, 
according to ([3]), radiation of a photon (anomalous Doppler effect) with the 
wave vector directed at a smaller angle relative to the direction of the particle 
motion than the angle of parametric (quasi-Cherenkov) radiation is possible 
for one of the dispersion branches. 
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Equation ([3]) can be analytically treated for a specific case when 



(5) 



where 



2 I l^ol 



/3i - 1 - A(3, 



\9o\ 

















9o 





(6) 



Let us represent ([3]) in the form: 



IjJ 



(7) 



In view of ([5]), ([7]) splits into two independent equations corresponding to the 
upper and the lower radiation branches in the absence of diffraction. 



col 
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(9) 



Neglecting the dependence of /3i on w, one can obtain the frequency solutions 
of (IHD and (E]) clS cl function of a radiation angle: 
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where 



X 



UJr. 



n 



2n' 1-/3 cos 9 

|- 1 7]| I COS ^ — I T_L I sin 6 cos (p 



T 

~2 



||f]|| COS0 — |rj_| sin0cos<y9| 



7]l and tl are the projections of the reciprocal lattice vector onto the direction 
of the particle mean velocity and onto the plane perpendicular to particle 
velocity, respectively; is the angle between r_L and k^. The dependence 
iij = ijj(Q) for the case of symmetric diffraction /3i = 1 is shown in Figured] 




S, ^ tfr ^ 

Figure 1. The dependence uj = uj{6) for the case of symmetric diffraction /3i = 1 

According to Figure [1] and (|TOl) and (|TT1) . the spectrum of radiated photons 
essentially modifies under diffraction. In the absence of diffraction, each ra- 
diation branch splits into two subbranches. Thus, diffraction results in the 
excitation of an additional branch in the complex Doppler effect with the fre- 
quency close to the Bragg frequency ujb, and in the formation of a radiation 
non-transparency range in the angular distribution AO = 9^~9_ = lO^^^lO"^ 
rad. When the angle 6 changes, the solution is first realized for one dispersion 
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branch and then for the other. Figure [T] shows the angular range in which 
|tt| < Ifi'ol- It should be noted that the radiation frequency of the additional 
diffraction branch changes a little with the change in the radiation angle. 
As a result, the angular range, in which |a| < \go\, may considerably exceed 
the ordinary angular interval, characterizing diffraction of an X-ray external 
monochromatic wave, when A6 is of the order of several angular minutes (in 
Figure [H this is the angular interval A6 = + A6^). As it is rather com- 
plicated to obtain the analytical solution, the numerical calculation of the 
dependence of lal/l^fol on the radiation angle 9 near the Bragg angle 9b was 
made for the oscillator moving along the crystal direction (110) and photon 
diffraction by crystallographic planes (400) in Si. According to these calcula- 
tions, the magnitude of A9 weakly depends on the energy and eigenfrequency 
of the oscillator and retains within the interval 10~^ — 10^'^. In Figure [21 the 
magnitude of lal/l^fol is shown as a function of the angle 9 at the following pa- 
rameters of the oscillator: Q = 1 eV, 7 = 2- 10^. As one can see, A9-. ~ 4-10^^ 
rad, A9+ = 5.2 ■ 10"^ rad and the total interval \a\ <\go\, A9 = 9.2 ■ lO'^ rad. 



'J 

-15 



■T - - ^- ^- 





Figure 2. 



The values of the angles 9± at a = are determined from the equation 



= - Wo f 1 ± — I • 
V 90 



'12) 
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A similar picture is also observed for the second branch (see (fTT]) ). So, if 
uo 7^ ujb, and the second term in the radicand of (fTTl) is negligible, then, for 
example, un+ = <^b and = Uq at Uq < ub- This means that, as in the 
previous case, we have the excitation of the wave with the frequency close 
to the Bragg frequency in addition to the solution far from diffraction at a 
given radiation angle 6. As the calculation shows, the magnitude of a does not 
depend on the frequency of the diffracted wave and keeps practically constant, 
being determined by the oscillator eigenfrequency within the whole interval of 
radiation angles corresponding to the fulfilment of the condition . Although 
the parameter a depends on the eigenfrequency of the oscillator fl and equals 
zero at 



it remains small (|a| < \go\) within rather a wide interval of eigenfrequencies. 
For example, in the case considered above Q± = 0.092 eV and \a\ < \go\ for 
the interval i^+min = 0.08 eV and fi+max = 0.15 eV. 

It should be noted that the equality (fTTj) also has the solution for the negative 
eigenfrequency of the oscillator (at the frequency close to w^) that corresponds 
to the anomalous Doppler effect, i.e., the radiation of the oscillator is accom- 
panied by its excitation. Such a process is possible because under diffraction 
condition the index of refraction can be greater than unity. 

The analysis of the dispersion expression for radiation propagating at a large 
angle and for arbitrary geometry was made in [24|. Here we discuss the spectra 
of DRO and PXR only in the two- wave diffraction case. However, due to crystal 
symmetry, the diffraction condition can be satisfied for many waves, that is. 




(13) 
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the case of multi-wave diffraction can be realized. In this case several indices of 
refraction ni{u) corresponding to different dispersion branches can be greater 
than unity. It appears that the possibility of new effects in radiation, such 
as the effect of excitation of radiation in a roundabout way takes place. This 
means that the PXR intensity in a diffraction peak may differ from zero even in 
the case when a given reflection is forbidden because of the lattice symmetry. 
The particular properties in the angular distribution of radiation are observed 
in the vicinity of the point of degeneration of dispersion branches. 



2.1 General expression for spectral- angular distribution of radiation gener- 
ated by a particle in a photonic crystal 



Both the spectral-angular density of radiation energy per unit solid angle W^u) 
and the differential number of emitted photons dNfi^^uj = l/hu ■ Wfi^^ can be 
easily obtained if the field E{f,u) produced by a particle at a large distance 
r from the crystal is known [T7] 



er 
4^ 



E(f,u) 



(14) 



The vinculum here means averaging over all possible states of the radiating 
system. In order to obtain E{f,u), Maxwell's equation describing the interac- 
tion of particles with the medium should be solved. The transverse solution 
can be found with the help of Green's function of this equation, which satisfies 
the expression: 

G = Go + Go^ie-l)G, (15) 
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Go is the transverse Green's function of Maxwell's equation at £ = 1. It is 
given, for example, in [6]. 

Using G, we can find the field we are concerned with 

Enir,u) = J Gne{r,r,u)^joeir,uj)(fr', (16) 

where n,e = x,y, z,joe{'f^,u}) is the Fourier transformation of the e-th compo- 
nent of the current produced by a moving beam of charged particles (in the 
linear field approximation, the current is determined by the velocity and the 
trajectory of a particle, which are obtained from the equation of particle mo- 
tion in the external field, by neglecting the influence of the radiation field on 
the particle motion). Under the quantum-mechanical consideration the cur- 
rent Jo should be considered as the current of transition of the particle-medium 
system from one state to another. 

According to [T7|85] . Green's function is expressed at r — > oo through the 
solution of homogeneous Maxwell's equations E^\f,uj) containing incoming 
spherical waves: 

„jfcr 

limG„e(r,r-,a;) = — E <4:^^*(^" (17) 
^ s 

r — )■ oo 

where e** is the unit polarization vector, s — 1, 2, ± ± fc. 

If the electromagnetic wave is incident on a crystal of finite size, then at r — )■ oo 

pikr 

El^-\r,uj) = e^e^'^'V const , 
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and one can show that the relation between the solution El^~^ and the solution 
of Maxwell's equation E^^\k,uj) describing scattering of a plane wave by the 
target (crystal), is given by: 

E'S-> = E^(j-) (18) 

k —k ^ ' 



Using flTBl) . we obtain 



En{fM = '—'-^T.<j Ef~>{f,u)Ur-,co)d'r'. (19) 



As a result, the spectral energy density of photons with polarization s can be 
written in the form: 



s 

n,Lu 



Et^ ■** (r , u)jo{f, u)d^r 



jo{r,uj) = / e'^%{f,io)dt = eQ / e'^'v{t)S{r- f{t))dt 



(20) 
(21) 



where eQ is the charge of the particle, v{t) and r{t) are the velocity and the 
trajectory of the particle at moment t. By introducing fl2T|) into f l20|) we get 



dm 



E\r^'*{r{t),uj)v{t)e^H 



{22) 



Integration in fl22|) is carried out over the whole interval of the particle motion. 
It should be noted that the application of the solution of a homogeneous 
Maxwell's equation instead of the inhomogeneous one essentially simplifies the 
analysis of the radiation problem and enables one to consider various cases of 
radiation emission taking into account multiple scattering. 
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3 Parametric X-ray radiation (PXR) 



Using equations (120|) -( 122|) . one can easily obtain the explicit expression for the 
radiation intensity and that for the effect of multiple scattering on the process 
under study [2|17|85] . 




Consider, for example, the PXR radiation. Let a particle moving with a uni- 
form velocity be incident on a crystal plate with the thickness L being L ^ Lc, 
where L^. = (uq)"^''^ is the coherent length of bremsstrahlung q = 9 /A and 9 
is the mean square angle of multiple scattering. The latter requirement allows 
neglecting the multiple scattering of particles by atoms. A theoretical method 
describing multiple scattering affect on the radiation process is given in ^125j . 

According to ( l22l) . in order to determine the number of quanta emitted by 
a particle passing through the crystal plate, one should first find the explicit 
expressions for the solutions As was mentioned above, the field -E^"'*'^ 

can be found from the relation E^^ = (E^^^^)* if one knows the solution 
Ej^'''^ describing the photon scattering by the crystal. 

In the case of two strong waves excited under diffraction (the so-called two- 
beam diffraction case [3]), one can obtain the following set of equations for 
determining the wave amplitudes (see |146] ): 



Here k.f = k + t, t is the reciprocal lattice vector, XO) Xt are the Fourier 
components of the crystal susceptibility. It is well known that the crystal is 





(23) 
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described by a periodic susceptibility (see, for example, [3]: 

X{r) = Xr exp(iff). (24) 

Cs = e*e^, e*(e|f) are the unit polarization vectors of the incident and diffracted 
waves, respectively. 

The condition for the linear system ( 123|) to be solvable leads to a dispersion 
equation that determines the possible wave vectors A; in a crystal. These wave 
vectors are convenient to present in the form: 

k^s = k + Kfj^gN, = ——e^g, 

where /i = 1, 2; is the unit vector of a normal to the entrance crystal surface 
which is directed into the crystal, 

1 If 2 

£^1(2). = - [(1 + f3i)xo - Attfi] ± - |[(1 - /3i)xo + Atti?] 

+4P,C'^XrX-ry'^' ■ (25) 

aB = {2kf + T'^)k~'^ is the off-Bragg parameter {as = if the exact Bragg 
condition of diffraction is fulfilled), 

7o = ■ N, = -, Pi = — , 71 = n^r ■ N, n^r 



k 7i \k + T\ 

The general solution of (123!) inside a crystal is: 
2 

= E [^A^exp{tk^sf^ + e^A.^exp(2fc^,.f)] . (26) 



By matching these solutions with the solutions of Maxwell's equations for the 
vacuum area, one can find the explicit form of ■**(r*) throughout the space. 
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It is possible to discriminate several types of diffraction geometries, namely, 
the Laue (a) and the Bragg (b) schemes are most well known. 

(a) Let us consider the PXR in the Laue case. 

In this case, the electromagnetic waves emitted by a particle in both the 
forward and the diffracted directions leave the crystal through the same surface 
{kz > 0, kz + Tz > 0), the z-axis is parallel to the normal (where is the 
normal to the crystal surface being directed inside a crystal). By matching the 
solutions of Maxwell's equations on the crystal surfaces with the help of (!23|) . 
(BSl) . (12^ . one can obtain the following expressions for the Laue case: 



e ^0 



Akr 



ik-rf 



e{-z) 



Akr 



T* —'i-—£Z,{L—z) 



+el/3i _ 

x^(L - z)e{z) + e'e'''''e{z - L), 



^ik-rT 



(27) 



where 



2s 



2s 



2£2,ls - Xo 



2(e2s — £^ls 
CsX— r 



2(£2s - £ls)' 



[^) = 1 if 2 > and e{z) = if ^ < 0. 



(2J 



Substitution of (127|) into (122|) gives for the Laue case the differential number of 
quanta of the forward directed parametric X-rays with the polarization vector 
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2 



X [e^("-^^-^^)^ - 1] , (29) 
where T = Ljc^^ is the particle time of flight; ei || [kr]] 62 \\ [kei]. 



One can see that formula fl2I?l) looks like the formula which describes the 
spectral and angular distribution of the Cherenkov and transition radiations 
in the matter with the index of refraction n^s = kz^s/kz = 1 + i^^slkz- 

The spectral angular distribution for photons in the diffraction direction kr = 
k + f can be obtained from fl29]) by a simple substitution 



Cr u XtCs 

Sl(2)s ■ 



2(£:is - e2s) 
k y kq-i k^g y krj-^g k^g -\- T. 



(b) Now let us consider PXR in the Bragg case. In this case, side by side with 
the electromagnetic wave emitted in the forward direction, the electromagnetic 
wave emitted by a charged particle in the diffracted direction and leaving 
the crystal through the surface of the particle entrance can be observed. By 
matching the solutions of Maxwell's equations on the crystal surface with the 
help of ( l23!l . (l25l) . (I26l) . one can get the formulas for the Bragg diffraction 
schemes. 

It is interesting that the spectral angular distribution for photons emitted in 
the forward direction can be obtained from (!29|) by the following substitution, 



22 



7i(2)s = 2£:2(i)s - Xo (2£:2{i)s - Xo) - (2£:i(2)s - Xo) 

X exp l — {€2(l)s -€i^2)s)L 

L 7o 



(30) 



The spectral angular distribution of photons emitted in the diffracted direction 
can be obtained from (l29l) by substitution 



Cg y CgT, k y k^j kug y k 



^/.s exp 



T 

7o 



^7, 



where 



71(2)5 = -A[CsXt] (2£2(1)s - Xo) - (2£l(2)s - Xo) 



X exp 



7o 



-1 



The angular distribution for the photons emitted at large angles in the the 
Bragg case was derived in |125] . From ( l29l) . ( I30l) we can obtain the angular 
distribution for the photons emitted in the forward direction: 



+ exp 



UJB {T^ + ^'-X^f-mrs 



{T' + ^^-Xof-mrs 



(7-' + ^2-Xo)^ 



sin^ 



27oC 



7-2 + _ 



-2 



(31) 



According to ( !3T|) . the PXR angular distribution for this case oscillates as a 



function of {}, L, ub- If "(9^ ^ 7 ^, Xo, the oscillation period is 'Oqs = \JcIujbLq. 
For kB = {(jJb/c) = 10^ cm~^, Lq = L/^q = 10^^ cm, we have 'dos = 3 x 10^^. 



For low energy electrons, oscillations in N^^ disappear. Let us note that the 



PXR photon number is proportional to Q^. As a result, the PXR intensity is 
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very high for heavy nuclei. 

For example, for Pb the photon number may be 1 per nucleus for L = 1 cm. 
It can be used for the detection of particles and for the precise measurement 
of their energy. 

4 Surface parametric X-ray (quasi-Cherenkov) radiation (SPXR) 
and DRO 

When a particle travels in a vacuum near the surface of a spatially peri- 
odic medium, new kinds of radiation arise - surface parametric (quasi- 
Cherenkov) X-ray radiation (SPXR) and surface DRO (see Figure [3]). This 
phenomenon takes place under the condition of uncoplanar surface diffrac- 
tion, first considered in [1]. 




Figure 3. 



The solution of Maxwell's equation E^^'' (r) in this case of uncoplanar surface 
diffraction was obtained in P]. It was shown that the surface diffraction in 
the two-wave case is characterized by two angles of total reflection (several 
angles in the case of multi-wave diffraction [27J). The solution obtained in [27] 
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contains the component, which describes the state that damps with growing 
distance from the surface of the medium, both within the material and in the 
vacuum, and which describes a surface wave, i.e., a wave in which the energy 
flux is directed along the boundary of the surface of a spatially periodic target 
(see review [28j). According to [4j, this solution, which describes scattering 
of a plane wave by the target under the surface diffraction geometry, can be 
written in the form: 

Ef^' = e.e*^*" + As{k, u)e'''''' + B,{k, u)e'^^\ (32) 
where the wave vector in a vacuum k = {kt,k±), ki = {kt — k±)\k2i.\ = 



k"^ — fcif, k2 = {k2t — k2±), k2t = kf + 2ttt, kt is the component of the 
wave vector that is parallel to the surface, r is the reciprocal lattice vector, 
u is the photon frequency. The amplitudes Ag and Bg are given in [5|26] . 
Substituting the solution E^^' = (E^^^')* into we can find the spectral- 
angular distribution of SPXR and DRO. For example, in the case of DRO, 
the differential number of emitted photons for a particle moving parallel to 
the crystal surface is 



d'^Ng e'^cuT ,2 



dojdVL 2t: he? 



uBs{k,uj) 6{ktU + 2'KTU-uj)e-^'^^''^^\''K (33) 



Here we assume that a particle moves parallel to the target surface at a dis- 
tance Zq at a constant velocity u; T is the flight time. The argument of 6- 
function in fl33|) is equal to zero for the frequencies 
1271™! 



ntu/c 



where fit is the component of the unit vector in the direction of fc, which is 
parallel to the surface. 
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After integrating ( 133|) over the frequencies, the angular distribution of radia- 
tion takes the form: 



-2lmfc2x|2o| 



(34) 



The spectral-angular distribution of SPXR generated by a particle incident 
on the crystal at a small angle relative to the crystal surface was obtained in 



It should be noted that, according to the analysis made in [25], the formation 
of circularly polarized quanta under the surface diffraction is possible and, as 
a consequence, such circularly polarized quanta can be produced in the SPXR 
process [32] . 

Let now an oscillator with the vibration frequency Q in the laboratory system 
move along the surface. In this case we also obtain the expression similar 
to equation fl33|) with the following replacements made in the exponent: u by 
uj±Q and v by the velocity amplitude at frequency Q. As a result, the integrals 
appearing in equation fl55]) will give 5-functions of the form 5{ktV — a; ± 
5{k2tV—uj±VL). Now all three integrals are non-zero. The first one describes the 
normal Doppler effect, the second one, the same effect including the influence 
of a mirror reflected wave on radiation. Of greatest interest is the third integral. 
In this case the 5-function leads to the equality uj = {2TTTtV±Q){l — ntv/c)~^. 
If 2'kttv > Q, then both signs of the frequency Q are allowed for radiation. 
From the quantum viewpoint, one sign corresponds to the emission of quantum 
by the oscillator (or atom) when it drops to a lower energy level, and the 
opposite sign corresponds to the inverse process: of the quantum emission 
when the oscillator (or atom) rises to a higher energy level. In other words. 
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the phenomenon of surface diffraction results in the appearance of the vacuum 
anomalous Doppler effect |25j . 



Let now a beam of charged particles (or oscillators) move along the surface of 
a natural or a photonic crystal in a vacuum. The phenomenon of spontaneous 
radiation causes the beam's instability relative to the photon emission and 
the formation of the charge density wave in the beam. The processes consid- 
ered above also leads to the appearance of such instability. According to [25] , 
multi-wave diffraction of emitted photons in this case also leads to a different 
reduction of the generation threshold similar to that appearing when a beam 
passes through a crystal |47ll48f49f50f51fl22] . 



Due to such instability the beam radiates photons collectively and its lon- 
gitudinal energy decreases. The presence of the external field (in our case - 
excitation of surface diffraction by the external field) can accelerate the beam. 
Note that the instability studied in |47f48f49f50f51fl22] is a particular case of 
instabilities caused by the processes of the emission of waves (instability may 
be caused by, for example, the parametric process of the pump wave splitting 
into two waves, the Mandelstam-Brillouin effect, four-wave processes) |159j . In 
all these cases, the power of the root dependence of the instability increment 
changes in the periodic medium if at least for one of the waves, the diffraction 
conditions are chosen according to the requirements of |122j of the coincidence 
of the roots of the dispersion equation characterizing the periodic medium. 
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4-1 Parametric X-ray radiation in crystals under action of high frequency 
ultrasonic waves 

According to [5|T7] . in the presence of an external variable field (for example, 
an ultarsonic field), a crystal is characterized by an effective index of refraction, 
depending on external field parameters. By varying these parameters, one can 
change the properties of parametric radiation. 

As the characteristics and yield of PXR depend on the solution E^~\f,uj) 
of the homogeneous Maxwell's equation describing the diffraction process in 
crystals, the investigation of the influence of an external ultrasonic (US) field 
on diffraction of X-ray points to the strong influence of a US external field 
on the PXR process. In [33] it was pointed to the essential modification of 
scattering process and X-ray radiation process under the diffraction condition 
in crystals. Due to a dynamical character of PXR formation, according to 
[3l] . the influence of a US wave on this process will be maximum when the 
US wavelength coincides with the period of extinction beatings. 

The theory of PXR under the action of an external US wave on a crystal target 
was derived in [35|36f37f38f39] . The boundary problem of diffraction of X-rays 
by a crystal target subjected to an external US wave was solved for the case 
of two- wave diffraction in P0|41f42] . Here we do not give the expressions for 
photon wave functions E^~\k,r) because they are very clumsy. The spectral- 
angular PXR distribution in the presence of a US wave was obtained in [36f37] 
(for detail see |139j ). Experimental observation of the effect of the US wave 
on parametric radiation was performed in |160fl61j . Experimentally observed 
features of the PXR in the US wave were explained in [162j . 
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5 Diffracted X-ray radiation from channeling particle (DCR) 

As we discussed above, the X-ray radiation of a relativistic oscillator in a 
crystal essentially modifies under diffraction conditions of emitted photons. A 
new diffracted radiation of oscillator (DRO) appears as a result of coherent 
summation of two processes - photon radiation and photon diffraction, but it 
cannot be reduced to a sequence of these two processes. The relativistic oscilla- 
tor itself can be a relativistic atom or a relativistic charged particle channeled 
in the potential well of averaged crystallographic potential of axes (planes), or 
an oscillator formed by an external electromagnetic field (ultrasonic, laser). It 
was shown in [T7|19|20|21|22] that the DRO spectrum is rather complex and 
is determined by the complex and anomalous Doppler effect (see Sections [H 

ED. 

It is known that the transverse energy of channeled electrons (positrons) is 
discrete and state-to-state transitions result in radiation, i.e., in this case a 
channeled particle is like a one-dimensional or two-dimensional oscillator with 
the eigenfrequency in the laboratory frame flnf = — where and 
Ef are the eigenvalues of corresponding one- or two-dimensional Schrodinger 
equation, in which the particle rest mass is replaced by the total energy m'y. 

For the analysis of DRO characteristics, it is necessary to obtain the spectral- 
angular distribution. The description of the channeled particle motion with 
the help of one- or two-dimensional Bloch functions was given in [TO]. The 
expressions for spectral-angular DRO distribution for different cases of photon 
dynamical diffraction were obtained in [TT|12f20j . For example, in the case of 
two-wave Laue diffraction the spectral-angular distribution of DRO can be 
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written in the following way [TT|12f20] : 



dudVL n'^hc^ 



nn 



nf 



-1 — IQ" fl 



5Z ''MS ^1^^ 

11=1,2 'Jznf 



(35) 



where 



COS ( 



o _ i^A 

7o 



(36) 



6' is the angle between the photon wave vector directed at a small angle 
relative to the particle velocity and the z-axis. In the dipole approximation, 
which is true for the X-ray radiation, we have 



9nf 



—I 



in an arbitrary nondipole case Qnf is defined in p!T][T2] 



Pnf= I = fnkip)f*fkip)(^^P^ 



ip^j: and (pj^j: are the two-dimensional Bloch functions satisfying the equation 
similar to the Schrodinger equation (see [TT]|12ll20j ). L is the crystal target 
length, Qnn is the population probability of the particle transverse energy 
state n. 



According to (135|) . the maximum intensity should be observed at the angles 
and frequencies that satisfy the equation 



U! 

uj{l - cos 9) - Qnf = 

7o 



(37) 



The solutions of this equation were obtained above. 
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In the case of rather thick crystals, the angular distribution of DRO was 
obtained in [T71I12II20] . For example, the angular distribution of radiation gen- 
erated by plane-channeled particles can be written as [14] : 



271 
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for r-polarization: 



Fr{0,ip) = < pioj'i^^f sm COS ip— — 



Tz sin 6 sin <f — Ty cos 9 



and for vr-polarization 



/3iu;^J sin 6 cos v9[cos 6{n, f) - T2] [sin^ 9 cos (p{nif) - r. 



n 2 



\T±\ 



where 



(^nf = ^nf{.l - P cosO - ^Re(5^s(a;^'p) \ 



is the matrix element obtained from one-dimensional Bloch functions, Lgfj 
is the effective length (at L < Labs, L^^ = L, L :^ Labs, L^k = Labs, where Labs 
is the absorption length). The term in square brackets takes account of the 
influence of the dispersion of the medium on the angular distribution. As the 
frequency, satisfying dispersion equation fl37|) . goes over from one dispersion 
branch to another with changing the radiation angle 6, the summation over 
fi means that we select the corresponding root of the dispersion equation for 
each definite radiation angle 6; ni is the unit vector directed along the wave 
vector of the photon propagating at a small angle relative to the mean velocity 
of a channeled particle. 
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Numerical calculation of the angular DCR distribution taking into account 
the dispersion characteristics of the medium under diffraction conditions, was 
made in |25 • According to Figure IU the angular distribution has a fine struc- 
ture which corresponds to the region of transition from one dispersion branch 
to another. 




Figure 4. 

One can see that the DCR distribution looks like two narrow rings - one of 
which corresponds to the solution /i = 2, the other, to /i = 1 (in the case 
of PXR generation, only one solution satisfies the Cherenkov condition). The 
angular position of the DCR distribution maxima can be estimated in the first 
approximation as: 



6*1,2 



61 ± 



l + 2sin2^B/^ 



(39) 



where 6d is, the angle that satisfies the following equation 



9' + j' + 



2Q 



nf 



B 



Ob and ujb are the angle and the frequency satisfying the exact Bragg con- 
dition. For example, according to Figure HI the values of these angles are 
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01 = 6.1 ■ 10~^ rad and 62 = 6.2 • 10~^ rad. The ratio of the angular width to 
the value of the angle 9 is about A9/9d = 0.1. The expressions for the an- 
gular distribution of radiation are simplified essentially if the particle energy 
is rather small (1 — /3 ^ l/7o,iRe5^s). In this case we can consider that the 
frequency corresponding to the maximum intensity does not depend on the 
dielectric properties of the crystal, being determined by the radiation angle 
alone. The DRO characteristics for this case were considered in [13]. 

In [H] the possibility of experimental observation of the DRO by measuring 
the angular distribution was analyzed. It was shown that for such experiments 
particle beams of high quality are required because the radiation characteris- 
tics are very sensitive to the parameters of a particle beam. Indeed, the DRO 
angular distribution shown in Figure H] takes place only for a particle beam 
whose characteristics satisfy the following inequality: 



where 7 is the Lorentz factor, A6'efr is the angular spread, A7/7 is the energy 
spread, AQ = 27rc/L is the divergence of the oscillator eigenfrequencies. 

In the opposite case, the angular width of the maxima is equal to 



As an example, the dependence of DRO angular distribution characteristics 
on the energy divergence of the particle beam is shown in Figure E] one can 
see that this dependence is rather sharp indeed. 

The dependence of the angular density of radiation on the energy of a rela- 
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tivistic oscillator, which has a "resonance" character at a given angle of ra- 
diation observation, was also considered in [H]. If the observation angle is 
equal to zero, the maximum of the angular distribution sharply increases with 
7 — )■ = {ub I'^^Y^'^ ■ When the frequency ub equals ujb = (^mdx = 2^7^, 
the maximum value of the radiation density is observed at 7 being a little 
larger than 7^?. The angular distribution in this case looks like a bell and its 
width decreases sharply with 7 — )■ 7_r (see Figure E]). 

In the range of 7 > 7/? (w^ < Wmax)? the single narrow maximum splits into 
two peaks is fixed), which shift to the range of larger radiation angles 6 
with increasing particle energy E. In [44|, the relative estimation was given 
for the contributions from different radiation mechanisms to the total radi- 
ation angular distribution which can be observed in a definite reflex. It was 
shown that at A7/7 ~ 1%, the ratio of the DRO angular density of diffracted 
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Figure 6. 

bremsstrahlung at ^ = is 

^ ^ ^DRO ^ Qnngi(sin^ + cos^ 29b cos^ V') 

where the estimation is given for a channeled electron (positron), 6^ is the 
Lindhard angle, ip is the angle between the particle oscillation plane and the 
diffraction plane, 9^ is the mean square angle of multiple scattering per unit 
length, 1/4 6'! is the classical estimation of the magnitude of \xnf\'^^'^c~'^. For a 
Si crystal and the channeled electron with the energy E = 23.6 MeV (7 ~ 7^?,; 
planes of channeling (100), ^ = 0, diffraction plane (220)) the value of the ratio 
is Ri = 25, that is, the DRO intensity is 25 times larger than the intensity 
of the diffracted bremsstrahlung at A7/7 = 1% and Aip'^'j'^ < A7/7. If the 
diffracted radiation is observed at the angle ^ 7^ 0, we should compare it with 
the contribution from the parametric (quasi-Cherenkov) radiation. In this case 
the analogous ratio is estimated as [mil39j 



35 



7 

;i + ^d7' + 7^7,^')' (sin' + cos' cos' 29 b, ^^^^ 



(sin (y9 + cos^ 9? cos^ 26'^) 

where (p is the angle between the wave vector k and the diffraction plane, 
7„ = lob/ojl is the Lorentz factor corresponding to the threshold magnitude 
of the energy E = mc'((7o)^^^'- see that this ratio essentially depends 

on the value of the azimuthal angle (p. For example, for a Si crystal this ratio 
is estimated as -R2 — 5 when the electron with the energy E = 34 MeV is 
channeled between the planes (100) and the diffraction plane is (220). 

Thus, the experimental observation of the diffracted radiation of oscillator is 
possible with the help of the particle beams of high quality. 

A relativistic oscillator can be formed not only by an unperturbed crystal 
channel but also by an external ultrasonic or laser field [T3|25fl23] . A more 
detailed treatment of particle radiation in crystals under the laser wave was 
given in [?5]. In [l5j the radiation of electrons (positrons) in a crystal subjected 
to a laser wave, which forms an oscillator, was considered. The intensity of such 
radiation was estimated. A relativistic oscillator can be a channeled particle, 
which moves in a plane channel bent by a variable external field (ultrasonic 
or laser wave), i.e., in some crystal undulator [13j. In this case, the oscillator 
frequency in the laboratory frame is Q' = — ^, where k is the wave vector 
of an external wave in a crystal, f2 is its frequency (the z-axis is chosen along 
the direction of the average particle velocity u). Diffracted radiation of the 
oscillator formed by an external ultrasonic wave was considered in 



According to [I3], the trajectory of a particle moving in the dynamic ultrasonic 
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undulator is written in the form (see Figure [7]) 
r{t) = fch(t) + f^{t) = fch(t) +dcos{n't + 5), 



(44) 



where fch{t) is the radius vector describing the motion of an ordinary high- 
frequency channeled particle, and P(t) is the radius vector describing the 
motion of a particle in the dynamic undulator. 




Figure 7. 

Assuming that the frequency of an ultrasonic wave is much smaller than the 
frequency of particle oscillations in a crystal channel, we can consider these two 
kinds of particle motion independently: the motion of an ordinary channeled 
particle and the motion of the equilibrium trajectory center of particle gravity 
inside the bent channel formed under the action of the external variable field. 
a and 6 are the amplitude and the initial phase of particle oscillation in the 
ultrasonic channel. It should be noted that if the amplitude of the ultrasonic 
wave satisfies the condition a ^ uU{EdK^)^^ {E is the particle energy, d is 
the width of the crystal channel, and V is the depth of a potential well for a 
crystal channel), then the radius of the crystal curvature due to the action of 
the ultrasonic wave is much larger than the radius of the trajectory curvature 
for the channeled particle incident on the crystal at the Lindhard angle. In this 
case, the equilibrium trajectory of a positively charged particle gravity center 
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corresponds to the trajectory of a stable channeling regime, and the curvature 
of the crystal channel caused by the action of the ultrasonic wave leads only 
to the displacement A of the equilibrium trajectory center of gravity during 
the particle passage through the crystal. That is why for positively charged 
particles, for which + A < one can take into account the dechanneling 
effect, which is due to the channel curvature, by considering the mean square 
angle of multiple scattering in this bent channel in the same way as in an 
amorphous medium [23] (aj is the amplitude of particle oscillation for the 
ordinary channeling regime). 

In the case under consideration, an essential difference arises in comparison 
with the case of diffracted radiation from the oscillator caused by a channeled 
particle. This is that the atomic (nuclear) oscillations, resulting in the for- 
mation of the ultrasonic undulator, will simultaneously lead to the dielectric 
constant modulation in a crystal and, consequently, can change the diffrac- 
tion process itself. As a result, the photon wave function changes. Maxwell's 
equations describing this situation are given in [23]. The case when the influ- 
ence of an ultrasonic wave on the X-ray diffraction process can be reduced to 
the change of the magnitude of the Fourier components of the crystal dielec- 
tric susceptibility alone was considered in detail. The spectral-angular distri- 
bution was obtained and the contributions of parametric (quasi-Cherenkov) 
radiation and of DRO itself were separated. The spectral and angular charac- 
teristics were analyzed and the total number of photons in a diffraction peak 
was estimated. 

It was shown that if the following inequality [aVt'Y > {o-f^fY is fulfilled, the 
diffracted radiation from a particle in the external field will be more intensive 
than the DRO from an ordinary channeled particle (aj and VLf are the am- 
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plitude and the frequency of particle oscillations in an unperturbed channel). 
According to [13], this inequality can be realized for a standard ultrasonic 
field source and, as shown by the estimations, the influence of this wave on 
dechanneling process can be ignored in this situation. 

In conclusion it should be noted that the diffracted radiation can also be forced 
under the motion of the oscillator over the crystal surface, by analogy with the 
surface parametric (quasi-Cherenkov) radiation [25]. Radiation from particles 
moving in crystal undulators is now being actively studied both theoretically 
and experimentally (see 



6 X-ray radiation from a spatially modulated relativistic beam in 
a crystal ( fL47\ ) 



As was mentioned in the Introduction, in our works we suggested a new type of 
the free electron laser based on volume multi-wave distributed feedback. It was 
shown that multi-wave distributed feedback enables reducing the generation 
threshold from the values of the current densities of a relativistic beam, such 
as 10^^ A/cm^, which are practically unavailable, to acceptable values, such 
as lO'^ — 10^ A/cm^. This reduction can be achieved even in the X-ray range 
(10-100 keV). 

One of the most important ways of reducing the generation threshold in free 
electron lasers consists in prior modulation of the density of a beam of rela- 
tivistic particles. In particular, modulation gives rise to radiation coherent to 
the beam, which appears alongside with spontaneous radiation of electrons as 
the beam gets into the undulator. 
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According to |25f49f51f59fll9fl37j . one of the schemes of the X-ray laser is 
based on using parametric quasi- Cherenkov radiation, diffracted radiation of 
a relativistic oscillator in the crystal (or radiation in a crystal undulator). In 
this case the relativistic electrons (positrons) move not in the vacuum, but in 
the crystal. Interaction of particles with the atoms of matter result in multiple 
scattering, which will apparently diminish the beam's modulation amplitude 
and, hence, the intensity of coherent radiation. 

In |147] . the spectral angular distribution of the intensity of parametric X- 
ray (quasi-Cherenkov) radiation was obtained as well as the radiation of a 
relativistic oscillator formed by a spatially modulated beam under multiple 
scattering. The conditions were determined under which multiple scattering 
does not affect the intensity of coherent radiation. It was shown that the 
number of quanta coherently emitted by the beam in the X-ray range with 
reasonable requirements to the value of the current density of the beam of 
relativistic particles, such as ~ 10® A/cm^, is too small in comparison with 
the number of spontaneously emitted quanta. For this reason, the experimental 
observation of the effect of coherent radiation of quanta in the X-ray range 
should be carried out in a rather narrow spectral angular range. In this case, 
in the wavelength range of 50 - 100 A, it is possible to observe the effect 
when the degree of beam modulation is as small as /i ^ 10^^. 

In our case of relativistic particles, their energy is much greater than the en- 
ergy of emitted 7-quanta. So one can consider the radiation process within 
the framework of classical electrodynamics. Let a beam of charged particles 
(electrons or positrons) traverse the area occupied by matter. The spectral 
density W^^ of the radiation energy per unit solid angle of photons character- 
ized by the polarization vector {n = k/k, k is the photon wave vector) can 
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be found when one knows the Fourier transform in time of the electric field 
strength in the produced electromagnetic wave E{f,u). According to [7ll5] . at 
long distances from the target 



where E^^^{f,oo) is the solution of the homogeneous Maxwell's equations, 
which describes scattering by this wave of the photon with wave vector k and 
polarization s. The scattering asymptotics has a form of an incident plane 
wave plus a converging spherical wave. In this case Et^-> = E'}t\ where 
E is the ordinary solution of the homogeneous Maxwell's equations with 
asymptotic behavior containing a diverging spherical wave; 



is the Fourier transform of the beam's current, which equals 
]{f,uj) = eY,Vi{t)6{f-r{t)), 

i 

ri{t) is the coordinate of the z-th particle of the beam at moment t, Vi{t) is 
the particle velocity; e is the particle's electric charge. Using fHS]) . we obtain 
the following expression for spectral-angular energy distribution of radiation 
generated by the particle beam in the target: 



j Ei {f,uj)j{f,uj)d^r 



(46) 



The vinculum in (146|) denotes averaging over the distribution of the coordi- 
nates and velocities of the particles in the beam (allowing for multiple scatter- 
ing in the target). With the help of the distribution function w{fi, vi, t; r^, v'^, t') 
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defining the joint probability density of finding in the l-th particle the coor- 
dinate fi and velocity vi at moment t, and finding in the m-th particle the 
coordinate and velocity -u^ at moment t', (H^ can be written as follows: 



According to P7|) . the spectral-angular distribution of energy Wi^ can be 
represented as follows: 



I l^m 

where 



has a meaning of the amplitude of photon emission by the /-electron; = v/c. 

In the case when correlations in the positions of different particles can be ne- 
glected, w{fi,vi,t; fm,Vm,t') = Wi{fi,Vi,t)wi{f.m,Vm,t') (/ ^ m) . As a result, 
( 148|) can be recast as 



WL = Y.i\^^l?-\^^l?) + 



I 



(49) 



The first sum in fl49|) describes spontaneous incoherent radiation of photons, 
the second one - which is proportional to the squared modulus of the sum 
of averaged emission amplitudes - describes coherent across the beam photon 
emission. To define the conditions under which the distribution of coordinates 
and velocities in the beam has no effect on the process of coherent radia- 
tion of a photon, let us take into account that when a photon is emitted 
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in a homogeneous medium (e.g., due to the Cherenkov effect), Et^~^ has a 
form of a plane wave, while in the case when parametric X-ray radiation or 
diffracted radiation of the oscillator is generated in the crystal (in a medium 
with spatially-periodic dielectric permittivity), the expression for Et^~^ has a 
form of a superposition of plane waves [5|17f20f85j . In particular, in the space 
region occupied by the target, E^ can be represented in the form: 

4(-)(n(t),a;)=E^ne*-^'(*), 

n 

fi{t) = fio + uit + 6fi{t). (50) 

Here fiQ is the electron coordinate at time t = 0; ui is the electron velocity 
in a vacuum; 6fi{t) describes the variation of the electron trajectory under 
the action of the forces (in particular, those leading to multiple scattering) 
that affect the particle in the area occupied by the target. Let us assume that 
Ui = u + Aui, where Aui ^ u. To be more specific, take the 2;- axis as directed 
along u. 

Now, let us determine the conditions under which the distribution of the co- 
ordinates and velocities in the beam has no effect on the intensity of coherent 
radiation. 

Let the crystal surface through which the electrons enter the target be located 
at point z = 0. At the moment of time t = 0, their coordinates were tqi 
and were located in the area z < 0. Consequently, the Z-electron reaches the 
surface z = at time ti = \ziq\/uzi. During this time, the electron will move 
over the distance 6ri± = ui±\ziq\/uzi in the transverse direction. As a result, 
the electron's transverse coordinate at entering the target will be 

^ , ^ \zio\ ^ ^ ZiQ 

ri± = rio± + ui± = rio± - Uij_ — . 

Uzl Uzl 
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The major contribution to the amphtude fii comes from integration over the 
interval corresponding to particle motion inside the target [8]. We shall con- 
sider this very contribution. In the integral over time involved in the amplitude 
fii, we shall shift the zero-time position corresponding to the particle motion 
in the target area to point t = 0. 

As a result, the amplitudes fii can be represented in the form: 

= A;/-*^^^^'o^e"'^''V'^"'^^, (51) 



where Ai is the amplitude of radiation in the crystal of the Z-electron, which 
entered the target at point r = at time t = 0. If the additional phase shift 
occurring during the time T of the electron motion in the crystal is small due 
to velocity distribution Aui, i.e., kAuiT < 1, then multiple scattering has no 
effect on the amplitude of radiation. Here, in the case of parametric (quasi- 
Cherenkov) radiation, all the amplitudes Ai equal one another, while in the 
case of diffracted radiation of the oscillator (radiation of channeled particles), 
the amplitudes Ai depend on the the particle entry point into the channel, 
which should be taken into account in averaging the amplitudes. It is known 
that the root mean square angle of multiple scattering (-(9^) = g^j, where 
Eg =21 MeV, E is the particle energy, / is the traveled path length, L is 
the radiation length, g is the coefficient for the difference between (i?^) and 
('(9^)am in an amorphous medium. For positrons moving in the regime of planar 
channeling, the magnitude of g can be much less than unity. When move at 
small angles relative to the axes, the magnitude of g can become much greater 
than unity. Two requirements follows from inequality kAuiT < 1: k±^Jl^^d^, 
I < 1 and kli^d"^) < 1. The former can always be chosen with the help of the 
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observation angle. The latter yields the requirement / < where is 

the root mean square angle of multiple scattering over the unit length (Z ~ 7, 
7 is the particle Lorentz factor). This leads to the fact that even for quite 
hard radiation with the quantum energy of the order of 100 keV, the target 
thickness required to avoid the effect of multiple scattering on the radiation 
process, appears rather large (/ < 10^^ cm) for electrons with the energy of 
1 GeV. We shall further assume that this condition is fulfilled. So from ( H9|) 
we have the following expression for the spectral angular distribution of the 
number of emitted quanta cPN/dudQ = W,^^/hu: 

2 AT A2 ] 



dudVt dudVt dudVt 



where 



ri = ro±i - ui± — , — zoi , K = [k±,- 

Uzl Uzl / \ u 



in the case of parametric radiation; if radiation is produced by a relativistic 
oscillator (radiation from channeled particles, diffracted radiation of the oscil- 
lator) with the vibration frequency VL in the laboratory system of coordinates, 
vector K = (k±, • From flS^ follows that to achieve the conditions under 
which the velocity distribution Aui in a beam has no effect on ( l52l) . the lon- 
gitudinal dimensions Lf, of the bunch of particles should satisfy the relations: 

Lb < 1 and k± < 1, 

u u u 



I.e. 



Lb < 7-^ and Lb < 
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where v is the characteristic angular distribution of particles' velocities in the 
beam; is the quantum emission angle. When the longitudinal dimensions of 
the bunch satisfy these conditions, the magnitude of vi in (132]) can be taken 
equal to roz- Now averaging in f l52|) is easy to perform. Since the double sum 
in ( 152|) does not depend on velocities of particles, averaging is reduced to 
averaging over the initial distribution of the coordinates of particles in the 
beam. Let us introduce the beam density p(r); Jv^p{r)d^r = (H is the 
bunch volume, N^. is the number of particles in the bunch). Upon averaging 
with the distribution p(r), we obtain 



dudil dudVt dudVt 



2 



(53) 



where d'^Ni/dcudfl is the spectral angular distribution of quanta formed as a 
result of spontaneous radiation by a single electron. Let us consider a beam 
modulated as p(r) = po + Picos(rr). As a result, (154|) will include rapidly 
oscillating integrals, which can be replaced by 5-functions with good accuracy. 
This enables writing fl53l) in the form 

+ -r-^-N,i2^po6{K - r), (54) 



dujdfl dudfl dudVt 2 

where p = pi/po- From (!54|) follows that the total number of quanta emitted 
by the bunch equals 

d^N^ 



N = NANi + c 



dudVt 



K=f ^ '^O 



Here N\ is the number of quanta of incoherent spontaneous radiation, which 
are produced by a single electron traversing the target under study; A^o = |^| = 
1^, Aq is the radiation wavelength equal to the spatial modulation period d of 
the beam. 
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Note that (155|) is easy to obtain from the following quantitative considerations: 
In performing integration over dudfl in (|53l) . the characteristic range of values, 
where the second term in fl53p is nonzero is: Akj_ ~ 1/Lj_h for the transverse 
dimension and Au/c ~ ^/Lh, for the longitudinal dimension. As a result, we 
have 



d'^N, r ^ 2 ^2 



dudQ 
d^Ni 



dujdVt ' ^ 



dudVt 



, , A?k\Ak\\ 
NIp'c- 



dujdVt 



k'^ 

k"^ ' 
K=T '^0 



(56) 



L\i, Lb fco dojdil 

i.e., the expression appearing in f l55|) . The ratio of the total number of coher- 
ently emitted quanta to that of incoherently emitted quanta is 

(57) 



Here Ak{AQ) is the characteristic range k (of solid angles), where the quantum 
is emitted through incoherent spontaneous radiation; Ak/k ~ I/7, AQ ~ 
1/7^; 7 is the particle Lorentz factor. In obtaining ( 157|) . the estimate A^^i ^ 
^^^^A/cAf2 was used; can also be recast as 



dkdQ 

,2 



A^coh At P0.3 



A^ incoh ^0 



7^ (5^ 



Note that according to [5], in the case of parametric radiation when 1/7^ < 
^eff — 1 (^eff is the effective index of refraction in the radiation area), the char- 
acteristic range of radiation angles (the frequency distribution) A??, Au/u ~ 
y/uefi — 1. As a result, with growing 7, the ratio ( l58l) tends to the limit 
A^coh/A^incoh — (At^Po/^o)('"'cff-l)~^^^- Recall that in the X-ray range of 10^100 
keV, neff — 1 ~ 10^^ ^ 10^^. From (1571) . (!58|) follows that with other condi- 
tions being equal, the ratio A'coh/A^incoh decreases rapidly with growing k (with 
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the decrease in the radiation wavelength). Thus for example, even at current 
densities in a beam as high as 10^ A/cm^, which corresponds to po — 10^'' 
and at = 1, the ratio A^coh/A'^incoh < 10^^ for quantum energy of 10 keV 
and 7 = 10^. At the same time, the ratio of the spectral densities of radia- 
tion within the range of the emission angles At? ~ l/ZcL^i, ^ I/7 of coherent 
radiation is fi'^N^. 

Thus, within this range of emission angles, the intensity of coherent sponta- 
neous radiation exceeds that of incoherent spontaneous radiation when the 
modulation depth is /x > ly/T^. For example, when the number of electron 
in the bunch is A"e — 10^^, it is sufficient that fi > 10^^, which enables one 
to appreciably simplify the problem of experimental observation of coherent 
radiation in the X-ray range of 10-100 keV. Note that the formulas derived 
here are also applicable when the target is irradiated by the laser pulse of 
length greater than Lb. Indeed, let us divide the whole pulse length into the 
bunches, each of length Lb, i.e., we have ( 15^ . ( I55ll . where A'e stands for the 
number of electrons in the entire pulse. 

Now let us estimate the number of quanta which can be produced coherently, 
e.g., in the case of parametric (quasi-Cherenkov) radiation mechanism. From 
general formulas (!52|) for parametric radiation in the Laue case follows the 
below expression for spectral-angular distribution of quanta emitted in the 
direction of diffraction: 



dkdVt 



e 



2 



K=T 



Xr 



B 



2 



koN\ (59) 



where L is the crystal thickness; Xn is the Fourier component of the crystal 
susceptibility for quantum diffraction by the system of planes defined by the 
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crystal reciprocal lattice vector ri, is the Bragg diffraction angle. 

According to ( 159|) . ( 155|) . for the number of coherently emitted quanta per 
electron we have 



/Icoh 1^ 10 ^ 



B 



2P0L 



For radiation generated in S'i, Xr ~ 10 ^ for (400) plane, ko = 10^, 'ds — 45°, 
po — 10^'^, L = 10~^ cm, we have A^'icoh — 10"^^- 

In the case of PXR generation in a layered medium (i.e. in NiC) in the range 
of wavelengths ^ 50 ^ 100 A, Xt ~ 10^^, = 10^, L ~ 10"'' cm, we have 
A^ieoh ^ 10" V^- At the same time, in the case of generation of surface PXR, 
when the quantum absorption length in the medium does not restrain the 
radiation intensity, which is proportional to the path length traveled over the 
grating [2^, A'^icoh — 10^/i^ for L = 1 cm. When PXR is generated in the 
optical range fco = 10^, Xt = 10"^, L ~ 10"^ cm, we have A^'icoh — lO^p^. Note 
here that the obtained estimates increase by two-three orders of magnitude in 
the case of Bragg diffraction because of the increase in d'^Ni in narrow spectral 
ranges |138] . 

For the oscillator mechanism of radiation [20] 

□ 2 , \ 2 



dkdVt hen \ c J 



koL\ (61) 



f_L is the velocity amplitude of the transverse vibrations of the particle. For 
this reason, to obtain the estimates, \XtMb\'^ in (I6OII is replaced by (f^/c). In 
particular, for positrons channeled in Si, the characteristic values of v±/c ^ 
10"^ and of rfA^iosc ^ ^A^ipxr- 
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Thus, it follows from the above analysis that using a spatially modulated beam, 
one can observe coherent parametric X-ray radiation. Moreover, in a soft X-ray 
range, under the condition when the surface parametric radiation is generated 
[25], coherent radiation can be observed even when the degree of beam modu- 
lation is not very high (e.g., for ~ 10^^, we have A^'icoh — 10^''). As a result, 
when the transmitted pulse contains A^e — 10^^, we obtain ~ 10^, which 
is quite acceptable. The possibility to use modulated beams for generating 
coherent radiation in crystal undulators has recently been considered in ^148j . 

7 Crystal X-ray Free Electron Lasers on the basis of PXR and 
DRO (DCR) 

High spectral and angular densities of parametric (quasi-Cherenkov) and diffracted 
radiation of the oscillator as well as narrow spectral and angular widths of radi- 
ation reflex give a basis for application of considered spontaneous mechanisms 
of X-ray radiation for the construction of an X-ray coherent radiation source 
by using beams of relativistic particles in crystals. Such a system can be con- 
sidered as a crystal X-ray free electron laser (FEL). The idea of X-ray FELs 
based on spontaneous parametric and DRO (DCR) radiation in crystals was 
first expressed in [lS|47f48] . In |15|47f48f49f5Uj . the dispersion equation for the 
eigenstates of the system consisting of electromagnetic radiation, a beam of 
relativistic oscillators and a crystal was obtained. The increment of the beam 
instability was also analyzed. The possibility in principle of obtaining X-ray 
coherent radiation with the help of a beam of relativistic oscillators in crys- 
tals was shown. In |51j, the parametric (PXR) relativistic beam instability 
in a crystal was considered and the corresponding increment was obtained. 
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Radiative instability caused by spontaneous radiation in crystal undulators 
and by a laser wave propagating through a crystal was studied in [60]. Thus, 
in |46f47f48f49f5Uf51f(jU] . a new kind of the X-ray FEL - the solid X-ray free 
electron laser (SXFEL) was suggested. As we have pointed out above, sev- 
eral mechanisms of spontaneous X-ray radiation generated by a relativistic 
electron beam in crystals can constitute the basis for such SXFELs: paramet- 
ric ( quasi- Cherenkov) X-ray radiation and diffracted radiation of oscillators 
formed in crystals, for example, by channeling [I6|47f48f49f50f56f57f58f59] or 
under the action of an external field [SSpU] . 

The main feature of such an X-ray generator is that the crystal target, in this 
case, not only forms the mechanism of spontaneous radiation, but also acts 
as a three-dimensional resonator for X-ray radiation which produces a dis- 
tributed feedback (DFB). The construction of the X-ray generator by using 
channeled electron beams in crystals was also considered in |611I62II63II64II65] 
and the construction of the X-ray generator on the basis of resonance tran- 
sition radiation was discussed in [66|67] . The possibility of using the crystal 
as a resonator that produces a one-dimensional distributed feedback for an 
X-ray coherent generator was first expressed in [68]. This idea was used for 
the formation of a one-dimensional DFB in a solid X-ray FEL on the basis 
of channeled particles in [SlfSS] . However, in all these works, the DFB was 
traditionally considered in a one-dimensional geometry when the radiated and 
diffracted waves propagate along one line in opposite directions. The authors 
of [611165] obtained a low generation threshold for such a system with a one- 
dimensional DFB only due to ignoring the radiation self-absorption inside the 
crystal. The correct consideration of absorption, as it was shown in [59] , leads 
to the threshold beam density of the order of j*^ ~ 10^^ A/cm^ for this DFB 
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geometry. 



In the solid X-ray free electron laser, suggested 

a three-dimensional DFB that allows one to optimize the system and to essen- 
tially decrease the generation threshold. The analysis showed that the process 
of amplification and generation in such a crystal (natural or photonic) solid 
resonator essentially modifies and, under definite conditions, develops more 
intensively. It was shown that the interaction between the particle beam and 
the electromagnetic field is the strongest near the region of degeneration of 
roots of diffraction dispersion equation, particularly, in the case of multi-wave 
diffraction. 

Let us consider in detail PXR and DRO (DCR) crystal X-ray FELs 
8 Parametric (Quasi-Cherenkov) X-ray FEL 

The quasi- Cherenkov instability of a relativistic electron (positron) beam in a 
three-dimensional periodic medium in the X-ray range was first considered in 
[5T|139] . The authors formulated the problem of X-ray parametric radiation 
amplification in an infinite medium caused by quasi- Cherenkov instability of a 
beam of relativistic particles. The dispersion equation for the case of two-beam 
diffraction and the increment of instability was obtained. It was shown that the 
strongest interaction between the particle beam and radiation was close to the 
region of degeneration of the dispersion equation roots. The boundary problem 
of amplification of radiation in a finite parallel-plane crystal target was solved, 
and the generation threshold for the particle density was obtained. It was 
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assumed that a relativistic particle beam with a mean velocity u was incident 
at a definite angle V'o on the parallel-plane crystal target with the length L. The 
orientation of the particle beam relative to crystallographic planes was made in 
such a way that spontaneous photons radiated by a particle beam were under 
diffraction conditions for planes with low indices. The fulfillment of diffraction 
condition not only brings about the possibility of quasi-Cherenkov radiation 
in the X-ray range itself but also produces a three-dimensional distributed 
feedback. 

The closed set of equations describing the interaction of a radiating beam 
with a crystal, in the general case, consists of Maxwell's equations for the 
electromagnetic field and the equation for particle motion in the field (for 
a "cold" particle beam, 0%l) < {kL)~^, where 9 is the radiation angle, ip is 
the angular spread of particles in a beam, k is the photon vector) or the 
equation for the distribution function (in the case of a "hot" particle beam). 
For example, in the case of a "cold" beam we have: 



V X V X E{r,uj) = — —j{r]Lj) + —D{r,uj), 



at m7 [ 

Va{t) ( Va{t) 



C 

Va{t) 



H{fa{t),t) 



(62) 



where j{f,t) = eJ2va{t)d{f — fa{t)) is the microscopic current density of 

a 

particles in a beam, n{f,t) — eJ2^{'^~'^a{'t)) is the corresponding charge 

a 

density, fa{t) and Va{t) are the trajectory and the velocity of the a-th particle 

— * — * 

in a beam, E{f, t) and H{f, t) are the electric and magnetic strengths of the 
field, D{r,u) — e{r,u)E{r,u), e{n,u;) — J2 ^t{^)^~^'^^ is the crystal dielectric 
constant, Eq = l + go = 1 — cjI/cu^, cu^ = Ane'^no/nie, uq is the electron density 
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in a crystal, Qr = Sr is the Fourier component of the dielectric constant, r is 
the reciprocal lattice vector. 

For the case of two-wave generation, with the trajectory and velocity of a 
particle represented as ra(t) = foa + ut + 6fa(t) and Va(t) = u + 6va(t), 
where toq, is the position of the a-th particle in a beam at the moment of 
intersection of the crystal boundary, the system fl62l) can be written as a 
system of Maxwell's equations for electromagnetic fields E{k, to) and Er{k, co) 
in the following way [5^53f54|[55] 




equations fl63l) is written for a-polarization of radiation because it is excited 
with maximum probability at parametric ( quasi- Cherenkov) radiation, a)£ = 
47re^no/me, no is the mean density of the unperturbed particle beam. 

Comparison of ( l63l) and the ordinary set of Maxwell's equations describing X- 
ray dynamical diffraction in a crystal allows one to conclude that the boundary 
problem of X-ray amplification (generation) under penetration of a particle 
beam through a periodic medium can be reduced to the problem of X-ray 
diffraction by an "active" periodic medium, which consists of the crystal + 
radiating particle beam and is characterized by the following dielectric con- 
stant: 





(63) 



where k^ = k + f, = E{k,uj)e^, E^ = E{kT-,u)eo-, 



[kr]. The set of 
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As the electron density in a beam is much smaller than that in a crystal, the 
second term on the right-hand side of (IMl) can be neglected. The last term has 
a resonance behavior under the fulfillment of synchronism condition between 
the particle beam and the electromagnetic field u — ku = 0. In the X-ray 
range, the fulfillment of this condition for a diffracted wave is impossible, that 
is why it is possible to consider eo{kr,uj) = Sq. 

In the case of a "hot" beam, the dielectric constant of such an "active" medium 
is represented as: 



'^ij^ {ipi cos V9 + '?/'2 sin (y9 + '?/'|| / 726*)^ ' 

where x = u — ku, ipi, 1^2, "^y are the transverse and longitudinal divergences 
of the particle velocities in a beam, if is the azimuthal angle of a photon, 6 is 
the angle between the photon wave vector and the z-axis, directed inside the 
crystal clS cL normal to the crystal surface. 

Thus, the reduction of the analysis of amplification and generation processes 
in an X-ray FEL to the solution of a boundary problem of X-ray dynamical 
diffraction by an " active" crystal target of length L enables one to find gener- 
ation thresholds of such a system for different regimes of FEL operation and 
to perform the optimization of parameters [511155] . 



The dispersion relation determining the solutions for the electromagnetic wave 
vector inside the "active" medium in the case of two-wave generation is as 
follows: 
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= -5^i^(A;V - co'Kky - uh,). (65) 

and the general solution of the set of equations fl63l) can be represented in the 
form 

E = j2 e.C^S^\\ + S^e^'^y'-^ (66) 



where 



are the roots of dispersion equation (165|) . The wave amplitudes are found 
from the matching conditions for the electromagnetic field (1661) on the bound- 
aries of the crystal target and are given in 



As a result, the generation conditions for the case of a two- wave distributed 
feedback and different cases of "cold" and "hot" beams as well as for different 
regimes (weak and high gain regimes) were obtained. It was shown that in 
all cases the generation threshold has a simple meaning: on the left-hand side 
of equality there is always the term describing the "production" of radiation 
inside the crystal; on the right-hand side, there are two terms describing the 
radiation losses. Particularly, the first of these terms refers to the radiation 
losses through the crystal boundaries, while the second term corresponds to 
self-absorption of radiation inside the crystal. 

The analysis showed that the conditions of generation threshold are optimal 
near the region of degeneration of roots of dispersion equation (!65|) . This region 
corresponds to the edge of a nontransparency region in the dynamical diffrac- 
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tion theory, and the interaction of the electromagnetic field with a particle 
beam and with a crystal is most effective here. 

It should be noted that the condition of the degeneration of dispersion equa- 
tion roots leads to the requirement for the photon radiation angle 6"^ — 
(— /3i)^^/^|(7t-| — Ififol — 7^, and this, in turn, gives the restriction for the possible 
geometry of Bragg diffraction in the X-ray range — < /3i < 0, where 

r = grQ-T- 

The estimations of the threshold magnitude of the beam current density 
showed that the case of a two- wave solid distributed feedback is not an optimal 
case for achieving the generation regime in the X-ray range. If for a "cold" 
particle beam in LiH and il)±_ < 10~^ rad, ipw < 10~^ rad, the threshold 
current is of the order of 10^ A/cm^ at i ~ 0.1 cm, then multiple scattering 
of electrons makes the beam "hot" and leads to the increase in the threshold 
density of a particle beam up to jth > 10^° A/cm^. 

Because of the destructive influence of multiple scattering process on the qual- 
ity of a beam of relativistic particles and, consequently, on the threshold con- 
ditions, the crystal target length should be to reduced as much as possible. 
As is known from the optical laser theory, the mirror resonator, similar to a 
Fabry- Perot resonator, is used for this purpose. In the X-ray range, mirrors can 
be replaced by crystal plates that are oriented in such a way that the radiation 
wave vector is under the Bragg condition. Due to radiation, the generation pro- 
cess takes place in a narrow angular and spectral interval Auj/uj ~ < 10~^; 
high effectiveness of the radiation reflection under diffraction conditions by an 
external crystal resonator can be obtained under the corresponding coordina- 
tion of the resonator and the crystal target ("active" medium). This allows 
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one to essentially reduce the radiation losses through target boundaries. In 
|54j . the generation threshold for the system with the external Braggmirrors 
was derived, and it was shown that the term connected with losses through 
the boundaries could be reduced by a factor of (1 — |-R|), where R is the re- 
flection coefficient of a Bragg mirror. As a result, we can shorten the crystal 
length to the size necessary for achieving the generation threshold. However, 
the estimation showed that the threshold magnitude holds rather high. 

As the analysis in showed, the transition to the distributed feedback under 
the surface uncoplanar diffraction when the radiating particle beam is incident 
on a crystal at a small angle ip ~ relative to the crystal surface (see 
Figure E]), allows one to step down the generation threshold. First of all, the 
destructive influence of multiple scattering on the particle beam is suppressed. 
Besides, the behavior of dispersion equation roots changes, which modifles the 
process of radiation amplification. 

The disadvantage of the case of two-wave diffraction distributed feedback is 
that the coordination between the degeneration condition of dispersion equa- 
tion roots and the requirement of Cherenkov synchronism hardly fixes the ge- 
ometry of distributed feedback and leads to small magnitudes of the diffraction 
asymmetry factor This, in turn, leads to the enhancement of self-absorption 
of radiation inside the crystal target. In [28] was pointed out that the transition 
to the multi-wave diffraction allows one to modify the functional dependence 
of the increment of the particle beam instability and, consequently, to step 
down the threshold density of a beam as well. The dispersion equation for the 
three- wave coplanar diffraction geometry of distributed feedback was obtained, 
and the rule for writing the dispersion equation for an arbitrary multi-wave 
diffraction distributed feedback was formulated. In [54j, the expression for the 
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generation threshold in the case of three- wave coplanar diffraction was also de- 
rived. It was shown that in this case the Cherenkov condition was fulfilled for 
two dispersion branches, which provided the possibility of the coincidence of 
diffraction roots with Cherenkov synchronism condition near the exact Bragg 
condition and, consequently, the possibility of optimization of the threshold 
magnitude. In the case of the Laue-Bragg diffraction geometry, the threshold 
density of the beam can be reduced to jth ~ 10*^ A/ cm^, at ~ 5 ■ 10~^ rad in 
the vicinity of the double degeneration of dispersion equation roots. It should 
be noted that even in the case of three-wave generation, it becomes possible 
to apply the phenomenon of anomalous X-ray penetration under diffraction 
condition and, as a result, to step down self- absorption of radiation inside the 
crystal. 

Thus, we can conclude that the most suitable geometry for the achievement 
of the generation regime of quasi- Cherenkov X-ray radiation with the help of 
relativistic electron (positron) beams in crystals is the grazing geometry of the 
particle beam incidence on a target with the distributed feedback formed by 
multi-wave surface diffraction. 

The spectral-angular distribution of the coherent PXR near the generation 
threshold was obtained in [55] . 

In [52], the spectral-angular distribution of coherent radiation far from the 
generation threshold was derived within the framework of the perturbation 
theory, and the possibility of experimental observation of the coherent PXR in 
existing accelerators was analyzed. It was shown that observation of coherent 
parametric (quasi-Cherenkov) radiation far from the generation threshold was 
a very complicated problem for the X-ray range, but it was possible to observe 
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the coherent radiation in an optical range nowadays. 

9 The X-ray generator on the basis of diffracted channehng radia- 
tion (DCR) 

The second type of a crystal generator is based on the application of diffracted 
radiation by the oscillator (DRO) as a spontaneous radiation mechanism 
[5^57|58|59|60] . As stated above, the radiating oscillator can be formed in dif- 
ferent ways. This can be electrons channeled in an averaged crystallographic 
potential of planes or axes, or electrons moving in an electrostatic wiggler [SH] 
or, for example, an oscillator formed by an external ultrasonic (optical) wave 
in a crystal [60j . It is obvious that the general approach to the consideration of 
the generation problem with the help of a relativistic oscillator beam does not 
depend on the formation mechanism of the oscillator itself. Since the oscilla- 
tor is a quantum system, it is more accurate to perform the calculation of the 
polarizability tensor of a particle beam within the framework of quantum elec- 
trodynamics. Reduction of the problem of radiation application (generation) 
by a particle beam in a finite crystal target to the problem of diffraction of 
X-rays by an " active" medium, consisting of a crystal and a beam of radiating 
oscillator, holds true in this case as well. 

The expression for the polarizability of such an "active" medium in the case 
of channeled particles in unperturbed averaged crystal potential was obtained 



in [57]: 




2 



(67) 
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where d2i{k) is the matrix element of the operator aexp{ikr), which in the 
dipole approximation takes the form 

The axis is chosen so that it hes along the transverse particle oscillations in a 
channel, u\\ is the longitudinal velocity parallel to the channeling planes, {k^ = 
kfix, uiifix = 0), ^21 is the frequency of the transition, \\ [fcf], Wi and W2 
are the population of the states 1 and 2, F is the phenomenological constant, 
taking into account inelastic oscillations; its order of magnitude estimate gives 
(Ld)^^, where is the dechanneling length. In obtaining ( 1671) . it was taken 
into account that the synchronism condition could be fulfilled only for the 
wave propagating at a small angle relative to the longitudinal velocity of the 
particle. The fulfillment of the synchronism condition for the diffracted wave 
is impossible in the X-ray range. As the analysis showed [57], although there 
were a lot of zones (states) of transverse energy of a channeled particle, the 
main contribution to the polarizability tensor was made by a certain transition 
with the frequency ^21- This means that the consideration is reduced to the 
two-level problem. 

Indeed, the contribution to the beam polarizability from the transition be- 
tween the levels m and n is determined by the deviation from the exact syn- 
chronism condition of the radiation field with the oscillator, i.e.. Re (a; — uk — 
^mn) = 0. This contribution should be taken into account only if 



Re (a; — uk — fimn) < 



I'Oiiku — uj) + Vmn (6J 



If the magnitudes of Rew and the angle between u and k are fixed, the number 
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of transitions contributing to the polarizability depends on the relationship 



between AQ and lm{ku — u) + Tmn , where AQ determines the typical value 
of the difference Qn+i,n ~ ^n,n-i which characterizes the unharmonism of the 
averaged potential (for a harmonic potential AQ = 0). The analysis of the 
magnitude of AQ for different kinds of averaged crystallographic plane poten- 
tials shows that 



An > 



lm{ku — w) + 



and, consequently, the synchronism condition can be fulfilled only for a cer- 
tain transition Q2i- Other terms in a polarizability tensor can be neglected as 
nonresonant. It was shown that the most effective interaction between the os- 
cillator beam and the radiated wave takes place near the degeneration region 
of roots of the dispersion equation determining the eigenstates of the field in 
an "active" medium. But, as contrasted to the parametric (quasi-Cherenkov) 
generator, for which the radiation condition is realized only at large deviation 
from the exact Bragg condition, now there is a possibility to overlap the syn- 
chronism condition with the exact diffraction condition. As a result, in the case 
under consideration the manifestation of the effect of anomalous X-ray pen- 
etration through the resonator under dynamic diffraction (Borman effect) is 
possible. This circumstance is very important because of strong absorption of 
X-rays inside a crystal target. In [57] the boundary problem of X-ray diffrac- 
tion by an "active" medium of a finite size was solved, and the generation 
condition was obtained. It was shown that the beam can be in synchronism 
with one of the modes of the "active" medium. These modes correspond to 
the waves with the wave vectors being the roots {6i and 62) of the dispersion 
equation. According to [57J, the generation condition can be realized in two 
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cases: for the wave corresponding to the root 82 at the positive magnitude of 
a = a+ and for the wave corresponding to the root 61 for the negative devia- 
tion from the exact Bragg condition a = a_. It was shown that the solutions 
of the generation equation for different modes are identical in structure. All 
of them lead to the phase condition 

i(»)-4») = ??^ (69) 

where Sf''* and 62^^ are the solutions of the diffraction dispersion equation. The 
conditions of generation are written for the case of a channeled particle in [57] 
and for the case of electrostatic and magnetostatic wiggler in [58j . 



If the condition fl69|) is fulfilled, the longitudinal structure of the modes turns 
out to be close to the structure of a standing wave. That is, l-Ep and l-ErP are 
proportional to ~ cos^ ~ L). This condition is similar to a well known 

phase condition of the stand wave appearance in a mirror resonator of an ordi- 
nary laser [60]. The meaning of amplitude conditions is the same as in the case 
of the quasi- Cher enkov X-ray generator. The field amplification, due to the 
radiation process, should be equal to the radiation losses caused by absorption 
inside the crystal and the output of radiation through the boundaries of the 
crystal target. Because the gain in the weak-gain regime is proportional to the 
current density of a beam, the formula for the threshold gain sets requirements 
for the current density. The invariant characteristics of a particle beam are of- 
ten used instead of the current density, that is, the current /, the normalized 
emittance e„ = 7?" (■?/') and the normalized brightness i?„ = J/vr^e:^, where r is 
the beam radius, (ip) is the angular spread. The angular spread corresponds to 
the divergence of the longitudinal velocity as = u{ip^)/2, and the correspond- 
ing divergence of the particle energy is = 7|f For ^ LiH crystal, the 
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diffraction plane (220) the values of the threshold normalized brightness of 
the beam, which correspond to the generation threshold for magnetic, optical 
undulators and channeled particles, are given in Table [TJ According to Table 
[H the value of brightness in the case of two-wave distributed feedback is rather 
high. But, as it was shown for a parametric quasi- Cherenkov generator, using 
the surface multi-wave diffraction for the formation of disturbed feedback, one 
can decrease the threshold characteristics of a particle beam and provide the 
achievement of the generation regime. 

In [57] the underthreshold spectral-angular distribution of radiation was an- 
alyzed and it was shown that the observation of collective radiation by rela- 
tivistic oscillators was a very complicated problem in the X-ray spectral range. 

Thus, the parallel consideration of two kinds of crystal three-dimensional X- 
ray generators, which are distinguished by the mechanisms of spontaneous 
radiation, shows that three-dimensional distributed feedback allows one to 
decrease the current density of the particle beam by several orders of magni- 
tude in comparison with the results obtained in pT|62f63f64f65] . This enables 
one to consider the construction of the FEL in the hard X-ray range as a 
scientific problem of nowadays, which can be analyzed both theoretically and 
experimentally. 

10 Crystal X-ray FEL based on a natural crystal or an electromag- 
netic undulator 

Crystal X-ray FELs based on a natural crystal and on an electromagnetic 
undulator were first considered in [60j. The English version of this paper is 
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Table 



1 



Parameters 


Magnetostatic 


Optical 


Channeled 




Wiggler 


Wiggler 


Particle 


Accelerator 








Energy 


= 5 GeV 


= 290 MeV 


= 500 MeV 


Normalized brightness 


= 3.5 • 10^ 


= 1.7- 10^° 


= 5-109 


Energy spread 


= 2.4 • 10-^ 


1.2 • 10-^ 




Density of current 


= 5.3 • 10^ 


= 1.3 • 10^ 


= 3.3 • 10^ 


Wiggler 








Wavelength 


= 1 mm 


= 5/xm 




Magnetic field strength 


= 17.5 kG 






Laser energy 




= 0.75 gW 




Crystal 








Wavelength of radiation 


= 0.05 A 


= 0.15 A 


= 1 A 


Asymmetry parameter 


= 9 


= 1 




Diffraction plane 


(220) 


(220) 


(100) 



given below. 

Let a relativistic electron (or positron) beam of velocity Uq and the beam ve- 
locity distribution Au move in a spatially periodic medium (e.g., in a crystal). 
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Let a linearly polarized laser pump wave be incident onto the beam along the 
direction Hp. Let the wave have the wave vector kp = kpn, frequency Up, and 



phase. The z-axis is chosen directed along the beam's average velocity uq. The 
presence of the electromagnetic pump wave induces radiation leading to vari- 
ous sorts of instability. The most known of them is the so-called the three- wave 
parametric instability, which emerges due to the conversion of the pump wave 
into the Doppler-shifted electromagnetic wave and the charge density wave: 
Up = u + Uch, kp = k + kch k are the frequency and the wave vector of 
the high-frequency wave, respectively, Uch, kch are the frequency and the wave 
vector of a charge density wave). This process is scrutinized for the case of 
homogeneous matter |12U] . 

In the case of a 3 D-periodic medium that we are considering here, the situation 
is basically different from that studied earlier because for co and k satisfying 
the Bragg condition, the wave scattered by the beam is diffracted. 

Maxwell's equations and the equations of particle motion, which describe the 
process under study, have the form: 



the field strength Ep = 



coslkpf — ujpt + (5), where 5 is the initial vibration 



^ IdD A'K^ 

lotH = -— H j 

c at c 




oo 



— oo 



lOtE 



IdH 
c~dt' 



divE 



47rp, 





(70) 



p(f, t) 



dt 



rri'-yi mcji 



—E{fi{t),t) + — 



e 



{v,{t){vi{t)E{r,{t),t))), 
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— * — * 

where E{ri{t)),t), H{ri(t),t) are the electric and magnetic field strengths at 
moment t at the location point of the i-th electron fi{t), D is the electric 
induction, e{f, t) is the dielectric permittivity of the crystal at point r, e{f, t) = 
at t < 0, Vi{t) is the velocity of the i-th electron of the beam at time t, j{r, t) 
is the beam's current density at point r at time t, p(r, t) is the density of the 
beam's electron charge, e is the electron charge of the particle (e = ±|e|), 
7i = (1 — f^/c^)^^/^ is the particle Lorentz factor, c is the velocity of light. 

For the Fourier transforms -D(r, lo) and E{r, cj), we have the relation -D(r, uj) = 
e{f, u)E{f, u), where the crystal's spatially periodic dielectric permittivity can 
be represented as the Fourier series: e(r, cu) = J2^t{^)^^'^'' , ^ is the reciprocal 
lattice vector. Studying the possibility of occurrence of induced radiation, one 
should first of all determine the gain coefficient or, which is the same, the beam 
instability increment. To find them, in solving (TfOl) we may confine ourselves to 
linear approximation, finding the dielectric permittivity of the beam+crystal 
system in the presence of the given pump wave. Take the Fourier transform in 
time of ( iTOl) . To be more specific, let us further consider the case when a single 
diffracted wave T appears through diffraction. Let us also assume that 

the scattering plane of the wave, i.e., the plane defined by vectors kp and k 
coincides with the diffraction plane defined by fc. A;,-. In a similar way as in the 
case of dynamical diffraction of X-rays, the generation of tt- and a-polarized 
waves in the crystal in the absence of the beam can be considered separately. 
Recall that the polarization vector of cx-polarized waves is orthogonal to the 
diffraction plane 




We also assume that the polarization vector Cp of the pump wave is parallel 
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to Co-. In this case, upon taking the Fourier transform, (1701) takes the form: 

{k'^c^ — eQu'^)E„{k, u) — erU}'^E„{KT-, u) = A-niujj„{k, w), 
' {ky - eoUJ^)E„{kr, oo) - e-rOJ^E^{k, u) = 0, C^^) 
ujE\\{k,u) = — 47rij|| (/c, cj), 

where || denotes the component of the field and current strength parallel to 
the wave vector k, E^j = (Ee^). In writing (I7T!) . the components j^{kr,co) and 
j\^{kr,uj) of the current, for which the synchronism condition does not hold 
in the X-ray range and whose contribution is therefore small, are discarded. 
Without channeling, the particle velocity in the beam can be represented as 
a sum: 

Vi(t) =Ui + Vpiit) + Svi„(t) + Svi\\(t), 

where Ui is the particle's longitudinal velocity in the absence of the electro- 
magnetic wave. 

e^p - _ ^ - _ 

Vpiit) = — sm{kpfoi - Upt + 6) = vz sm(A;„foi - Upt + 6) 

is the velocity of the i-th particle under the external electromagnetic pump 
wave, ZJp = ujp — {kpUi), 5vi^{t) is the velocity disturbance of the z-th particle in 
the field of a transverse scattered electromagnetic wave, 5vi\\ (t) is the velocity 
disturbance induced by a longitudinal wave. The position of the i-th particle 
at time t is written accordingly rj(t) = + Ui{t) + Sfi{t), where r\)i is the 
position of the i-th particle (e^) in the beam at time t = 0, Sfi{t) is the time 
change of the position of the i-th (e^) in the beam. 
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The Fourier components of the transverse and the longitudinal current in the 
linear approximation with respect to perturbation can be presented as follows: 



Uk, to) = i^^^{kE{k, CO)) + C-^^iik - kp)E{k -kp,u;- co^)) 



+e{Y, 6vi^{u - ku^)e-'''''°^ 



(72) 



e{v^nk)no, 



+ { 



{{k — kp)E{k — kp,u — Wp)), k 7^ kp. 



(73) 



( ) denotes averaging over the velocity distribution in the beam. 

To obtain a closed system of equation, let us take the Fourier transform of the 
equation of particle motion (170]) and separate the transverse components of 
the velocity and the perturbation field from the longitudinal ones. The linear 
approximation to the disturbance of the particle velocity gives 



^'^^oA-/':^SE^ik,u)) + (^ 



X 



+ (i?e.) 

Up — CO 



2(^(fc-M) 



+2{v^ea){{nkVi))6E\i{k - kp,uj - Wp)| 



+ {ui + ea){uinK)SE\i{k,u) 



5Efj{k — kp,co — cOp) (74) 
[uie„f5E^{k,uj) 



mc^to'yi 
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leriQ 



teuQ 
mcZJ'ji 



(75) 



X 



Up — u 



{^t^e^)6E^{k - kp,uj - Up), 



where u = u — kui, (...) is averaging over u. 



Substitution of the obtained expressions for the velocity disturbance into f[72|) . 
(1751) and then into flTT]) gives a closed system of equations for transverse and 
longitudinal fields under two- wave dynamical diffraction of scattered radiation: 



{k'^c^ - eouj'^ - {ujlj_))SEaik,uj) - erCO^SE^{kr,uj) 
= Gi6E„{k,uj) + G26E^{k,u), 

{ky - Eou^ - {ul^))5E,{K,uj) - e^r^HE,{k,u) = 0, 



(76) 



where 



47re^no 



(7) me{-i) 



hq is the unperturbed particle density in the beam, me is the mass of 



Gi = { . 2/-2 -/cc 



(77) 



\ 4c2[(a7-cJp)2-u;2^] 



{kpnkp-k)c - Up{f3nkp-k) 



The determinant of fITBl) defines the dispersion equation that specifies the 
relation between k and u and enables finding the instability increment of the 
beam. We shall describe the beam by Maxwell's velocity distribution f{u) 
with the effective temperatures Ti, T2, T3 over the x, y , z axes, respectively 
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or, which is the same, with thermal-induced velocity dispersions characterized 
by average thermal velocities VaT = \f^- '■ 

Y -(ux-UQ^)'^ _{^v-noy)'^ _{iiz-tiOz)^ 

f{u) = -— e e e "ir . (78) 

Upon averaging ( !77|1 with the function ( !78|) . we obtain 



(^?e,)^w2^((^-M^fc) 

Oro — 



+ {kc - u0nk))] 1 1 + ) ^ , (80) 



where W{^) = ^1 + ^ / e^^dxj e is the Kramp function (it is tabulated 
in |121] ). uj = CO — kuo, oJp = Up — kpUQ, = klvfj. + fc^t>|y + k^v^rp, = 

{kx kpx)'^vfrp + (^ky kpy)^'^2T ~^ (^2 kpz)'^V^rp. 

Equations fl76|) together with functions fl79|) . flHOl) determine the parametric 
Cherenkov instability and parametric decay instability of the thermal beam. 
From (!79l) . (!80|) follows that in the case when the parameters a; and k are such 
that C ^ 1, the velocity distribution in the beam appears insignificant, and so 
one can analyze an appreciably simpler case of a cold beam. Since in radiation 
of a relativistic particle quanta are emitted in a narrow angular range in the 
direction of the particle motion, then k^, ky kz. As a. result, at the given 
values of Q and Qi, the beam's velocity distribution in transverse direction can 
exceed that in the longitudinal direction. Of particular interest is the radiation 
pattern for parametric decay instability, in which case there is a possibility in 
principle to make the differences k — x — kpx and ky — kyp, appearing in the 
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expression for Qi, go to zero. The transverse velocity distribution in this case 
appears to be insufficient at all. 

In a real situation, due to the finite dimensions of the system and non- 
monochromaticity of radiation, these differences cannot vanish, though they 
can be reduced appreciably. Let us estimate the maximum possible magnitude 
of the increment. Suppose that the parameter C ^ 1. As stated above, in this 
case we go over to the approximation of a cold beam and can use the expres- 
sion f{u) = 6{u — Uq) for the distribution function. The parametric Cherenkov 
instability for this case was considered in |51j (see Section above). The term 
proportional to G2 leads to parametric decay instability of the beam in a 
periodic medium. In the general case, the Cherenkov and decay instabilities 
develop at different frequencies and can therefore be considered separately. So 
here we shall focus on the decay instability, omitting the term containing Gi 
in (1761) for corresponding k and u. In the case of a cold beam, using (!76|) . (!77|) . 
we have the following equation for decay instability: 

[{k'^c^ - uj^eo){klc'^ - uj^eo) - a;^e^£_^](cJ - cJp)^ = ^2(^:^0^ - (81) 

where A2 is the factor of the term [{uj — UpY — w^/,]^^ in G2', (jJpi, is dropped 
as it is small for the particle densities uq ^ 10^^ ^ 10^^ cm^^ in the beam that 
are of interest to us. 

The solution of dispersion equation flHTj) can be found in the weakly-coupled 
waves approximation, which is applicable in this case because a non-linear 
right-hand side of (IHT]) is small. In this approximation, the solution of disper- 
sion equation (IHTj) is sought near the intersection of solutions of the dispersion 
equations for coupled waves into which flHTl) is split when its right-hand side 
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equals zero, i.e., 



(82) 

((w — kuo) — UpY = 0. 



Using weakly-coupled waves approximation of the perturbation theory and 
assuming that Eq, e±r are real values (i.e., neglecting the intrinsic absorption 
of the medium), the upper estimate of the increment of the beam's parametric 
decay instability when a pump wave is scattered in a 3D periodic medium can 
be obtained from (IHTl) : 



4 



c \ Auj"^ ujb 



Erl, (83) 



where -^^^ = i i, \^ , cjb is the Bragg frequency. Equation ( l83ll is obtained 
in the assumption that the synchronism condition is fulfilled exactly and that 
D{kz,k±,u) = in the degeneration point of the eigensolutions of the disper- 
sion equation of diffraction. In this case the multiplicity of the root, determin- 
ing the increment, is greater by one than in the case of modified decay in the 
absence of diffraction of a scattered wave. This results in the decrease in the 
instability threshold, depending on the beam density and the intensity of the 
external electromagnetic pump wave |118fll9] . According to (!83|) , the numer- 
ical estimate of the maximum increment of the parametric decay instability of 
a cold relativistic beam is as follows: rj ~ 0.5 at ~ 10^ in CGS, 7 = 2- 10^, 
co^ = 1.6-10^^ s-\ Up = 10^^ s-^no-lO^^ cm-3, £0-1 = 5■10-^ ~ 10-6. 

In the above theoretical analysis we did not actually use the explicit form of 
Vp and the direct relation between ujp and kp. This theory is therefore fully ap- 
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plicable to the case when, instead of the electromagnetic wave, crystal planes 
bent under the ultrasonic wave act as a dynamic wiggler. A detailed treat- 
ment of the particle motion in such an ultrasonic wiggler was given in [T3] . 
The corresponding dispersion equation is obtained from fl76|) . f l77|) . flHTj) by 



ultrasonic wiggler, fi^ is its frequency, r%j_ is the amplitude of the particle's 
transverse vibrations in the ultrasonic wiggler [13]) and by substituting k for 
kp] Qg, for Up] Q'g, for uJp. Upon this substitution, for the instability increment 
in an ultrasonic wiggler, we obtain the expression similar to fl83|) . This en- 
ables readily obtaining the following ratio of the instability increment in an 
ultrasonic wiggler to that in a light wiggler: 



The analysis shows that ( l84l) can be greater than unity. In a real situation, 
absorption can always reduce by one or two orders of magnitude the upper lim- 
iting estimate given above, but the conclusion about a significant magnitude 
of the gain coefficient holds true. 

11 Volume Free Electron Laser 

As has been stated, the features of radiation from relativistic particles in 
crystals, which were considered earlier in this review, have a general character 
and are also manifested when radiation is generated in artificial crystals (at 
present termed "photonic crystals") [25]l60f87ll88|l89|lll7yi64|ll81|ll76] . 




substituting for Vp the corresponding expression for the case of the ultrasonic 
wiggler \Vp\ = (eo-r^j_)fi'^ (where Q'^ = k^Uq — Qs, is the wave vector of the 




(84) 
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In contrast to generation in the X-ray range, generation of radiation in mi- 
crowave and optical ranges does not require such high particle current densi- 
ties. As a result, application of photonic crystals made it possible to develop 
a new type of generator, called the Volume Free Electron Laser (VFEL). Its 
main features will be considered below. 

12 Generation equations and threshold conditions in the case of 
two-wave diffraction 

Using multi-wave diffraction in a VFEL for the formation of a volume dis- 
tributed feedback enables one, on the one hand, to appreciably reduce the 
length of the interaction area, and, on the other hand, to employ electron 
beams with a large transverse cross-section for pulse generation, which im- 
proves the electrical endurance of the generator and prevents burning of the 
glass in the discharge tubes during high-power lasing. Besides, multi-wave 
diffraction provides the selection of radiation modes in oversized systems j^lOSj . 

When radiation is generated in an FEL, the electrons interact with the elec- 
tromagnetic wave in a finite spatial region, and release energy into the wave. 
Depending on the length of the spatial region, different generation regimes are 
realized. At llmfc^Ll ^ 1, generation occurs in a strong (exponential) ampli- 
fication regime, which is mainly employed in amplifiers. The regime of weak 
single-passage amplification [\ImkzL\ < 1) is chiefly used in generators. In 
order to determine the structure of the fields and to describe the evolution 
of instability in such systems, in addition to the knowledge of the dispersion 
equations and their solutions, one should match the fields at the boundaries 
of the regions (joining solutions). This procedure gives the field distribution 
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in the system. 

Now we shall formulate the boundary problem. Let an electron beam with 
mean velocity u be incident onto a plane-parallel spatially periodic plate of 
thickness L. The electron beam is oriented so that the radiation generated 
by the beam will be under diffraction conditions. Under two-wave diffraction, 
two fundamentally different geometries are possible. In the first case (Laue 
geometry (see Figure [8]) both waves are emitted through one and the same 
boundary of the periodic structure (7071 > 0, where 70 = -^^-^ and 71 = 
■^^^ are the cosines made by the wave vectors k and kr with the normal 
vector to the surface of the periodic medium). In this geometry the amplifying 
regime is possible. 

T 

Figure 8. The geometry of two-wave Laue diffraction k, kr are the wave vectors of 
the incident and diffracted waves, r— is the reciprocal lattice vector of the periodic 
structure. The projections of both wave vectors onto the direction of the normal to 
the surface have the same sign. 

In the latter case (Bragg geometry (see Figure [9])), the incident and diffracted 
waves leave the plate through the opposite surfaces (7071 < 0), when the 
electron current emits photons, positive feedback appears, and the generation 
regime is available. 

Let us write the expression for the fields appearing in the system described 
here 
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Figure 9. The geometry of two-wave Bragg diffraction, k, k^- are tlie wave vectors of 
the incident and diffracted waves, respectively; r is the reciprocal lattice vector of 
the periodic structure. The projection of the wave vectors onto the direction of the 
normal to the surface are opposite in sign. 

ae exp (ikorj + hCr exp (ikoT^j (/) 

Cj exp (iki-rj (e + SjcV exp (irr)) (//) 

i 

/e exp (ikQ-rj exp {—ik^zL) + eV exp (ikoT-rj exp {—ikorzL) (HI) i 

(85) 



aeexp (ikof^ + gcr exp (^zA^q,- (/) 

Cj exp [ikiT^ (e + SjcV exp (irr)) (//) (86) 

i 

/e exp (ikofj exp {—ikozL) + 6eV exp (^zA;ot- exp [ik^rzL) {m) 

In (jHSD, fISB]) (/) and (///) denote the fields before and after the periodic 
structure. In Laue geometry (1851) . in the general case two waves with specified 
amplitudes are incident on the system: one along the direction with the wave 
vector k and amplitude h and the other one along the direction with the wave 
vector k-r and amplitude h. They propagate from one side (in this case to the 
surface ^ = 0). In Bragg diffraction geometry fl86|) . the waves are incident on 
the system from opposite sides (the incident wave of amplitude a falls onto 
the surface z = 0, and the diffracted wave with amplitude 6, onto the surface 
z = L). Wave vectors k^ and fcor satisfy the standard dispersion equation 
describing propagation of electromagnetic waves in vacuum (/c^c^ — cu^ = 0). 
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(//) denotes the field in a periodic medium. The dispersion relation in the 
medium is defined by f llOOp . 



is the coupling coefficient between the diffracted and the incident waves for 
the i-th mode determined from equations ( j85l) . ( !86l) and dispersion equation 

(unnD. 

In the case of a "cold" electron beam, dispersion equation (llOOp takes the 
form: 



= -d. ' {k\^ _ - c^.f)) (87) 

As a result, sum (//) consists of six terms because equation (!86|) is of the sixth 
order with respect to k^. But if the radiation geometry is not plane and the 
dielectric susceptibility is small, then the two waves mirror-reffected from the 
target surface may be ignored. In the case of a "hot" electron beam dispersion 
equation fllOOp reads: 



\^ujf I ue 





2 



(2)A , A 



00 r 



- Xiexp(-x2) {ec'-co'){ky~co¥^^) (88) 



4S \ 

Sum (//) now contains only four terms because the corresponding dispersion 
equation is of the fourth order. With the mirror-reffected waves being ne- 
glected, only two terms remain, and the system of the boundary conditions 
for deffning the unknown coefficients takes a very simple form: the field in the 
periodic structure is completely determined by the two boundary conditions 
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for the fields of the incident and diffracted waves at the boundaries of the 
periodic structure. 



Cl + C2 



SlCi + S2C2 



(89) 



Cl + C2 



b = SiCi exp (ikizL) + S2C2 exp {ik2zL) 



(90) 



Equation ( 189|) gives the continuity conditions for the incident and diffracted 
waves at the boundary 2; = for the Laue diffraction case. For the case of Bragg 
diffraction, equation ( |90|) gives the boundary conditions for the incident and 
diffracted waves, respectively when 2 = and z = L. In the regime of a "cold" 
electron beam, provided the mirror-reflected waves are neglected, dispersion 
equation fl86l) has four solutions. To define the structure of the field appearing 
in a periodic plate, four boundary conditions should be used. Two of them 
have the same form as those in the case of a "hot" beam: the continuity of 
the incident and diffracted waves at the target boundaries. Two additional 
conditions are the continuity of the charge and current densities at the input 
boundary. As a result, the system of equations defining the fields in this case 
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has the form 
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E Ci 



a 



1=1 
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(91) 
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1)E 



or 



2) E -SjCiexpli/cj^i^} = 6. 



i=l 



The first equation in ( 19T|) expresses the continuity of the incident wave at the 
front boundary, the third and second ones describe the continuity of the charge 
and current densities at the front boundary, 5i = The fourth equation 

expresses the continuity of the diffracted wave value at the front boundary in 
Laue geometry (denoted by figure 1 in hne four of fl9T|) or at the back boundary 
in the case of Bragg geometry (figure 2 in hne 4 of ( 19T1) ). It is easy to see (see 
( 189|) . ( 19T|) ) that in Laue diffraction geometry, the electromagnetic field in the 
interaction area will also be absent if the amplitudes of the incident waves 
equal zero. 

The situation is different in Bragg geometry. From (!90|) . ( 19T|) follows that in 
this case, under certain conditions it is possible that the field in the medium 
exists at nonzero amplitudes of the incident fields. The generation equations 
defining these conditions are obtained by equating the determinants of the 
systems fl90l) and fl9T]) to zero. For a "hot" electron beam, the generation 
equation takes the form: 

S2 exp {ik2zL) — si exp {ikizL) = 0. (92) 
A similar equation for a "cold" electron beam reads 
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Si exp{ikuL}Sl (S2 - S3) {S2 - S4) {S3 - S4) 
-S2 eyip{ik2zL}Sl {Si - S3) {Si - S^) {S3 - S4) 

+S3 exp{ik3MSl {Si - S2) {Si - S,) {S2 - S,) (93) 
-S4 exp{ik4^L}Sl {Si - ^2) {Si - S3) {S2 - S3) = 0. 



Upon solving these equations, one may determine the threshold generation 
conditions, i.e., the electron current and other parameters of the beam, mark- 
ing the starting point from which radiation prevails over the losses. Besides, 
radiation generation takes place when some phase relations are fulfilled: phase 
shift between the two diffraction modes traversing the interaction area should 
be a multiple of 2n (the field at the output should be similar to that at the 
input): 



Re {kiz — k2z) L = 27rn. 



(94) 



In the region, where these conditions are fulfilled, the solution of fl93l) . fl94l) 
has the form: 



2(1-/3) 



\Xt\Xo 



7oC 



IQtt n 



2^2 



nuJ -/3{kxTL^) kU 



bj" in (1951) is the increment of absolute instability, which describes the increase 
of the field amplitude in time in the field linear regime. 



.VZ^(a!^!4-ixWexp(-xW2 



for a "hot" beam 



^ {^S^- (^(Pc2 - co^)k^Llf{y) for a "cold" beam 



. (2?/ + vrn) siny - y {y + irn) cos y 
y-^ {y + nn) 



(96) 
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is the profile function depending on tlie detuning from synchronism conditions 




Figure 10. The profile generation function in the root degeneration region versus 
the detuning from synchronism conditions y 

The profile function of detuning is plotted in Figure HOl As follows from ( l96ll 
and Figure [TOl one of the peculiarities of generation in the region of degen- 
eration of the roots of the diffraction equation is that stimulated radiation 
is nonzero at zero detuning. This peculiarity appears due to the interference 
of contributions to radiation coming from the two branches of the diffraction 
equation. In a standard classical FEL, generation does not occur at the exact 
fulfillment of the synchronism condition. This is due to the fact that the gain 
coefficient is proportional to the difference between spontaneous radiation and 
absorption. In the case when these branches overlap each other, the gain co- 
efficient becomes proportional to the derivative of the spectral function, while 
at zero detuning from synchronism conditions the spectral intensity of sponta- 
neous radiation reaches its maximum, and the derivative at this point equals 
zero. 
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Equation ( 195|) has an obvious physical meaning: the first term between the 
braces is proportional to the intensity of radiation produced by the electron 
beam (in the unit length). The next two terms describe losses of radiation 
appearing through absorption and due to the fact that radiation leaves the area 
of interaction with electrons. The values of the electron beam parameters, at 
which the radiation generation equals radiation losses [u" ) mark the starting 
point for the generation process. From the obtained expression ( p5|) follows 
that at given values of the size of the interaction area and absorption, there 
is an optimal diffraction geometry in which the losses Fioss are minimum. In 
the general case, the optimal geometry is not at all plane. Figure [11] plots the 
relation T\oss iP) /rioss {P = ~1) against the asymmetry factor for diffraction 
in the millimeter wavelength range (A ~ 4 mm). 
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13 Generation equations and threshold conditions in the geometry 
of three-wave Bragg diffraction 

Above we have considered the theory of a volume distributed feedback genera- 
tor in the geometry of two-wave Bragg diffraction. It has been shown that the 
transition to a volume distributed feedback opens up wider perspectives for 
the optimization of the system. Transition to a multi-wave distributed feed- 
back provides additional possibihties. Here we shall mention some of them. For 
example, in a hard X-ray range, where the susceptibility has a negative value 
Xo < 0, the condition of fulfillment of Cherenkov synchronism imposes restric- 
tions on the asymmetry factor of diffraction \/3\ > 
X-ray spectral range, where the absorption is large, radiation losses under the 
condition of two- wave dynamical diffraction will be rather large at such values 
of the asymmetry factor. Moreover, strict requirements are set for the param- 
eters of the starting generation current. One of the ways to diminish the losses 
is using supplementary external mirrors. Transition to multi-wave diffraction 
also enables one to reduce losses because the parameters of the system re- 
quired to initiate the generation process change due to the rearrangement of 
the field structure in the interaction area. Synchronism conditions in this case 
contain additional parameters as compared to the two- wave dynamical diffrac- 
tion, which enables one to match the Cherenkov conditions to the parameters 
corresponding to the regions with smaller absorption of radiation. 

The possibihty of generation at the point of degeneration of several diffraction 
roots is another important feature of multi-wave diffraction. At this point the 
functional change of the threshold characteristics occurs (e.g., the coincidence 

of the two roots gives the doubly degenerate case and s = 2) and, as a result. 



Xo 



+ 7-')/IXr|'. In the 
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it becomes possible to reduce the length of the generation area (at a given 
operating current) or the operating current (at a given length of the generation 
area) . 

Application of multi-wave diffraction for generating in a microwave range has 
one more remarkable feature - the possibility of selecting modes in oversized 
waveguides and resonators. Production of microwave power pulses requires 
high electric strength of the generator and radiation resistance of the output 
window. To reduce load on these elements, the transversal (with respect to the 
direction of the electron beam velocity) dimension of the resonator should be 
large (much larger than the wavelength). As a rule this leads to a multi-mode 
generation regime and low efficiency. In the presence of n-wave diffraction, 
selection of modes can be effectively carried out due to the requirement of 
fulfillment of n — 1 Bragg condition. 




Figure 12. Figure Three- wave diffraction in Laue-Laue geometry. Projection of wave 
vectors k, kri, kr2 onto the normal to the surface have the same sign. ti,T2 are the 
reciprocal lattice vectors of the periodic structure. 

To illustrate the potential of the multi-wave distributed feedback, we shall give 
a more detailed consideration of three-wave diffraction. In this case, volume 
distributed feedback (VDF) can be realized in three different geometries: 
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(1) Laue-Laue diffraction, when the three waves exit through the same surface 
(70,71,72 > 0, Figure [121 

(2) diffraction in Bragg-Bragg geometry, when 70 > 0, while 71,72 < 0. That 
is, the wave vector satisfying the synchronism condition is directed along the 
normal to the surface of the slow-wave structure, and the two vectors corre- 
sponding to diffracted waves point in the opposite direction; 

(3) Bragg-Laue diffraction, when 70, 71 > 0, 72 < 0. Vector k corresponding to 
the synchronous wave and one of the vectors ki = k + Ti corresponding to the 
diffracted wave point in the same direction, while vector k2 = k + T2 points in 
different direction (see Figure [T^ . 



T2 




Figure 13. Figure Three-wave diffraction in Laue-Bragg geometry. The projection 
of wave vectors k, k^i onto the surface normal and the projection of vector kr2 are 
opposite in sign, -fi, T2 are the reciprocal lattice vectors of the periodic structure. 

Similarly to the two-wave case, the problem of the beam interaction with the 
target field may be reduced to the problem of three-wave diffraction of an 
electromagnetic wave incident onto the active medium. The active medium 
here is the system "spatially-periodic structure -|- electron beam". 

The system of equations describing the process of three- wave coplanar diffrac- 
tion by such an active medium takes the form: 
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Z^(o)^(o) 








(97) 



^ Xl-2 



Ei'^ + = 0, 



where the terms 



; 2 2 2 (i) I 

k-c - ijj e\ + 



7 



+ Tj ( A;, 



- (w - a;„)' 



(9J 



that contain the contributions due to resonance interaction of electrons with 
the wave, 



r, = <^ 



7 (uj—kiU—uJo')'^ 



in the case of a " cold" electron beam 



7 <5W 



in the case of a "hot" electron beam 



(99) 

The values of Fj in equation (|99|) are given for two generation regimes: in the 
regime of a " cold" electron beam all the electrons have the same velocity and 
the velocity spread is small (^7 ^ 1, where Am is the thermal straggling 
of the electron beam, Aw is the width of the emission line). In this case all 
the electrons participate in the interaction with the electromagnetic wave. In 
the regime of a "hot" electron beam the relation > 1 is fulfilled. In this 
case only a part of the electron beam participates in the interaction process. 
The explicit form of f2 and depends on the generation mechanism. For the 
Cherenkov mechanism = (^^^ , where Cj is the polarization vector, = 0. 
In (197 



2 47rne 



UJ — km — u, 
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are the thermal stragghng for a "hot" electron beam in the distribution func- 
tion: 



is the distribution function describing the velocity spread of the electrons in 
a beam before the interaction with the electromagnetic wave. 

In the case of two-wave Bragg diffraction, when only two wave vectors sat- 
isfy the Bragg condition |3] and two strong waves are excited, the dispersion 
equation for Cherenkov instability takes the form 



Here r = XtX-t 

As it was noted in |19|50] . the functional relationship between the imaginary 
part of the solution and the electron beam density may change apprecia- 
bly at the points of degeneration of roots kz of the diffraction equation (the 
dispersion equation without the electron beam describing the dispersion of 
electromagnetic waves in a periodic medium). This occurs, in particular, at 
the point of s-fold degeneration of roots: Imk^ ~ ny^^^+s) the Compton 
generation regime, this quantity significantly exceeds in magnitude a similar 
parameter for Irak^z ~ nl^^ beyond the root degeneration region. From this 
fact the authors of [I9f50] concluded that the instability increment in the 
degeneration region grows sharply. 

The dispersion equation corresponding to the system ( 1971) for the three-wave 





(100) 



case of coplanar diffraction reads: 



(101) 



where 



fci, k,) = (k\^ - ujh,) (klc^ - uj\,) {klc' - u'eo 

-uj'^ (klc^ - xiX-i - uj^ (X1X-2X2-1 + X2X-1X1-2) , 
(^1, ^2) = {klc' - uho) {klc^ - ujh,) - a;^Xi-2X2-i. 

In deriving fllOip . we based on the assumption that the electrons are syn- 
chronous to the wave with wave vector k. The field in the three-wave system 
is written as follows: 

(1) in Laue-Laue geometry 

a'-^-'e exp (^ikorj + a^^^ei exp (ikoirj + a^'^''e2 exp (iko2rj 
Cj exp (iki-rj (e + SjcV exp {irr)) 

i 

beexp (ik^rj exp {—ikQ^L) + b^^^ei exp (ikoivj exp {—ikoizL) 
+}P'^e2 exp (iko2rj exp {-iko2zL) , 



(2) in Bragg-Bragg geometry 

a'-^-'eexp ei exp 62 exp (^iA;^2 

Cj exp (jkivj (e + SjcV exp (zrr)) 

j 

6e exp (iko-f^ exp {—ikozL) + a'-^'^^ei exp (^A^oi ^'^) ^xp (ikoizL) 
+a^^^ei exp (iA;Qi exp (ikoizL) , 
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(3) in Bragg-Laue geometry 

a'-^-'e exp (jkovj + a^^^ei exp (ikoirj + }P'^e2 exp (^^02''^) 
Cj exp {ikiT^ (e + SjeV exp (zrf)) 

i 

beexp (ikorj exp {—ikozL) + b^^^'ei exp (ikoirj exp (ikoizL) 
+a^^^ei exp (iA^oi ''^) ^^P (^^oiz-^) 5 



In the above expressions the waves incident on the structure have the am- 
phtudes a^^\ while the exit waves have the amphtudes b^^\ In the expression 
for the field in a periodic structure, summation is made over the diffraction 
modes. In the case of a "cold" beam, the sum contains five terms, and in the 
case of a "hot" beam, three terms. 



(1) AiAi2 — r2 (2) AjAii — ri 

Si = and s — 



2 2 



Ai2Xi + X2X1-2 " AiiX2 + X1X2-1 



are the coupling coefficients between the diffracted waves and the incident 
wave; their form is obtained by the system of equations ( 1971) . Xia = {{ki + 

By analogy with a two-wave case, we shall join the fields at the system bound- 
aries and write the system of equations defining the field in the interaction 

(1) in the case of Laue-Laue geometry: 



55 5 

i=l i=l i=l 

Ef = o E| = o, (102) 

i=l i=l "i 
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(2) for Bragg-Bragg geometry 



5 



= a 



i=l 



i=l 



i=l 



y^ = o V — = 

i=l i=l 



(103) 



(3) for Bragg-Laue geometry 



55 5 



1=1 i=l i=l 



i=i i=i 



(104) 



In all three of these systems fll02l) - fll02p the first three equalities are obtained 
from the requirement of the continuity of fields at the boundary. In the Laue- 
Laue case the boundary conditions are written for the boundary z = 0. In 
Bragg-Bragg geometry one condition for the incident wave is written for the 
boundary z = 0, two others, for diffracted waves for the boundary z = L. In 
Laue-Bragg geometry two conditions are defined at z = 0, and one, at z = L. 
The last two equations of the obtained systems follow from the requirement of 
continuity of the current and charge densities at the boundary 2; = 0. In the 
three-wave case of Laue geometry, as well as in the case of a two-wave VDF, 
only the amplification regime is possible. In Bragg-Bragg and Bragg-Laue ge- 
ometry both the amplification regime and the generation regime accompanied 
by radiation from spontaneous noise at zero amplitude of the incident field are 
possible. Using the systems (11031) and (11041) . write the generation equations 
for these two geometries. 

Bragg-Bragg geometry 
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(4'^4'^ - 4'^4'0 -^i^s (53 - S,) exp{z (fci, + k2.) L} 

+ [s^^h? - 4'^4'^) '5i'54 (^2 - ^3) exp{t {h, + h,) L} (105) 

+ {s2^s'^^ - -53^^ -s^^^) 52<^3 (<^i - <^4) exp{i {k2z + k3z) L} 

- ('5?^4^^ - s^^s'^?) ^2^4 (^1 - ^3) exp{i (A;22 + /i;42) ^} 

+ (4^)4^) - 4^)4^)) 53^4 (5i - ^2) exp{z (^3. + k,,) L} = 0. 

Bragg-Laue geometry 

4'^ exp (zki^L) 6i{sS^62 {S, - 6s) - s^^^6, {6, - 62) + 4'^'^4 (^3 - ^2)} 
-4') exp {zk2.L) 52{4'^5i (^4 - 53) - 4'^<^3 (^4 - S^) + 4^)54 (^3 - ^i)} 
+4'^ exp [tks^L) 53{4'^5i (^4 - 62) - 4'^52 (54 - 5i) + 4'^54 (52 - 5i)} 
-4') exp (tk^^L) 54{4'^5i (53 - 52) - 4'^52 (53 - 5i) + 4^)53 (52 - 5i)} = 

(106) 

Similar generation equations in these geometries may be written for a "hot" 
electron beam |105] . 

Before we proceed to the analysis of generation equations, let us examine the 
behavior of the roots of the dispersion equation for three-wave diffraction. It is 
known that the most effective interaction of the electron beam and the electro- 
magnetic wave takes place at the root degeneration point when the synchro- 
nism condition is fulfilled. Thus, the parameters for the generation threshold 
can be determined by the weakly-coupled mode method |106j . According to 
this method, we first find the solutions of the dispersion equation without 
the electron beam, then the additional synchronism condition is imposed on 
these solutions. To complete the procedure, one should substitute the value of 
k = ko that satisfies the condition u — k^u = into the dispersion equation 
describing three-wave diffraction. Making use of the form of Fi^) (k, /ci, A;2 ) in 
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(dni]), we get: 



llik - Iru - hr2 - hn - / = 

f^Mk + I {Pih + P2I2) - /3i/32ri2 - - P2r2 = 0. (107) 



In dinn!) 

/ = {klc^ - uj\^) lu? = sin2 e + (^-^ j 



2 

- Xo, 



fl,2 = f + ai,2, "1,2 = TT-^ 



are the two parameters of deviation from Bragg diffraction conditions, / = 
X1X-2X2-1+X2X-1X1-2, /3i,2 = 70/71,2 are the asymmetry factors of diffraction, 
70,71,72 are the direction cosines of diffraction. In fll07p the first equation is 
the requitement for simultaneous fulfillment of the synchronism conditions and 
the dispersion equation of diffraction F^^^ (^k, ki, k2^ = 0. The second equation 
in (11071) appears due to the requirement for root degeneration of the diffraction 
equation. Introducing new variables Zi = li/l and Z2 = h/i, we obtain from 
dMl): 



r2Zi + riZ2 ri2 / , 
Z1Z2 7^ - /3 = (108) 

f3i(32ZiZ2 + PiZi + P2Z2 - /3i/32^ - ^^^^ ^2 ^^^^ = 0- 



Solve (dOHD for z^: 



ri - P2XiX2^ ± Xi-2 (1 + ^) y'-|(/2 + /3,ri + /32r2) 



P + f32r2 



_r2- PiXiX2^ T Xi-2 (1 + g^) yZ-l jl' + + /32r2) 



(109) 
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Thus, at the point of intersection of roots, the parameters of deviation from 
the Bragg conditions ai 2 appear to be expressed in terms of the angle of 
the radiation wave vector with respect to the velocity vector (recall here that 
I = sin^ 6* + (c(m7))^ — xo), the asymmetry factors of diffraction (/3j) and 
the polarizability (xi) of the periodic structure. From fll09p follows that at 
the point of degeneration of diffraction roots and simultaneous fulfillment of 
Cherenkov synchronism, the following condition should hold: 

/9i/32(/' + /3iri + /32r2) < 0. (110) 



From (11 101) follows that there are no degeneration points in Laue geometry 
(ri,r2 > 0, while for the Laue case doth asymmetry factors f3i and (32 are 
positive values). In Bragg-Bragg geometry, a critical angle exists: 



sin6, 



thr 



\ 



\ 2 

C 



-/3iri-/32r2 + Xo- — • (HI) 
\ M7 / 



The degeneration region occurs at radiation angles 6 < 6thr- So, in the X-ray 
spectral range in Bragg-Bragg geometry at least one of the asymmetry factors 
should be large. In Bragg-Laue geometry there is the opposite situation, when 
the degeneration region occurs at 6 > ^thr- In three- wave geometry there 
is a possibility of threefold degeneration of roots. At the point of threefold 
degeneration, there is a strong relationship among the direction of the photon 
emission, asymmetry factors, and the polarizabilities of the periodic structure. 

The conditions for the threshold current values are obtained by through solv- 
ing equations (llOSp and fll06p . In a weak single-stage amplification regime 
{llmkzLl ^ 1), this conditions in the region of two-fold degeneration have the 
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form: 



G'" = W + ^ jrj^\nBUss,\ (112) 

|/3i/32| \nuj {klL^) kL^ 



In (I112p the following notations are used: 

^1^2 + ^1 + ^2 ^2 1 ^^^77 r 

^12^2 ^ 

r r -U r -U r n+r-2+ri2 , ri2+zir^^+Z2r{^ ^ 
Z\Z2 + Zi + Z2 72 1 —77 h TT- 



^ -y ''12 ' 
2^1^:2 

S.'^ SI — S'- 

rjBB = X- 



/ (1) (2) \ 
^1 — ^^2 + 1^S2 ^1 — Si ^2 J 


1 cos{k {d 


{ - 4) 


1^} 


1 


fJ2)„{l) Jl),(2)^ 




1 



^1 = — 72 



r, — ''1 1 X2 



/3l (^ + 



X1X2-1 



^ '•2 \ Xi 
^2 - 72 ) — 



^2 . w 

/9 / ^2X2 _|_ X2XI-2 ^ 

G(^) is defined in ([95D, ([96]). 

The threshold condition in equation f lll2p holds true in the region of the two- 
fold degeneration of roots kiz, ^22, and as the modes corresponding to these 
roots pass through the interaction area, their relative phase shift should satisfy 
the condition {ki^ — k2z) L = 27rn (here the roots are "almost" degenerated 
if {27m) / klxrlL) ^ 1). The third diffraction root is located at a distance 
1^32 — ^1^1 ~ IXrl, 1^32 — ^2x1 ~ \Xt\ from the degenerated roots. 

The dependence of the threshold conditions on the length of the interaction 
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area at the point of three-fold degeneration changes appreciably: 



5 / 27r 



) 



4 




(113) 



Realization of this regime requires that the following phase conditions be ful- 
filled: {kiz — k2z) L = 27m, {k2z — k^z) L = 27rm, n ^ m. The coefficients A, B 
in ( I113P depend on the polarizabilities, ^i, Z2 and the indices m, n. Since they 
are awkward, they are dropped here. From flll3p follows that in the region of 
three-fold degeneration of roots, the functional dependence of the losses at the 
boundaries on L changes significantly. According to (11131) . for a "cold" elec- 
tron beam, when absorption is not important, jthr ~ i^'^l) (^^ would be 
recalled here that in the regime of a "cold" electron beam G ~ jL"^ (see (!96|l ). 
Under dynamical diffraction, when inequality ^ ^ 1 holds, this dependence 
leads to an appreciable reduction of the threshold current. Under multi-wave 
(s-wave) VDF, the threshold current depends on L as 



so the transition to multi-wave diffraction enables one to significantly reduce 
the longitudinal dimension of the generating system. 

As follows from the above results, the volume distributed feedback (VDF) has 
a number of advantages that make its application beneficial for generating 
stimulated radiation in a wide spectral range (with the wavelengths from cen- 
timeters and millimeters to angstroms). Moreover, in a short-wave spectral 
range, the size of the radiating system reduces appreciably due to the change 
of the functional dependence under multi-wave VDF. A short-wave spectrum 
corresponds to greater values of the wave number fc, so at a given operating 




(114) 
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current, due to the presence of the factor yj^j , L can be appreciably 
reduced. In optical and X-ray ranges, where the requirements for the current 
density and the quality of the beam are very strict, there appears the pos- 
sibility to noticeably reduce the threshold values of the current for a given 
beam propagation area. In this case the VFEL is a unique system providing 
lasing at relatively small interaction lengths. In a microwave range, VFELs are 
beneficial for both the reduction of the generator size and selection of modes 
when producing high-power radiation pulses in oversized generators. 



14 Application of volume diffraction gratings for terahertz lasing 
in Volume FELs (VFELs) 

Generation of radiation in millimeter and far-infrared range with nonrela- 
tivistic and low-relativistic electron beams is a complicated task. Gyrotrons 
and cyclotron resonance facilities arc used as sources in millimeter and sub- 
millimeter range, but for their operation a magnetic field of several tens of 
kiloGauss (u ~ ^7) is necessary. Slow-wave devices (TWT, BWT, orotrons) 
in this range require application of dense and thin (< 0.1 mm) electron 
beams because only electrons passing near the slowing structure at distance 
d < A^7/(47r) can effectively interact with electromagnetic waves. It is dif- 
ficult to guide thin beams near a slowing structure with desired accuracy. 
And electrical endurance of resonator limits radiation power and density of 
the acceptable electron beam. Conventional waveguide systems are essentially 
restricted by the requirement for transverse dimensions of a resonator, which 
should not significantly exceed the radiation wavelength. Otherwise, the gen- 
eration efficiency decreases abruptly due to the excitation of plenty of modes. 
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Figure 14. General view of Volume Free Electron Laser formed by metal threads 
with several sheet electron beams. 

Most of the above problems can be overcome in Volume Free Electron Lasers 
(VFEL) 

In volume FELs, the greater part of the electron beam interacts with an elec- 
tromagnetic wave due to volume distributed interaction. Transverse dimen- 
sions of a VFEL resonator could significantly exceed radiation wavelength 
D ^ X.ln addition, the electron beam and the radiation power are distributed 
over the whole volume, which is beneficial for electrical endurance of the sys- 
tem. Multi-wave Bragg dynamical diffraction provides mode discrimination in 
VFELs. 

14-1 Amplification and generation in a photonic crystal 

Let us consider an electron beam with velocity u passing through a periodic 
structure composed of either dielectric or metal threads (see Figure [T4|) Fields, 
appearing while an electron beam passes through a volume spatially periodic 
medium, are described by the set of equations given in |105] Instability of 
electron beam is described by the dispersion equation |105] : 

{k^C^ - Uj\){ky - Uj\ + Xt'^) - ^\rX-r = 0. (115) 
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= k + T is the wave vector of the diffracted photon, r = x"^' T"^} 
are the reciprocal lattice vectors, a, 6, c are the translation periods, xi^) is the 
part of dielectric susceptibility caused by the presence of the electron beam. 
As synchronism conditions are incompatible with those of Bragg, in the in- 
stability range P ^ k^. At the same time, two different types of instability 
exist, depending on radiation frequency. Amplification takes place when the 
electron beam is in synchronism with the electromagnetic component k + f, 
which has a positive projection k^. If the projection k^ is negative and gener- 
ation threshold is reached, then generation evolves. In the first case, radiation 
propagates along the transmitted wave which has positive projection of group 
velocity = {k'f^ = y^w^e — kj_), and beam disturbance moves along 

it. In the second case, the group velocity has negative projection Vz = 
and radiation propagates along the back-wave and the electromagnetic wave 
comes from the range of the greatest beam disturbance to the place, where 
electrons come into the interaction area. For a one-dimensional structure such 
a mechanism is realized in a backward-wave tube. In amplification case, equa- 
tion f lllSp gives for the increment of instability: Imk'^ = —^f, where 




if the condition 2k'^f ^ 



^ xo 



is fulfilled. Here r = XtX-t, h{ue^) /c = I/7 



if the electron beam propagates in a strong guiding magnetic field, otherwise, 

2 " 

h = 1/7. In case 2k'zf ^ ^ J° a dissipative instability evolves. Its increment 
is 



Imky 



\ 
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If inequalities k'^ ^ 2kzk\^ and k'^ ^ " J° are fulfilled, the spatial increment 
of instability can be expressed as 



but the parameters providing such dependence correspond to the conversion 
from the amplification to the generation regime (for the Compton instability 
this situation takes place at kf^ ^ 0). The frequency of amplified radiation is 
defined as: 



TU 

uj = ■ (116) 



The instability in the generation regime is described by the temporal incre- 
ment and cannot be described by the spatial increment. The increment of ab- 
solute instability can be found by solving the equation lmk^^\ijj) = lmk^~\uj) 
with respect to the imaginary part of u. Calculated dependence of temporal 
increment on the parameter of detuning is presented in Figure [151 
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Figure 15. Calculated dependence of temporal increment on tuning out parameter 
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It follows from Figure [T5] that at a certain value of the parameter of detuning: 



/ 



-0.3 (117) 



the increment of instability has a maximum 



- + 



2u}U)"£0+Uj'^X0" 

' ^' ^0.98. (118) 



/ 

The increment of absolute instability can be found from (11181) . The absolute 
instability can evolve if current exceeds start value, which is determined by 
dissipation. The amplification regime has no threshold and decisive influence 
of dissipation causes dissipative instability. Frequencies, corresponding to the 
generation regime, are defined by an expression different from (I116p : 

TU 

uj = . (119) 

Thus, it follows from (I116p . (I119p that the change of a radiation angle causes 
smooth frequency tuning. As a result, the generation frequencies are less than 
those corresponding to the amplification regime. Hence, using the system as 
an amplifier, one should add dispersion elements in it to raise dissipation in 
the frequency range, where generation occurs. 

The use of Bragg multi-wave distributed feedback increases generation effi- 
ciency and provides discrimination of generated modes. If the conditions of 
synchronism and Bragg conditions are not fulfilled simultaneously, photonic 
crystals (diffraction structures) with different periods can be applied |107] . 
One of them provides the synchronism of the electromagnetic wave with the 
electron beam u — ku = fiu. The second photonic crystal (diffraction struc- 
ture) evolves distributed Bragg coupling \k\ ^ \k + Tj\, tj {j = 2 ^ n) are 
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the reciprocal lattice vectors of the second structure. Threshold conditions for 
s-wave diffraction are converted to: 



For dynamical diffraction, when k\x\L^, ^ 1, either the generation start cur- 
rent or the length of the generation zone at certain current value can be 
reduced. 

Each Bragg condition holds one of free parameters. For example, for certain 
geometry and electron beam velocity, two conditions for three- wave diffraction 
entirely determine transverse components of wave vectors and ky, and there- 
fore generation frequency (see flll6p . flllQp ). Hence, volume photonic crystal 
(diffraction system) provides mode discrimination due to multi-wave diffrac- 
tion. 

The above results affirm that a photonic crystal (volume diffraction struc- 
ture) provides both amplification and generation regimes even in the absence 
of dynamical diffraction. In the latter case, generation evolves with backward 
wave, similar to the backward-wave tube. The frequency in such structures is 
changed smoothly either by a smooth variation of the radiation angle (varia- 
tion of k^ and ky) or by the rotation of the diffraction grating or the electron 
beam (change of tu) (see (11161) ). For certain geometry and reciprocal lat- 
tice vector, amplification corresponds to higher frequencies than generation. 
Rotation of either the diffraction grating or the electron beam also changes 
the value of the boundary frequency, which separates generation and ampli- 
fication ranges. The use of multi-wave distributed feedback owing to Bragg 
diffraction allows one either to increase the generation efficiency or to reduce 
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the length of the interaction area (11201) . In this case, generation is available 
with both backward and following waves. In particular, the proposed volume 
structure can be used for generation of sub-millimeter radiation by accelerator 
LIU-3000. The parameters of this setup are: electron beam energy E = 800 
keV, beam current / = 100 -v- 200 A. To generate radiation with wavelength 
of 0.3 mm in such a system, a volume structure composed of strained threads 
should have a period ~ 2 mm, and a period of diffraction grating providing 
Bragg coupling ~ 0.16 mm. Generation of radiation in the terahertz range can 
also be obtained in a photonic crystal using higher harmonics of the Fourier 
expansion of the dielectric permittivity of the crystal [163] . 



15 Dependence of VFEL threshold conditions on undulator pa- 
rameters 

A sharp increase of amplification caused by volume distributed feedback (VDFB) 
yields a noticeable reduction of threshold currents necessary for the lasing 
start. This fact is particularly important for lasing in sub-millimeter and vis- 
ible ranges and for shorter wavelengths. Explicit expressions of the VFEL 
threshold currents were obtained in |105] . Here we shall consider the depen- 
dence of VFEL starting current on the undulator parameters |149j . 

The set of equations describing the interaction of a relativistic electron beam, 
which propagates in the spatially periodic structure of the undulator is [105] : 

DE - u'^XiEi - w^X2^2 - ... = 
-uj\^iE + D,Ei - uj\2-iE2 - ... = (121) 



-uPx-2E - uj^xi-2Ei + D2E2 = 



(122) 
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where Da = k'^c^ — u'^e + , ka = k + are the wave vectors of photons 
diffracted by the crystal planes with corresponding reciprocal vectors Tq,, eo = 
1 + Xo ; Xa cire the dielectric constants of a periodic structure; is the part 
of dielectric susceptibility appearing from the interaction between an electron 
beam propagating in the undulator and radiation. Setting the determinant of 
the linear system fll2ip equal to zero, one can obtain the dispersion equation 
for the system "electromagnetic wave + undulator + electron beam + periodic 
structure". In the case of two- wave dynamical diffraction, this equation has 
the following form: 

DD, - uj\,x-i = (123) 

For the system with finite interaction length, the solution of the boundary 
problem can be presented as a sum: 

E + El = ^ Cj exp{ikif}{e + eiSi '' exp irr) (124) 

i 

In (11241) Si '' = {kfc'^—uj'^eo)/{uj'^Xi) ^-iid ki are the coupling coefficients between 
the diffracted and transmitted waves and the solutions of dispersion equation 
(11231) . respectively. The coefficients Cj are defined by boundary conditions at 
the system ends z = and L. The part of the electron beam energy converting 
into radiation can be expressed by: 

/~7o|^(^ = L)\^ + hi\\Ei{z = 0)p (125) 
= (7o|a| +1711161 )(^j -mx.\L.nL.iT,,.„-Tr 
where L is the length of interaction in the undulator 
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f{y) is the function describing the dependence of generation on detuning from 
the synchronism condition, y = (to — kvu, — Vt)L/{2Uz) is the detuning factor, 
/3 = 70/71 is the diffraction asymmetry factor, 70, 71 are diffraction cosines, 
X" = Imxo, ©s = eH^/{mc^'jk^), 7"^ = 7-2 + 6^, = 2n/\^, is the 
undulator period, is the undulator field. It follows from fll25p that: 

1. the starting current in the case of two-wave diffraction is proportional to 

Jstart ~ {kL)-^{kxiL)-^] 

2. the non-one dimensional VDFB provides the possibility to decrease the 
starting current of generation by varying the angles between the waves. 

3. if the electron beam current is less than the starting value j < jstart? the 
energy of the electromagnetic wave at the system entrance can be written in 
the form: 

En,.\a? + = 1 - (^l^)V(y) (126) 

The conventional FEL gain is proportional to [kL)^ [104], but as follows from 
(11261) in the case of two-wave diffraction, the gain gets an additional factor 
~ (^4^) , which noticeably exceeds unity in conditions of dynamical diffrac- 
tion. Such an increase of the radiation output at the degeneration point can 
be explained by the reduction of the wave group velocity, which can be written 
as: 

"»--(i)/0-n(fe.-M (127) 
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It follows from fll27p that for multi-wave dynamical diffraction at the s-fold- 
degeneration point Vg ~ vq/ {kLy, the starting current jstart ~ {kL)~'^{kxiL)~'^' 
and the amplification are proportional to {kLY {kxiLY'^ . It should be noted 
that the considered effects take place in a wide spectral range for wavelengths 
from centimeters to X-ray, and the influence of the effect increases with fre- 
quency growth. 

The generation threshold in the undulator VFEL in the case of VDFB can be 
achieved at a lower electron beam current and a shorter undulator length when 
special conditions of degeneration of roots are fulfilled. Changing the VDFB 
conditions by varying the volume geometry parameters (for example, the angle 
between the wave vectors) gives the possibility to increase the Q-factor and to 
decrease the starting current (see Figure ??). Hence, the generation efficiency 
can be increased. 

16 Use of a dynamical undulator mechanism to produce short 
wavelength radiation in VFELs 

Here we shall consider VFEL lasing in a system with a dynamical undulator 
|15U] . In this system radiation of long wavelengths creates the undulator for 
lasing at a shorter wavelength. Two diffraction gratings with different spatial 
periods form a VFEL resonator. The grating with longer period pumps the 
resonator with the long wavelength radiation to provide the necessary ampli- 
tude of the undulator field. The grating with the shorter period is used to 
select the mode for the short wavelength radiation. Lasing of such a system 
in the terahertz frequency range is discussed below. 
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Numerous applications can benefit from tlie development of powerful elec- 
tromagnetic generators with frequency tuning in millimeter, sub-millimeter 
and terahertz wavelength range using low-relativistic electron beams. One of 
the ways to create such generators is to use VFEL principles. The main dis- 
tinction of VFEL in comparison with the traditional PEL is the use of the 
non one-dimensional distributed feedback, which allows wide range tuning of 
frequency, decreases starting currents of generation, and allows one to use a 
wide electron beam (or several beams) [55|105] . In the VFEL the generation 
evolves in a large volume, which increases the electrical strength of the res- 
onator (the electromagnetic power and electron beam are distributed over a 
greater volume). This peculiarity of the VFEL is essential for generation of 
power and super-power electromagnetic pulses. The mode discrimination in 
such an oversized system is carried out by multi-wave dynamical diffraction 
|105j . Low relativistic electron beams in the undulator system can be used 
for radiation of short wavelength radiation, but it requires manufacturing of 
undulators with a small period. For example, to obtain radiation with the 
wavelength of 0.3 — 1 mm at the beam energy E = 800 KeV— 1 MeV an 
undulator with the period ~ 0.3 — 1 cm is necessary. This is an extremely 
complicated problem. The use of a two-stage FEL with the dynamical wiggler 
generated by an electron beam |104j is a possible solution to this problem. 

The dynamical wiggler can be created with the help of any radiation mecha- 
nism: Cherenkov, Smith-Purcell, quasi- Cherenkov [55] , undulator. VFEL prin- 
ciples provide advantages of the two-stage generation scheme and, in partic- 
ular, allow one to smoothly tune the period of dynamical wiggler by rotating 
the diffraction grating. There is the possibility of smooth frequency tuning 
for both the pump wave and the signal wave either by variation of geometric 
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parameters of the volume diffraction grating or by their rotation. Moreover, 
the VFEL allows one to create the dynamical wiggler in a large volume, that 
is a great problem for a static wiggler. 



diffraction grating 




Figure 16. The schemes of dynamic wiggler. Rotation of diffraction grating (schemes 
a), and c).) changes the period of dynamic wiggler 

There are two stages in the generation scheme proposed above: 

(a) creation of the dynamical wiggler in a system with two-dimensional (three- 



108 



dimensional) gratings (in other words, during this stage, the electromagnetic 
field, which exists inside the VEFL resonator, is used to create the dynamical 
wiggler). Smooth variation of the orientation of the diffraction grating in the 
VEFL resonator provides means for a smooth change of the dynamical wiggler 
parameters; 

(b) radiation is generated by an electron beam interacting with the dynamical 
wiggler, which is created during the previous stage (stage a). 
Both stages evolve in the same volume. 

The idea of a two-stage tunable VFEL described above can be realized in 
different ways. Let us consider some examples (Figure [TBb-c). Figure [TB] a 
displays the two-stage device in which the dynamical wiggler is realized on 
the basis of a VEFL generator using the Smith-Purcell radiation mechanism. 
This generation scheme was considered in |107] , when experimentally realized, 
first lasing of a VFEL was observed |181j 

In this case only a part 6lx of an electron beam participates in the generation 
process during the first stage: ~ A„,/37/(47rZ^), here 51^ is the transverse 
size of the part of the electron beam participating in the interaction, is the 
transverse size of the electron beam, A^^, is the wavelength, /3 = u/c, u is the 
electron velocity, 7 is the Lorentz factor. The resonator in Figure fT6k consists 
of two diffraction gratings |164|181j . The lower grating provides the Smith- 
Purcell generation mechanism. The period di of this grating is determined by 
di ~ f3'y^\w cosip {di is the period of the Smith-Purcell grating, ip is the angle 
between the direction of the electron velocity and the direction of the grat- 
ing periodicity). The upper grating provides the distributed feedback [55|105j 
by multi-wave dynamical diffraction. The conditions of dynamical diffraction 
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\kw\ ~ l^u; + Til are fulfilled in this case (xj are the reciprocal wave (lattice) 
vectors of this grating). It should be noted that the period of the upper grat- 
ing does not coincide with that of the lower one. Radiation accumulated in 
the resonator during the described first stage creates the dynamical wiggler. 
Beam electrons oscillate in this electromagnetic wiggler and radiate just as 
in a conventional FEL (this is not necessarily the same electron beam that 
participates in the first stage. It can be another beam of higher energy). The 
field inside the resonator is a standing wave. Traveling waves, which form this 
standing wave are actually the pump waves. They are scattered by the electron 
beam according to the synchronism condition uj — (ku) ^ — (k^u). The 
resulting wave has the wavelength A ~ j^^^w (in the last estimation it is sup- 
posed that the wave vector ky, of the pump wave is antiparallel to the velocity, 
and the wave vector k is parallel to it. For a relativistic beam, this relation 
has the form A ~ A^/ (47^). It should be emphasized that in this case there is 
one more possible way to create the dynamical wiggler in the resonator shown 
in Figure [T6h . It is based on the excitation of a slow wave, which is diffracted 
by the lower grating (surface VFEL). The electron beam oscillates in this wig- 
gler and radiates. In this case the upper grating forms the distributed feedback 
which provides VFEL lasing at shorter wavelength. The change of the radiated 
frequency is provided by the rotation of both the upper and lower gratings. 
Figure [T6b presents the variant of the volume diffraction grating which can 
provide the generation mechanism and distributed feedback simultaneously 
[55|105j . Let us note that in these examples the generation mechanism during 
the first stage is based on the slowing of the electromagnetic wave and only a 
part of the electron beam participates in the creation of the dynamical wig- 
gler. The whole electron beam participates in the generation process during 
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the second stage. The larger the part of the beam that does not participate in 
the the first stage interaction, the more unperturbed electrons appear during 
the second stage and, therefore, they radiate more effectively. The dynamical 
wiggler in Figure [T6b uses the undulator radiation mechanism during the first 
stage. In this case the dynamical wiggler is formed during the first stage due 
to the interaction between the electron beam and a conventional magneto- 
static undulator. The frequency of the pump wave is lOw ~ ^ (i-^fos(6>)) ' where 
du is the period of the magnetostatic undulator. During this stage the lower 
diffraction grating is used to provide the distributed feedback for the wave 
with frequency and operation of the VFEL at this frequency. During the 
second stage, electrons oscillate in the field of this wave, which plays the role of 
the dynamical wiggler. As a result, during the second stage of the process the 
wave with the frequency u ~ 47^0;^ is generated and the distributed feedback 
is provided by the upper grating, the period of which corresponds to the wave 
with the frequency u. The rotation of diffraction gratings provides frequency 
tuning. 

It is clear that the time of the dynamical wiggler creation is a very important 
characteristic of the proposed system. This time is determined by ~ Q/wiu, 
where is the frequency of the pump wave. 

The Q factor of resonator for the frequency should be sufficient to create 
a magnetic field amplitude of about 100 G - 1 kG. It follows from the energy 
balance equation in the resonator, / Q)V {H^/ Sir) = Wq, {Wq is the power 
of the pump wave formed by the electron beam) that Q = {uw /Wo)V {H^/ Stc) , 
where V is the cavity volume, Hm is the amplitude of the magnetic field of 
the dynamical wiggler. It follows from the above that to create the magnetic 
field of about 100 — 1000 G, the time ~ 10"^° — 3 x 10"^ s is necessary 
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{V = 30 cm^, ly ~ 30 MW). As one can see, this time is small enough for the 
wiggler to evolve while the electron beam passes through the system. When 
the pump field achieves the necessary magnitude, stage (b) begins. 

Dynamics of the signal electromagnetic wave and the electron beam in the 
system (volume diffraction grating + pump electromagnetic wave) in this case 
is described by equations 

Do E -u'^xiEi - W^X2^2 - W^X3^3 - . . . = 
-uj\.^ E +D,E, - uj\2-iE2 - uj\3-iEs - . . . = 0, (128) 

-Lo'^X-2 E -L0'^Xl-2Ei + D2E2 - W^X3-2-£^3 - . . . = 0. 



In f ll28p Da = k^c^ —uj'^e + Xa\ ka = k + fa is the vector of the diffracted wave, 
Xa^ is the part of the dielectric susceptibility corresponding to the interaction 
of the electron beam with radiation 



- (Kvw) =F {(^w - (kwv))} 



(w) _ Q'w 



{k^,v) 



[ue. 



-{oJwU kyjC ){ka kyj^ 



k, u, e and k^,, Uw, are the wave vectors, frequencies and polarization vectors 
of both the signal and pump waves, respectively, V — (c, u), kyi] — (cu^/c, A;^), 
ttw = eHw/mcujw The dispersion equation corresponding to equation f ll28p 
has the following schematic form 



(129) 
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where the term F^^^ on left-hand side of (11291) corresponds to the n-wave Bragg 
dynamical diffraction (equation F^'^^ = is the dispersion equation defining 
diffraction modes in the n-wave case). The continuity of the current densities, 
charge densities, and transverse components of fields on the boundaries and 
dispersion equation fll29p give the equation for the generation threshold jl05] . 
From (11281) we obtain that for the n-fold degeneration point of the roots of the 
dispersion equation (when n + 1 roots of the equation D^™-) = at m > n + 1 
coincide), the equation for the generation threshold has the following form 



In (I130p k = oj/c, is the length of the interaction area of the electron 
beam with electromagnetic radiation, x" is the imaginary part of the dielectric 
susceptibility, which describes absorption, a„, 6„ are the parameters depending 
on the geometric parameters of the system (except L^). The equality ( llSOp 
has an obvious physical meaning: the left-hand side of (11301) contains the term 
describing generation of radiation by the electron beam, and the right-hand 
side includes the terms describing losses on the boundaries (the first term) and 
absorption losses (the second term) in the medium. One of the peculiarities of 
the VFEL with multi-wave distributed feedback is the possibility of a sharp 
decrease of losses at the boundaries (the first term on the right-hand side 
of (11301) decreases with the increase of s due to the condition ^ 1 

under dynamical diffraction). Let us express the synchronism condition for 
the above generation mechanism u — kv = Qyj, where = ~ vkw Then, 
the frequency of the signal wave is equal to (if the pump wave is oncoming) 

2a;^(fi, . . .f„,n„, 5')(1 - /3cos(e^)) 




(130) 



u = 



l-/3cos(e) 



(131) 
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In fll3ip the explicit dependence of the pump wave frequency on the geometry 
of the muhi-wave diffraction (ri, . . .7^) and the set of resonator parameters 
S is marked out (if a resonator is not oversized, then the dependence on S 
disappears). Smooth change of the VFEL geometry also varies the Q factor 
and, therefore, varies the generation efficiency. For example the dependence of 
Q on the diffraction asymmetry factor / = 70/71 is shown in Figure [T7] (70, 71 
are diffraction cosines j55]). 




-20 -15 -10 -5 



Figure 17. Calculated dependence of Q factor on diffraction asymmetry factor / 

It should be noted that the distributed feedback can be used for both the first 
and the second stages. To optimize the resulting radiation output, Q factor 
can be controlled at both stages. Thus, for radiation angle = the frequency 
u ~ 47^0;^, even the moderately relativistic electron beam {E ^ 1 MeV) gives 
a frequency multiplication ~ 35 times. If during the first stage the undulator 
mechanism is used (undulator period ~ 8 cm), then the wavelength of the 
pump wave is ~ 1 cm. Thus the signal wave is generated in the teraHertz 
range (Figure [T7|) . 

Thus, 

(1) the principles of VFEL can be used for creation of a dynamical wiggler 
with variable period in a large volume. 
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(2) the two-stage scheme of generation can be used for lasing in the teraHertz 
frequency range using low-relativistic beams, 

(3) the two-stage scheme of generation combined with the volume distributed 
feedback opens up the possibihty of creating powerful generators with wide 
electron beams (or system of beams). 

17 Formation of distributed feedback in a FEL under multi-wave 
diffraction 

We shall further consider the boundary problem for a quasi- Cherenkov FEL 
for the case in which a distributed feedback is provided by three-wave copla- 
nar diffraction of the emitted photons |105j . As compared with a two-wave 
distributed feedback, in the three-wave case the appropriate choice of param- 
eters enables reducing the interaction region between a particle beam and the 
electromagnetic wave necessary to reach the oscillation threshold as well as ra- 
diation losses inside the medium. For example, it allows a reduction in the size 
of an FEL and the construction of more compact coherent radiation sources 
in any spectral range. 

Here we show that the two-wave DFB (in which only two waves are strongly 
excited) is not optimal because the region of root degeneration often coincides 
with the region of strong radiation absorption inside the medium. This is es- 
pecially important for solid FELs. The most effective resonator is one with 
a multi-wave DFB, where the distributed feedback is formed by multi-wave 
dynamical diffraction. The advantages of a multi-wave DFB are analyzed in 
detail for the solid quasi- Cherenkov FEL with three-wave coplanar DFB be- 
cause this analysis can be conducted analytically. Consider, for example, the 
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specific case of a solid X-ray quasi- Clierenkov FEL wliere tlie crystal medium 
provides both a spontaneous radiation mechanism in the X-ray spectral range 
[H] and diffraction of the emitted X-rays by the crystal forming the three- 
dimensional distributed feedback. For the X-ray spectral range, the crystal 
target is a radiator and a resonator simultaneously. It should be mentioned 
that the analysis derived below will be appropriate for other spectral ranges 
as well. For example, a three-dimensional optical grating can be formed inside 
a solid target by a laser. Moreover, a multi-wave distributed feedback can be 
formed by surface dynamical diffraction if a particle beam moves over a three- 
dimensional periodic structure at a distance not larger than (where A is 
the radiation wavelength, 7 is the Lorentz factor). In this way the multi-wave 
DFB can be used with an ordinary undulator FEL. 

Let a relativistic particle beam be incident on a crystal plate {0 < z < L) 
at an arbitrary angle \l/o. The set of Maxwell's equations which describes the 
interaction of an electromagnetic wave with a crystal and with a particle beam 
passing through a crystal target can be written in the following form |55] : 

-co\.,E(^^ + D(J^E(J^ - c^^-ii??) = 0, (132) 
-u\.,E^J>^ - u\,.,E^^^ + Df^Ef^ = 0, 

where D'^f^ = k'^c^ — ui'^e'^^ + xi""*- We assume that a particle beam and a 
crystal plate are oriented so that the three- wave coplanar diffraction condition 
is satisfied for emitted photons. In this case only three strong waves with the 
cr-polarization are excited inside the crystal medium, (see [55] for two-wave 
diffraction geometry). The subscript a (a = — 2 ) labels the transmitted 
wave {a = 0) and diffracted waves (a = 1,2); i?^"^ are the a components 
of the amplitudes of electromagnetic waves, ki = ko + fi, k2 = ko + T2 are 
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the wave vectors of photons diffracted by crystal planes with corresponding 
reciprocal vectors ti and T2] e^f^ = 1 + Xa are the dielectric constants of a 
crystal for transmitted and diffracted waves. 



for the "hot" beam limit. Xb represents the part of the dielectric susceptibility 
produced by the interaction of a particle beam with radiation, = (w — 



Ku)/y/25o., 51 = (kl^^l + kl^^l + kl^^UDu^, and ^ = AV/\V\ is the velocity 



spread. As was shown in [55], comparison of a standard equation of X-ray 
dynamical diffraction with f ll32p leads to the conclusion that the combination 
of a crystal and a particle beam may be considered as an "active" medium with 
dielectric susceptibility equal to Xa + xi""* • It permits the boundary problem 
of X-ray amplification (lasing) due to the passage of a particle beam through 
a periodic medium to be reduced to the problem of X-ray diffraction by an 
"active" periodic medium consisting of a crystal plus radiating particle beam. 

The geometry of three-wave diffraction is shown in Figure [181 where V is the 
mean particle beam velocity, z = and z = L are the two surfaces of a crystal 
plate. 

The dispersion equation determining the electromagnetic modes inside the 
"active" medium can be represented in the following form: 




1 



- 00^ 




for the "cold" beam limit and 





(134) 
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Figure 18. Figure The geometry of three-wave diffraction. V is the mean particle 
beam velocity, and z = and z = L are two surfaces of a crystal plate 

where 



4/7 2 2 2 \ 4/7 2 2 2 

-u [k c -u eo)xi-2X2-i - w [k^c -u ')x2X-2 
-uj'^{klc^ - uj'^ei}'^)xiX-i - t^^(XiX-2X2-i + X2X-1X1-2); 

(fci; k,) = {ky - ujh'P){kl^ - ujh^^^) - uj\,_,X2-i 



From (11341) follows that the root degeneration (the strongest interaction is in 
the root degeneration region) and the fulfillment of the Cherenkov condition 
are possible simultaneously only under the following conditions: 

khh - kri2 - kr2 - kri - / = 0, 

^Mk + Wih - P2I2) - /3i/32ri2 - Ari - /^sra = 0, (135) 
where /q = 6*^ - Xo + 7"^ ^1,2 = + "m,2; "51,2 = (2/cfi,2 + 2)/^. 
It is straightforward to show that for the fulfillment of (I135p . the system 
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parameters should satisfy the following relationship: 



(136) 



where / = X1X-2X2-1 + X2X-1X1-2, Pi,2 = 7o/7i,2 are the asymmetry factors 
of diffraction ■ja = {kon)/\ka\, ri = XiX'-i aiid r2 = x'2X'-2- Expression (11361) 
is more restrictive than the equivalent relation for the two-wave case [55] . 

In the case of Bragg-Bragg diffraction in the X-ray spectral range there is a 
restriction on the radiation angle 

e^^B = {^J-f3in~P2r2 + Xo-rY'- (137) 

This relation leads to a large value of the diffraction asymmetry factor, which, 
in turn, leads to strong radiation absorption inside the medium in the vicinity 
of the root degeneration point - just as in the two-wave DFB case [53]. The 
situation changes in the Bragg-Laue geometry. In this case the inequality f ll36p 
can be satisfied at angles 6 > 6b-b, which makes it possible to reduce the 
asymmetry factors /3i and P2 and, consequently, radiation absorption inside 
the crystal. 

The solution of the corresponding boundary problem is presented as a sum: 
E = ^^Ci exp ikif{eo + eis^/'^ exp ifir + €28^^ exp ir2r), (138) 

i 

where s^^^ = (AA2 - ra)/ (A2X1 + X2X1-2), s^^) = (AAi - ri)/ (A1X2 + X1X2-1), 
and ki are the solutions of the dispersion equation and the coefficients of 
coupling between the transmitted and diffracted waves, 

s^^-'E' , Aq, = [{k + r)^c^ — a;^eo]/w^. To determine the unknown coefficients. 
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it is necessary to solve the boundary conditions for the waves on the crystal 
surfaces. These can be written in the following form for the Bragg-Laue case: 

3 3 3 

^ C, = 1, = 0' E ^C,/, = 0. (139) 

i i i 

The condition (I139p is written for the "hot" beam limit. For the "cold" beam 
limit the corresponding expression can be found in 



It is well known |116] that the oscillation threshold can be determined from 
the condition that A = 0, where A is the determinant of the system (I139p . 

Solving the equation A = 0, we obtain the threshold in the form 



G = axo + 



16 r7ncl^ vr^ri^ 



|/3i/32 



nu 



{kl^L^ykL, 



-Vb-l (140) 



with the phase condition {ki^ — k2z)L = 2nn {n is an integer), ri2 = X1-2X2-I5 
where 



G = {^}' k'Lli^, ± V^lx.l - 7-')(Xo ± V^IX.I) sin, 

X sm.y[{2y + Tin) siny — y{y + vm) cosy]y^{y + vm)"^ 



for the "cold" beam limit and 

r- (^0 + Xo - 7~^)(^o + X'o) sin 0^ ^-(^^)2 

for the "hot" beam limit, a and rjBL are smooth functions depending on the 
diffraction geometry and usually are of the order unity, y = k6L/2 and 6 is 
the deviation from the exact Cherenkov synchronism condition. 

The analysis shows that in the coplanar Bragg-Laue diffraction geometry of 
DFB radiation, absorption inside the medium can be reduced in the region of 
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root degeneration. For example, under dynamical diffraction by the (111) and 
(111) planes with a symmetry factor /3i = — /32 = P = 0.16, the current density 
required to achieve the threshold can be reduced by approximately one order 
of magnitude in comparison with a two-wave DFB PEL using LiH and an 
electron energy of 750 MeV with a transverse angular spread \E'x = 5 x 10~^ 
rad [55] . 



18 Distributed feedback under the multi-wave diffraction 



In principle, in the general case of the n-wave diffraction DFB, it is possible 
to achieve the degeneration of n roots. In this case the threshold condition 
has the form: 



where a„ and 6„ are the functions weakly dependent on the diffraction type. 

The energy loss through target surfaces can be reduced due to the fulfillment 
of the inequality klxrlL ^ 1 under dynamical diffraction. As the degree of 
dispersion root degeneration increases, radiation remains inside the target 
longer because the group velocity decreases and is proportional to ~ 
where n is the degree of dispersion root degeneration. As a result, the region 
of interaction between the particle beam and the emitted radiation can be 
reduced. This leads to the possibility of reducing the FEL's dimensions, i.e., 
of constructing a compact source of coherent radiation. 
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19 Parametric X-ray FEL operating with external Bragg reflectors 

This section is devoted to tlie application of external Bragg reflectors in a 
volume parametric X-ray FEL |151] . It is shown that external Bragg reflectors 
permit not only a reduction of electromagnetic energy losses through target 
boundaries but also a reduction of the radiation self-absorption inside the 
target by modifying radiation modes excited in the active medium containing 
a particle beam plus a crystal. As a result, the starting beam current can be 
reduced by more than a factor of ten. It is also shown that the best conditions 
for lasing are realized when the diffracted wave is reflected by the external 
Bragg reflector. 

19.1 Generation threshold for parametric X-ray FEL with external reflectors 

The scheme of an X-ray volume FEL using spontaneous X-ray parametric 
radiation has been considered in [55] . Following the remarks mentioned above, 
the scheme under consideration differs from the case described in ^55j by the 
addition of external Bragg reflectors which can be placed in the direction of 
both transmitted and diffracted waves. For example, in the case of a three- 
crystal scheme of diffracted wave reflection, the geometry of a volume FEL 
can be represented in the following way: 

— * — # 

Here V is the mean velocity of electrons, k and to are the wave vector and the 
frequency of the emitted photon, k and u are the wave vector and the frequency 
of the diffracted photon, a, b, c, d are the crystal targets. We assume that the 
crystal (a) in FigJT9] is oriented relative to the particle beam so that only 
two strong waves are excited inside the crystal under diffraction (two-wave 
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diffraction). 




Figure 19. 



The closed system of Maxwell's equations, including the equation describing 
the motion of electrons, and the method of its solution were described in [55] . 

The introduction of external Bragg reflections modifies only the boundary con- 
ditions where the reflection coefficients for diffracted and transmitted waves 
appear. Besides, the analysis shows that there appears the possibility of real- 
ization of lasing regime under a Laue diffraction scheme (see Figure [T9ll which 
is impossible, in principle, without external Bragg reflectors. 

Let us restrict ourselves to the analysis of Laue geometry for distributed feed- 
back. Solving the boundary conditions while taking into account the phase 
relation for the main crystal target (a) and for the external reflectors we can 
obtain the generation equations for two cases. 

Case (1) The "cold" beam limit and weak amplification: 
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G 



where 



(1 - |«|)(1 - f )^ + k[x'l + f = T^£, (1) 

(1 - |«|)(i - f + k[x'i + ^x'i]u = rg), (2) 



(142) 



fix) = (143) 



/[, = 1/2 [-as ± ^al + 4r], /; = l/2[+aB ± V«l + 4r], 

^1(2) - ^XO |i + ± j(;, _ )2 + 4^^]l/2 j ' 

the prime means the real part of the magnitude and the double prime means 
its imaginary part; equation (1) of the system fll42p describes the case of 
diffracted wave reflection, and equation (2) corresponds to the case of trans- 
mitted wave reflection. The quantity a = |a| expi(j) is the reflection coefficient 
for transmitted (2) or for diffracted (1) waves respectively. 

The value of (1 — |a;|) can be very small for a narrow interval of angles (1 — |q;|) = 
— Xo|Xr|~^(l — exp(— ty)), where W is the Debye- Waller factor. For example, 
for a LiH crystal (1 — |a|) < 10^'^ at ~ 1 — 10 keV photoenergies. 

As a result, the radiation losses through crystal boundaries can be essentially 
reduced. This means that the threshold condition will be determined mainly 
by radiation absorption. In this case (11421) can be represented in the form: 



(1 ^th = -^ A2 + ^Ai — —— = ——- ^c/i a,/3 , 

fn\ ■ ^ ( \" , ^0 \"\ bo + f3bi m-fc^ x'o , 

(2) Jth = 77 Ai + — -A2 TTs — — n^g2{aJ3), 



(144) 
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here 6o(i) = Ta+ [a^ + l]V^ a = ^, A'/(,) = ^, = a? = ^. 

The dependence of gi and (72 on the parameter of deviation from the exact 
Bragg condition is shown in Figure [20] (a) and (b) for a LiH crystal, the (111) 
diffraction plane, E = 700 MeV, P = 1, u = 3 x 10^^^ c^^ and the angular 
divergence of electrons in the beam A\l/x = 10~^ rad. 




T T - ■ '■ T 



Figure 20. 

From Figure 1201 (a) it can be seen that g2{0') has a minimum of ~ 200 at 
a definite deviation parameter a = —2. At the same time, (71(a) does not 
have a definite optimal value. It depends weakly on the asymmetry factor and 
approaches 2 as the absolute value of the deviation parameter a increases. 

A detailed comparative analysis of these figures shows that under the same 
conditions, the threshold current in the first case is 100 times larger than in 
the second case. It follows from Figure |20] (b) and fll44p that 

This can be explained by the fact that the growth rate of the electromagnetic 
field is maximum at the deviation parameter 0:^/2 |xt-| ~ —10. The analysis 
shows that in this case the main part of radiated electromagnetic energy turns 
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out to be accumulated in the diffracted wave. This means that for achieving the 
best condition for lasing, the Bragg reflector should be placed in the direction 
of the diffracted wave. 



Case(2). The "hot" beam limit. 

The boundary condition in this case is simpler, and for the diffracted wave 
reflection we can obtain the dispersion equation in the form: 



th ' 



(145) 



where 



7 % 

= , /o = e^-xo-7-^ 5o = uQA^^, (146) 

V2do 

O is the radiation angle, A'^± is the transverse angular divergence of electron 
velocities in the beam. 



At a small value of (1 — |a|), when the threshold condition is determined by the 
radiation absorption, the relation analogous to fll44p for the "cold" particle 
beam can be represented as: jth ~ ^dsio,, (3), where a is the function similar to 
that in (11441) . It weakly depends on the asymmetry factor, the particle energy, 
and the crystal parameters. The dependence of on the deviation parameter 
a is shown in Figure! 



We can see that in this case there is an optimal value of at which the 
magnitude of is minimum. At the same parameters as in the case of a 
"cold" particle beam and at A\&_l = 5 • 10"*^ rad jth = 2-10^ A/cm^. 
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Figure 21. 

So, from the above analysis we can conclude that the application of exter- 
nal Bragg reflectors allows one to lower the threshold electron beam current 
density of the electron beam necessary for achieving lasing in the X-ray spec- 
tral range by more than a factor of ten. The best situation is realized when 
the Bragg reflector is placed in the direction of the diffracted wave. In this 
case the value of the starting beam current can be jth — 2 ■ 10^ A/c for the 
"hot" particle beam and jth — 10*^ A/c for the "cold" particle beam and the 
generation frequency a; ~ 3 ■ 10^^ c~^ in LiH, E = 700 MeV, /3 = 1 and 
< 5 ■ 10^3 



20 Theory of induced PXR in a crystal plate (general formulas) 

Let an ultra-relativistic electron (positron) beam of velocity u enter, at a cer- 
tain angle, a crystal plane-parallel plate with length L (the 2;-axis is perpen- 
dicular to the crystal surface, and the plate lies in the interval {0 < z < L). 
The set of equations describing the interaction of an electromagnetic wave 
with the "crystal-beam" system consists of Maxwell's equations and those of 
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particle motion in the electromagnetic field. The dielectric susceptibility of a 
crystal has the form e{f;uj) = X^r £^t(i^) exp(— irr), where r is the reciprocal 
lattice vector. Perturbations of the current density and charge density in the 
linear field approximation may be written in the form: 

6j{k; cj) = e ^ exp{—ikfao) ^6va{oo — ku) 

a 

—iu[k6fa{uj — ku)]^6n{k-,uj) (147) 
= e exp(— iMrao) ^—i[kfa{ui — ku)]^ , 

a 

where Sj{k] u) and 6n{k] u) are the Fourier transformations of the expressions 
j{r;t) = eY,Vo.{t)6[f-f^{t)] 

a 

n(f;t) = ^5[f- rait)]. 

a 

u is the unperturbed electron (positron) velocity; 5va and Sr^ are perturba- 
tions of the velocity and radius vectors, respectively, due to the interaction 
with the radiation field: 

Va{t) = U + 5Va{t) 

ra(t) = faO + Ut + 5fa{t). 

The subscript a denotes the number of the particle. 

Deriving the perturbations of velocity and radius vectors by the particle mo- 
tion equation and using fll47p for the current density, one can obtain a set 
of Maxwell's equations, which describes the interaction of an electromagnetic 
wave with a crystal, and a particle beam penetrating through it, in the fol- 
lowing form: 

klE{kri,U]) — kri[kr'E{kri,Uj)] ^ Erikr' , Uj) E (kr+T' , Uj) 
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E{kT-i,uj) 




(148) 



X {uE{kr',uj)) 



( 



f = 0, Ti,T2, . . .. 

The set of equations fll48p describes the situation when a distributed feedback 
is formed by many diffracted waves (it is the analogy of the case of multi- 
wave X-ray diffraction in crystals [3]). However, the analysis of such a general 
situation is very complicated. So we only consider here a two- wave distributed 
feedback. This allows us to obtain all the main characteristics of X-ray FELs 
analytically and also to show the advantages of three-dimensional geometry 
of distributed feedback in comparison with the one-dimensional case. 

So, let us consider specifically the generation of a cr-polarized wave, i.e., the 
wave polarized in the diffraction plane. For the geometry of the so-called two- 
beam diffraction |190j . where two strong waves are excited and diffraction 
occurs by the set of crystallographic planes, determined by a reciprocal lattice 
vector. In this case one can obtain a set of Maxwell's equations describing 
two-wave diffraction in a crystal, having a beam penetrating through it, in 
the following way: 



In f lTl9|) E„ = E{k,u) ■ 4, E^ = E{k + T,u) ■ 4, e || [kr], w| = Ane^no/m, 



where no is the average electron (positron) density in a beam. Comparing 




e: = o. 



(149) 
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f ll49p with the standard equation of X-ray dynamical diffraction, one can see 
that the system of a crystal and particle beam may be considered as an active 
medium with dielectric susceptibility: 



eo{k,uj) - 1 = £0 - 1 



7^^ 7^^ c2 {io-kuY' 



eo{K,uj) -1 = Sq-I 



ul ul {ue^f klc^ - 



7U;2 ^^2 ^2 (^^_^^^^2' 



Further, we shall analyze the generation of the wave with a wave vector fc, 
which makes a small angle with the particle velocity vector u. In this case 
the wave vector k^. = k + t is directed at a large angle relative to u, and 
consequently the magnitude of {u — krU) cannot become small. As a result, the 
terms containing the expression {u — kru) in their denominators will be small 
and can be ignored. We shall also neglect the term ul/'y - this is justified 
for real beam densities. It is well known that in order to provide non-zero 
solutions for the equation set (11491) . its determinant should be equal to zero. 
This also defines the dispersion equation, and for the cr-polarized wave it can 
be written in the form: 



The dispersion equation in such a form was derived in [51]. To solve the bound- 
ary problem, we use field continuity, the beam density and the beam current 
density at the boundaries. For the last two conditions we apply the following 
expressions, obtained from the equations of particle motion and the expression 
for the particle beam current 



{u - kuf Uk'^c^ 




(150) 
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= e{ue„)no -^E„- (151) 

"^7 c^{uj-ku) 



The dispersion equation (11501) is sixth-order and, hence six solutions corre- 
spond to it. However, the two solutions corresponding to mirror reflected 
waves can be neglected due to the small value of er{uj) in the X-ray range, 
Eq — 1 = go; Er ~ 10^^. Small values of Er and Qq will also be utilized when 
performing joining (we only join electric field strengths at boundaries). 

The general solution for the field in a crystal is written as 

E = '^e^Ciexp{ikir)[l + Siexp(iff)], (152) 

i=l 

where ki is the i-th solution of the dispersion equation (11501) and kir = ki + r; 
T is the reciprocal lattice vector corresponding to the planes of diffraction 
reflection. 

Taking into account the remarks made, the boundary conditions may be writ- 
ten as: 

Cl + C2 + Cs + C4 = 1, 

/ici + /2C2 + /3C3 + /4C4 = 0, 

gici + 9202 + 9303 + 9iCA = 0, (153) 
sicie'^'^^ + sacae*^^^^ + s^c^e'^^^^ + SiC^e'^^'^ = 0, 

fc^c - (w - kiu) {uj - kiuf c2 

In (I153p . only those boundary conditions are written that determine the field 
inside a crystal. The first equation corresponds to the continuity of an incident 
wave at the boundary z = 0; the second and the third conditions correspond 
to vanishing of the beam density and the beam current density at the crystal 
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entrance. The last condition corresponds to the vanishing of the diffracted 
wave at the exit boundary z = L (we consider Bragg diffraction geometry), 
ki and ki^ (2 = 1-^4) are the solutions of the dispersion equation fllSOp . The 
linear system fll53p . defining the coefficients q placed ahead of z-modes in 
(11521) . has the solution q = Aj/A, where A is the determinant of the system 
(I153p . Aj is the i-th minor, obtained as a result of replacement of the z-th 
column by 



1 



Hence, at A — )■ 0, the field amplitudes inside a crystal will increase, and thus 
a condition occurs when the field is nonzero though the incident wave is equal 
to zero. The condition A = with Aj 7^ is called the generation threshold 
condition [116j . Substituting the expressions 



ki = kQ + kSifi, 
k = u/c, Si <^ 1, 



oj — k±u± 



(154) 



(where n is the normal to the crystal surface, /cq = (/eqz, ^±)) into the deter- 
minant A, we can represent the generation threshold condition A = as 



(^1 -^2)('^l ^ jkSiL _ i^l - - Sa){62 - 64) ^ ikSsL n c-c-N 

616161 ''^ 6j6l6j ''^ 



{6i-63){6i-6,){63-64) 



616161 



-S2e 



ik52 



L {^2 - 6^){62 - ^4) (1^3 - ^4) ^ ikS^L 



O2O3O4 



sie 
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In (11551) . the terms containing nonresonance fi and g-i (i = 1 -i- 4) were ne- 
glected. 

Further investigation will be based on a common consideration of (11501) and 
(11551) . We substitute (11541) into (I150p and transform the dispersion equation 
as follows: 



-(5^[47o7i5' + 2(7i/ + 7o/.)(5 + 11^ - r] 



C2 



■-(-) 

7 V w / 



2 

e'^sm'^ip{l + go)lr. (156) 



In (I156p 6 is the angle between k and u; (p is the angle between and tl^, 
where the symbol J.^ denotes the projection onto a plane perpendicular to the 
velocity, r = erS-r] I = 6'^ + go + 7~^; It = I + a', a = {2koT + T'^)/k'^ is the 
departure from the Bragg condition. 

(kofi) {korfi) 



7o = —r^, 7i 



K ' K 

are the cosines of the angles, made by the wave vectors of the transmitted and 
diffracted waves with the normal vector. 



21 Spontaneous and induced parametric and Smith-Purcell radia- 
tion from electrons moving in a photonic crystal built from the 
metal threads ( [117j ) 

Research and development of microwave generators using radiation from an 
electron beam in a periodic slow-wave circuit (traveling wave tubes, backward 
wave oscillators, etc.,) has a long history [69|70] . First generators operated in 
the centimeter wavelength range. 
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In 1953 Smith and Purcell [7T] made the next step and observed generation of 
incoherent radiation at visible wavelengths by using a finely-focused electron 
beam propagating close to the surface of a metal diffraction grating (at the 
distance 6 < 6 is the so-called beam impact parameter, A is the radia- 
tion wavelength, /3 = t>/c , f is the electron beam velocity, 7 is the electron 
Lorentz-factor) . Beam current densities were insufficient to produce signifi- 
cant amplification of the spontaneous emission. However, as it was shown in 
[72Ff3] . to change this spontaneous radiation process into a stimulated one, 
the generator should be supplied with a feedback cavity formed by a pair of 
refiectors. 

After the discovery of the Smith-Purcell effect, it soon became clear that it 
might be used as a radiation source in the millimeter to visible range, for 
which tunable sources were hardly or not available at that time [71|75] . 

The Smith-Purcell effect belongs to the general class of diffracted radiation 
effects induced by electron interaction with a medium. Diffraction of waves 
associated with the electromagnetic field of the electron by an obstacle leads 
to the so-called diffracted radiation |76j. Diffracted radiation in periodical 
structures is in the basis of operation of traveling wave tubes [69|70] and such 
devices as the orotron [77|78f79] and the ledatron [821 (see also [8T|82f83|l8^ ). 

All the above devices use feedback, which is formed by either two parallel 
mirrors placed on both sides of the working area or a one- dimensional diffrac- 
tion grating, in which incident and diffracted (reflected) waves move along the 
electron beam (one-dimensional distributed feedback). 

The conception of volume distributed feedback |19] was flrst originated in view 
of prediction and experimental study of spontaneous parametric and diffracted 
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transition X-ray radiation from charged particles in crystals (PXR) [85]. The 
detailed analysis of the induced PXR and channeling radiation revealed unique 
possibilities provided by volume distributed feedback |49J. It was shown that 
dynamical diffraction in a volume spatially periodic medium evokes the pecu- 
liar conditions, corresponding to degeneration of the eigenmodes that results 
in a new law of electron beam instability. Within these peculiar conditions the 
electron beam interacts with the electromagnetic wave more effectively |19]. 

Thus, for example, even in the X-ray range the electron beam current den- 
sity necessary for running up to the generation threshold in conditions of 
non-one-dimensional distributed feedback [19|25] appears significantly reduced 
(10^ A/cm^ for LiH crystal against 10^^ A/cm^ required in [HS]) that even 
makes possible to reach generation threshold for the induced X-ray radiation 
in crystals i.e. to create an X-ray laser. 

Moreover, all the conclusions are valid for a beam moving in vacuum close to 
the surface of the periodic medium either in the presence or absence of the 
undulator [25] . 

The originated law is universal and valid for all wavelength ranges (from X-ray 
to microwave) regardless the spontaneous radiation mechanism 
[25|49f55f87f88fl05fll0fl81fl82] . Radiation generators using non-one-dimensional 
distributed feedback, which is created with the aid of either natural or artificial 
(photonic) crystals, is called Volume Free Electron Laser (VFEL). 

Use of the volume distributed feedback makes available: 

1. frequency tuning at fixed energy of the electron beam in the significantly 
wider range than conventional systems can provide; 
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2. significant reduction of tlie tliresfiold current of the electron beam due to 
more effective interaction of tfie electron beam with the electromagnetic wave 
allows and, result, miniaturization of generators; 

3. reduction of limits for the available output power by the use of wide electron 
beams and diffraction gratings (photonic crystals) of large volumes; 

4. simultaneous generation at several frequencies; 

5. effective modes selection in oversize systems, in which the radiation wave- 
length is significantly smaller than the resonator dimensions. 

Studies of the two-dimensional distributed feedback application in millimeter- 
wave FEL-oscillators began in the 90s in p^95f96f97f98] . 




One of the VFEL types uses a "grid" volume resonator ("grid" photonic crys- 
tal) that is formed by a periodic structure built from either dielectric or metal 
threads (see Figure 1221) [5^11Ufll2j . 




\ from metallic threads 
electron beam 

Figure 22. VFEL "grid" resonator ("grid" photonic crystal) 

The "grid" structure formed by periodically strained dielectric threads was ex- 
perimentally studied in |110j . where it was shown that "grid" photonic crystals 
have sufficiently high Q factors (10^ — 10^). Lasing from VFELs with the "grid" 
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resonator formed by periodically strained metal threads was observed in [112j . 

Propagation of waves through a photonic crystal is the subject of numerous 
theoretical and experimental studies pU|91f92|l^ . 

In |117] the properties of a "grid" photonic crystal built from metal threads 
were considered along with its frequency characteristics in view of their im- 
portance for VFEL lasing. A challenge, which appears when considering inter- 
action of an electromagnetic wave with such a photonic crystal, is as follows. 
In contrast to the case of wave interaction with a thin dielectric thread, the in- 
teraction of the electromagnetic wave with a single metal thread of the "grid" 
(i.e., the diffraction grating unit cell) cannot be described in terms of the 
perturbation theory. 

The approach developed in |117] provides the description of diffraction in a 
"grid" grating in terms of the amplitude of the wave scattering by a single 
thread. Methods for calculation of scattering amplitudes are well-developed (TUI 
This approach enables developing the theory of diffraction in a " grid" photonic 
crystal similar to the dynamical theory of diffraction for X-rays |3|IT5] and us- 
ing the results obtained therein. The equations describing lasing of VFEL with 
such a resonator have been obtained. 



22 Scattering by a set of metal threads in free space 

Suppose a plane electromagnetic wave E = "^qc = e^^'^e falls onto a metal 
thread. Here e is the polarization vector of the wave, k is the wavevector, r 
is the coordinate. Suppose the wave falls perpendicular to the cylinder axis 
i.e. along OZ. The metal thread looks like a cylinder placed into the origin of 
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coordinates and the cylinder axis coincides with the axis x Figure (123|) . 




e H 

Figure 23. Coordinates and vectors 



Two orientations of e should be considered: e is parallel to the cylinder axis x 
and e is perpendicular to the cylinder axis x. The solution of the problem of 
the wave scattered by a cylinder is well-known |100lll01j . Below the case when 
the radiation wavelength A exceeds the thread radius R {X^ R) is considered. 

For clarity, suppose that e || Ox. In this case the expression, which describes 
the wave appearing due to scattering, can be expressed as a superposition of 
the incident \l/o and scattered "^sc waves as follows |100fl01j : 

^ = ^o + ^sc = e''^ + aoi/i'^(M, (157) 



here p is the transverse coordinate p = {y,z), Hq^^ is the Hankel function of 
zero order, k is the wave number, ao is the amplitude of the scattered wave. 



_ -MhR)j[,{kR) + ^tJoihR)MkR) 

''"^"^ UktR)Hi^\kR) - ^MktR)H^^\kR) ' 



-MktR)J'^{kR) + ^^J'^{ktR)UkR) 
~ MktR)H!,'^'{kR) - ^^J^{hR)Hi^\kR) ' ^^^^^ 



where St is the dielectric permittivity of the thread material, kt = y/slk. Note 
that for A ^ i? the amplitude does not depend on the scattering angle. 
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Let us now consider a wave, which is scattered by a set of threads placed 
at the points with the coordinates p„ = (?/„, The scattered wave can be 
expressed as a superposition of waves scattered by separate threads: 

^ = e^*^^ + ao^nH!,'\k \p-p^\)e'''- (160) 
or, using the integral representation for the Hankel function, H'^\kp) 



ik^\p-p-af-x'^ 

^ = e'^^^ + /loS„ / , =dxe'^'-, (161) 



oo 

e' 



-oo "4 \P - Pnf - 



where = |p - p^f = {y - Vnf + {z - z^f. 

Note that in (116 ip we neglected the contributions from rescattered waves. 

Let us consider a wave passing through a layer of cylinders whose axes are 
distributed in the plane xOy with the distance dy between them. Suppose 
that the transversal size of the layer L_l is much larger than both dy and 
the radiation wavelength (L_|_ ^ dy and ^ A). This assumption enables 
considering the ideal case when the layer is supposed to have an infinite size 
in the plane xOy. In this case summation over the coordinates yn provides the 
following expression for \E^: 

^ = + ^^e^^^^ = (1 + ^^y^\ (162) 

This expression reflects a well-known fact that a plane wave scattered by a 
plane layer of the scatterers is expressed as a plane wave with the modified 
amplitude. 

Thus, after passing m planes spaced dz apart from each other, the scattered 
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wave can be expressed as: 



\ V kdy ) \ kdy 



where 

/ 27r ReAp 

ip = arctg 



\ k dy 



m = Lz/dz inside the photonic crystal formed by threads and is the length 
of the photonic crystal. This expression can be easily converted to the form 
^ _ gifenz -yyiiere n is the index refraction defined as 



( \ ( i^R-6Ao 

n = n + in" = 1 + arctan ' — - 



STToiz I i_ >^ImAo 

"■y , 



-^2^1n(\/(tRe^o)" + (l-ilmAo 



(164) 



here A = 2-^ jk is used, v! and n" denotes the real and imaginary parts of n, 
respectively, RcAq and ImAg are the real and imaginary parts of Aq. 

Let us consider now the parameter |^^o| and suppose it to be small. When 
I^AqI ^ 1, fll64p can be rewritten as: 

n = l + ^^Ao. (165) 



The same expression can also be obtained for the index of refraction of the 
wave with polarization orthogonal to the thread axis (e _L Ox). 

When the parameter |^^o| grows with the growth of the density of scatterers 
(i.e., with decreasing d^^ the rescattered waves become important and the 
difference between the mean and local fields should be considered similarly to 
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Clausius-Mossoti (Lorentz-Lorenz) relation [99]. Calculations show that the 
parameter |^^o| behaves differently for different polarizations (see Figure [241) 
and for the wave with polarization parallel to the threads, this parameter can 
appear greater than 1. In this case formula (11651) for parallel polarization is 
not valid and the expressions considering wave rescattering should be applied 
(see [TT7] ). 
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Figure 24. Dependence of parameter l-j-^ol on kR for parallel and orthogonal po- 
larizations {dy = 0.5 cm, dz = 1.6 cm) 



It is important to remind that the conception of the refraction index is appli- 
cable even to considering waves of different types (X-rays, particles) passing 
through matter when the distance between scatterers exceeds many times the 
wavelength, i.e., kd ^ 1 pLj. At that, if the wavelength is either comparable or 



141 



smaller than the size of the scatterer, then the refraction index is determined 
by the amplitude of the wave scattering at zero angle. In the considered case 
of the wave scattering by the thread the amplitude Aq in fll65p should be re- 
placed by the amplitude of forward scattering of the wave by the thread. The 
explicit expression for this amplitude see in |100|101] . 

Radiation wavelengths of our interest are A < 3 cm. In this wavelength range, 
the skin depth 6 is about 1 micron for most metals (for example, 6cu = 0.66 
fim, 6ai = 0.8 /im, 6w = 1.16 /xm, and so on). Thus, in this wavelength range 
the metal threads can be considered as perfectly conducting. 

From fllSSp . f ll59p follows that the amplitude Aq for a perfect conducting 
cylinder for polarization of the electromagnetic wave parallel to the cylinder 
axis can be expressed as: 



_ 1 Jo (kR) No (kR) . 1 Jg (kR) 

TT J2 (^kR) + Ni (kR) ^ \ J2 (kR) + (kR) ■ ^ ^ 



Using the asymptotic values for these functions for kR = <^ 1 



Jo(x^O) ^ 1, iVo(x^O) ^ -Mn 



TT 1.781-2;' 



(167) 



4 {x^0) = -Ji ^ -|, iV^ (a; -> 0) = -N^ ^ \\, 
one can obtain: 



1 llnf ^ \ 1 1 



From fll68p it can be seen that for kR ^ 1, the values ReAodi) < 0, therefore 
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the real part of the index of refraction in this hmit can be expressed as: 

^(11) = 1- ^^^1^^6^0(11)1, (169) 



From fll69p follows that the growth of the wavelength (reduction of k) makes 
nj| — 0, and nj| can even become negative. However, it should be reminded that 
(11691) is valid only when \n'^^ — 1| ^ 1 (i.e ^ ^"^^2 |Re^o(||) I ^ !)• It is interesting 
to compare nj| calculated by formulas (I164p and (1169^ : in the range of large 
wavelength they yield different results because the condition — 1| ^ 1 is 
violated. In this case the expressions given in |117] should be applied. Note 
also that njj = corresponds to the threshold of the wave total reflection from 
the medium surface. 

Let us consider a particular example: suppose the radiation wavelength is 
A = 3 cm and the thread radius R = 50 fim, then 

i?eAo(||) = -0.096, /mAo(ii) = 0.032, (170) 

Suppose dy = 0.5 cm, = 1.6 cm, therefore, the parameter | j-^o| ^ 1 for 
both polarizations and ( 11651) gives: 

nil = 0.828 + 2-0.058, (171) 

The above analysis does not include the contribution to the index of refraction 
from rescattered waves. Their contribution was investigated in 
|89|103|112|117|152fl65j . It was shown, in particular, that the the index of 
refraction n± for a wave with the polarization vector orthogonal to the threads 
can be greater than unity, i.e., in such photonic crystals, the Cherenkov effect 
is possible. 
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23 Complex anomalous Doppler effect in the "grid" photonic crys- 
tal 



Let us now study Smith - Purcell (diffracted) radiation in a photonic crystal 
when an electron beam of velocity v passes through the "grid". 

The radiation condition can be expressed as 

uj — kn{k)v = TV , (172) 



where r is the reciprocal lattice vector and n{k) is the index of refraction. 
Suppose the electron beam velocity is directed along the axis OZ, then (11721) 
can be presented in the form: 

k-T^(3 = kn{k)/3 cos 9, (173) 



where /3 = |, the angle between k and the electron beam velocity is denoted 
by 9 and = ?2n!ih^ where rrih = 1, 2, ... is the harmonic number. From (11731) 
follows the equation 

i^-f^y = k' + V. (174) 
p cost^ 

which is similar to the equation for the complex and anomalous Doppler effect 

The roots of this equation give the spectrum of frequencies of diffracted 
(Smith-Purcell) radiation, which is induced by a particle moving in the above 
volume "grid" structure. In the case under consideration <^ 1. For /3 pro- 
viding <^ 1, the roots of (I174p can be easily found: 



144 



= '^^^ I l-/3cos9 



7] 1 - {(3 COS oy 



1 - (/3C0S^)2 y ' ^ T-2 13 



(175) 



1 + /3cos9, 



' 1- (/3cos^)2 \^ ^2 ^ 

where r? is taken at k = , ,1''^ . 

It should be reminded here that = , where m/j = 1, 2, ... is the harmonic 
number. From f ll76p follows that higher harmonics provide getting radiation 
with higher frequencies. For example, for the electron beam with the energy 
200 keV, considering 6 ~ 20° and dz = 1.6 cm, the first harmonic {mh = 1) 
gives radiation frequencies ~ 10 GHz and ~ 40 GHz for the roots 1 and 2 of 
equation f ll73p . respectively, the 30-th harmonic (m/j = 30) provides ~ 230 
GHz and ~ 1 THz. 

Let us now study diffracted radiation in a metal waveguide of rectangular 
cross-section with the "grid" structure ("grid" photonic crystal) placed inside 
it (see Figure l22il . Suppose that the lateral dimensions a and b of the waveguide 
noticeably exceed the distance between the next threads of the "grid" photonic 
crystal {a,b ^ dy,dz). This makes possible to consider a waveguide with 
the "grid" photonic crystal inside as a waveguide filled with matter with the 
effective dielectric permittivity e = n^. 

Remember that in the waveguide with the axis parallel to OZ, the x and y 
components of the wavevector are not continuous, but quantized. And only 
the wavenumber k^ changes continuously. The wave field inside the waveguide 
represents a wave standing in the transversal directions {OX and OY), but 
traveling along OZ. 
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In the waveguide there is the discrete set of the waves with the eigenval- 
ues Kmn determined by the waveguide transverse dimensions (width a and 
height b) 



mn \ 



,7m, 2 



(176) 



Therefore, in the waveguide the radiation frequency u is related to the wavenum- 
follows flGl: 



kl{m,n) = {-fe-^l^. (177) 
c 



Remembering that = e = 1 + x (e is the permittivity of the "medium" 
formed by the metal threads and x is its susceptibility) one can rewrite f ll77p 
as follows: 



kl{m,n) = iy-iKl,-v) , (178) 



here r] = ^X- According to the previous section n < 1, therefore both x and 
r] are negative. 

Therefore, addition of the described volume structure to the waveguide "re- 
places" their own eigenvalues by some effective K^n (-^mn — '^mn ~ v)- It is 
interesting to note here that, as r] < 0, the addition of the volume grid to 
the waveguide increases the waveguide limiting frequency (i.e., the effective 
waveguide appears as though having smaller transversal dimensions). 

Let us now consider the spontaneous Smith-Purcell radiation from a parti- 
cle moving along the waveguide axis. The radiation condition for the Smith- 
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Purcell radiation fll72p converts to 

UJ — kzV = TzV . 



(179) 



Combination of fll78p and fll79p enables one to find the equation for radiation 
frequencies similar to fll74p : 



(180) 



The roots of this equation for <C 1 can be obtained similarly fll76p : 



Ui (m, n) 



UJ2 {m,n) 



. _ (^mn - r]) 1 - 



(181) 



(/^mn - ^?) 1 - 



/32 



For a detailed treatment of the theory of VFEL lasing using electron beam 
radiation in a "grid" photonic crystal, see |117] . 



24 Generation of radiation in Free Electron Lasers with photonic 
crystals (diffraction gratings) with the variable spatial period 

In this section the equations providing the description of the generation pro- 
cess in PEL with varied parameters of a photonic crystal (diffraction grating) 
are obtained |152j . It is shown that applying photonic crystals (diffraction 
gratings) with the variable period, one can significantly increase the radiation 
output. It is mentioned that photonic crystals (diffraction gratings) can be 
used for creation of the dynamical wiggler with variable period in the system. 
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This makes possible to develop double-cascaded FEL with variable parame- 
ters changing, which efficiency can be significantly higher that of conventional 
system. 

Generators using radiation from an electron beam in a periodic slow-wave 
circuit (traveling wave tubes, backward wave oscillators, free electron lasers) 
are now widespread [Tllj . 

Diffracted radiation [70] in periodical structures is in the basis of operation of 
traveling wave tubes (TWT) [69f70] . backward wave oscillators (BWO) and 
such devices as Smith-Purcell lasers [7T|72f73f74f75] and volume FELs using 
two- or three-dimensional distributed feedback [571I88II181II182I. 



The analysis shows that during the operation of such devices electrons lose 
their energy for radiation, therefore, the electron beam slows down and gets 
out of synchronism with the radiating wave. These limits the efficiency of the 
generator, which usually does not exceed ~ 10%. 

During the first years after the creation of the traveling wave tube, it was 
demonstrated [70] that synchronism between the electron beam and the elec- 
tromagnetic wave in a TWT can be retained by changing the phase velocity 
of the wave. Application of systems with variable parameters in microwave 
devices significantly increases their efficiency [70|114] . 

The same methods are widely used for increasing the efficiency of undulator 
FELs [TT5] . 
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24 -i Lasing equations for the system with a photonic crystal ( diffraction grat- 
ing) with changing parameters 



In the general case the equations, which describe lasing process, follow from 
the Maxwell equations: 

IdD Att^ IdH 

rotH = -— H J, TotE = --, 

c dt c c at 

divD = Airp, ^ + divj = 0, (182) 

— * — * — * 

here E and H are the electric and magnetic fields, j and p are the current 
and charge densities, the electromagnetic induction Di{f,t') = J eii{f,t — 
t')Ei{r,t')dt' and, therefore, Di{r,uj) = eii{r,uj)Ei{f,uj), the indices i,l = 
1, 2, 3 correspond to the axes x, y, z, respectively. 

The current and charge densities are respectively defined as: 

J(f,t) = eJ2Ut)6{f-f^{t)), p{f,t) = eJ26{r- r^{t)), (183) 

a a 

where e is the electron charge, Va is the velocity of the particle a {a numerates 
the beam particles), 

dt [ c c"^ J 

2 1 — * — * 

here 7q, = (1 — ^)~^ is the Lorentz-factor, E{fa{t),t) {H{fa{t),t)) is the 
electric (magnetic) field at the point of location of the particle a. 

It should be reminded that (11841) can also be written as [187] : 

dpa djaVa f f^/ - /^N , 1 i 



dt 



= e [E{Ut),t) + -iUt) X ^(r.(t),t)]} , (185) 
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where pa is the particle momentum. 
Combining the equations in fll82p . we obtain: 

The dielectric permittivity tensor can be expressed as e{r) = 1 + x(r^, where 
x(r) is the dielectric susceptibility. When x ^ 1, fll86p can be rewritten as: 

AE(f,t) - 1^ j i{r,t- t')E{r,t')dt' = 4n (^^^ + Vp(f,t)^ .(187) 

When the grating is ideal x(0 = Z^t XTi.r)^^'^'^ , where r is the reciprocal lattice 
vector. 

Let the photonic crystal (diffraction grating) period be smoothly varied with 
distance, which is much greater then the diffraction grating (ptotonic crystal 
lattice) period. It is convenient in this case to present the susceptibility x{r) 
in the form, typical of the theory of X-ray diffraction in crystals with lattice 
distortion |113j : 

x(r) = Ee^*^^'^X.(r1, (188) 

where $r(0 = / T{f')dl', T{f') is the reciprocal lattice vector in the vicinity of 
the point r '. In contrast to the theory of X-rays diffraction, in the case under 
consideration Xt depends on f. Moreover, depends on the volume of the 
lattice unit cell fl, which can be significantly varied for diffraction gratings 
(photonic crystals), as distinct from natural crystals. The volume of the unit 
cell Q{r) depends on coordinate and, for example, for a cubic lattice it is 
determined as il{r) = (if)d2{f)d3{r) ' ^^ere di are the lattice periods. If Xr(^) 
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does not depend on r, the expression fllSSp converts to that usually used for 
X-rays in crystals with lattice distortion |113j . 

It should be reminded that for an ideal crystal without lattice distortions, the 
wave, which propagates in the crystal can be presented as a superposition of 
plane waves: 

oo 

E{f, t) = J2 ir-e*^'^^"^-"*) , (189) 

r=0 

where kr = k + f. 

Let us now use the fact that in the case under consideration the typical length 
for the change of the lattice parameters significantly exceeds the lattice pe- 
riod. Then the field inside the crystal with lattice distortion can be expressed 
similarly to f ll89p . but with depending on r and t and changing noticeably 
at the distances much greater than the lattice period. 

Similarly, the wave vector should be considered as a slowly changing function 
of a coordinate. 

According to the above, let us find the solution of (11871) in the form: 

E{f, t) = Re jf] l^e'^-^-^"^-"*)! , (190) 

where 0t-(^) = lo k{r)dl + ^ri"^), where k{r) can be found as a solution of the 
dispersion equation in the vicinity of the point with the coordinate vector r, 
integration is made over the quasiclassical trajectory, which describes motion 
of the wavepacket in the crystal with lattice distortion. 

Now let us consider the case when all the waves participating in the diffraction 



151 



process lie in a plane (coupled wave diffraction, multiple- wave diffraction), i.e., 
all the reciprocal lattice vectors r lie in one plane [31IT5] . Suppose the wave 
polarization vector is orthogonal to the plane of diffraction. 

Let us rewrite fllQOp in the form 

i(f, t) = eE{r, t) = e Re {i^e'^'^i^"^-"*) + ^e^C-^^Cr)— t) + ^^^| ^ (ig^) 

where 

r 

= J Hr')dl 



r r 

02(r) = j h{r')dl + J T{f')dL 



Then multiplying fll87p by e, one can get: 

AE(f,t)-i^ J e(f,t-t0i5^(f,O^^t' = 4vre(^l^^ + Vp(f,t)j .(194) 

Applying the equahty AE{f,t) = V(V-E) and using (11911) . we obtain 

A(iie'('^^("^)-"*)) = e*(^i("^-"*)[2iV0iVAi + iVA?i(f)Ai - kl{r)Ai], (195) 



(192) 
(193) 



Therefore, substitution of the above expression into (I194p gives the following 
system: 



^e'^Mr)-u^t)^2ik^{r)^Ai+iVki{r)Ai - kl{f)Ai 

+-eo{u:, r)A, + ^^^^^ + ^^-.(-, r)A, 

. 1 du'^E-riuj, r) dA2, 
+1 



dio dt 

1 d](f^t) 
dt 



+ conjugated terms = ine i — — h ^p{r, t) 
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^ 1 duj'^eQ{(jj,r) dA2 uJ^ 



^. 1 9a;^£:r(w,f) SAi^ 



du dt 

1 (9j(r^t) 
9t 



+ conjugated terms = 47re — — ^ h Vp{f,t) ) , (196) 



where vector k2{r) = ki{f) + T, eo{ijj,r) = H-Xo(^)) here the notation Xo(^) = 
Xr=o{f^ is used, er{uj,f) = Xr(^)- Note here that for a numerical analysis of 
(11961) . if xo ^ 0, it is convenient to take vector ki{r) in the form ki{r) = 



n#' + ^Xo(r1. 

For better understanding, let us suppose that the diffraction grating (photonic 
crystal lattice) period changes along one direction and define this direction as 
axis z. 

Thus, for a one-dimensional case, when k{{r)) = {k±, kz{z)), the system (11961) 
converts to the following: 

le^(.^x.^.+0..(.)-.t)[2,fci.(;^)^ + ^^^^1 - (kl + kUz))A 

, , , .1 dco^eoico, z) dAi , , ^ 

+-.o(c., z)A, + ^-2^^^ + -le-A^. AM 

1 du?e_A^, z) (9^2, 



(? du dt 

1 d](f^t) 
dt 



+ conjugated terms = 47re I — — h Vp{r, t) ) , 



l^.(.^,.-,+0.4.)-t)[2,fc2,(^)^ + ^^^^2 - {kl + kUz))A2 

, / 1 duj^£o{uj, z) dA2 uj^ , . 

+-eo{., z)A2 + ^-2^^^ + -^eAu:, z)A, 

^ . 1 du'^Eriui, z) dAi ^ 



(? doj dt 

1 d]{r^t) 
dt 



+ conjugated terms = Anel — — h Vp{r,t) ) , (197) 
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Let us multiply the first equation by Q-^i^±-^^+'t>i-Az)-^^t) ^y^q second by 
^-t(krr^+<f>2z{z)-u)t) _ This procedure enables neglecting the conjugated terms, 
which appear fast oscillating (when averaging over the oscillation period they 
become zero). 

Considering the right-hand side of f ll97p . let us take into account that micro- 
scopic currents and densities are the sums of terms, containing delta-functions, 
therefore, the right-hand side can be rewritten as: 



c2 dt ^i^s) 



C2 



eJ2Utmr) - (r),(t))e-*('=^^^^+'^-(^)^'^*) e{t - t^) 9{T^ - t) 



here ta is the time of entrance of particle a to the resonator, is the time of 
particle leaving the resonator, 6'— functions in f ll99p indicate that for the time 
moments preceding ta and following T^, the particle a does not contribute to 
the process. 

Let us suppose now that a strong magnetic field is applied for beam guid- 
ing through the generation area. Thus, the problem appears one-dimensional 
(components and Vy are suppressed). Averaging the right-hand side of f ll99p 
over the particle positions inside the beam, points of particle entrance to the 
resonator r^oa and time of particle entrance to the resonator ta one can obtain: 



e-^i^^r^+Mz)-<^t)^^^( ]_d3i^ ^ Vp(r,t) 

at ^ 

A7riup-&iu{t) e 1 /" ,2^ ^ 



c2 SJ 'T 







2 --^u,,, (199) 
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where p is the electron beam density, u{t) is the mean electron beam velocity, 



which depends on time due to energy losses, t?! = y 1 ~ 1^'^ — ^ ~ 

(( )) indicates averaging over the transversal coordinate of the point of particle 
entrance to the resonator r_|_OQ, and time of particle entrance to the resonator ta- 

According to |102j . the averaging procedure in (11991) can be simplified, when 
consider that random phases, appearing due to random transversal coordi- 
nate and time of entrance, presents in (11991) as differences. Therefore, double 
integration over (fr±o dto can be replaced by single integration |102] . 

The system (11971) in this case converts to: 

, ^dAi dkiz(z) , ,,9 ,9 / XX , 
2ikr,{z)^ + ^^^^1 - {kl + kUz))A^ 

, . , .1 du'^Eoiu, z) dAi u? , . , 

+— £:o(w, z)Ax + %— -— + —e.^[u, z)A2 

c c OiO cit c 



1 du'^e^ri^, z) dAo 2uj ^ , , ^ 



2zk,z{z)^ + ^^^A, - {kl + kl{z))A, 

^2 1 du^eoiuj, z) dA2 . . 

+-eo{u., z)A, + ^-1^^^ + ^eAu:, z)A, 

. 1 duj'^er{uj,z)dAi .2uj 
where the currents Ji, J2 are determined by the expression 



J„ = 27Tj^m f ^L^(e-*<^™(*'^'P) + e-i<^-(t,^,-P)) dp, m = 1, 2 
0^2 







(201) 



u r 



2 



j = eriQV is the current density, Ai = A^-^q, A2 = A^-, ki = A;^=o, ^2 = ^1 + ^• 
The expressions for Ji and ki independent on z was obtained in |102j . 
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When more than two waves participate in the diffraction process, the system 
(120 ip should be supplemented with equations for waves Am, which are similar 
to those for Ai and A2. 

Now we can find the equation for the phase. From (11921) . (11931) follows that 



(P(f>m 1 d(j)m dkm km d'^ Z 

dz^ V dz dz dz dt"^ ' 



(202) 



Let us introduce a new function C{z) as follows: 



— — = Cm{z)e Jo -d.' = -^Cm{z), 

dz v{z) 

z 

——Cm{z')dz' 
^ ' 



(203) 



Therefore, 



dCmiz) v(z) f dkm kmd'^z^ 



dz 



+ 



vq \ dz ' dt^ 



(204) 



In a one-dimensional case, equation (I185p can be written as: 



d'^Zn 



dt'^ m^{za,t,p) 
therefore. 



ReE{zo,, t), 



(205) 



dCmiz) v(z) dkm. ^ k 



e-d 



dz Vq dz Vov{z) m'y^{z,t{z),p) 

d(j)mit,z,p) 



Re{Am{z, t(z))e^*'"("'*(^)'^')}, (206) 



dz 



z=0 = kmz - -, 4>m{t,Z,p)\z=o = P, 
V 



A,l=L = El A 



2\z=L — ^21 



Am\t=o = 0, m = 1,2, 
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t > 0, z ^ [0, L], p E [—2tt, 27r], L is the length of the photonic crystaL 



These equations should be supplied with the equations for 7(2;, p). It is well- 
known that 



mc'^ = evE. (207) 
dt ^ ' 



Therefore, 



d'j{z,t{z),p) ^ ei? 



J2^Re{Y,Mz,tiz))e'^'^'''^'^^P^}. (20J 



dz ^ mc 

The above obtained equations (l2nTl) . (IM . (12061) . (IM) enable describing 
the generation process in a FEL with varied parameters of diffraction grat- 
ing (photonic crystal). The analysis of the system fl206p can be simplified by 
replacement of the 'y{z,t{z),p) with its averaged by the initial phase value 

< 7iz,t{z)) >= — y -f{z,t{z),p)dp. 


Note that the law of parameters change can be both smooth and stair-step. 

Using photonic crystals provide the development of different VFEL arrange- 
ments (see Figure |25|) . 



It should be noted that, for example, in the FEL (TWT, BWO) resonator 
with changing in space parameters of grating (photonic crystal), the electro- 
magnetic wave with a spatial period depending on z is formed. This means 
that the dynamical undulator with a period depending on z appears along the 
whole resonator length, i.e., a tapering dynamical wiggler becomes settled. It 
is well known that a tapering wiggler can significantly increase the efficiency 
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Figure 25. An example of photonic crystal with the thread arrangement providing 
multi-wave volume distributed feedback. Threads are arranged to couple several 
waves (three, four, six and so on), which appear due to diffraction in such a structure, 
in both the vertical and horizontal planes. The electronic beam takes the whole 
volume of photonic crystal. 

of the undulator FEL. The dynamical wiggler with varied period, which is 
proposed, can be used for development of double-cascaded FEL with param- 
eters changing in space. The efficiency of such a system can be significantly 
higher that of a conventional system. Moreover, the period of a dynamical 
wiggler can be made much shorter than that available for wigglers using static 
magnetic fields. 

It should also be noted that the compression of the radiation pulse in such 
a system is possible because the phase velocity of the electromagnetic wave 
depends on time. 

Thus, the equations providing the description of the generation process in 
FEL with varied parameters of diffraction grating (photonic crystal) are ob- 
tained. It is shown that applying diffraction gratings (photonic crystal) with 
the variable period one can significantly increase radiation output. It is men- 
tioned that diffraction gratings (photonic crystal) can be used for creation of 
the dynamical wiggler with variable period in the system. This makes possible 
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to develop double-cascaded FEL with variable parameters changing, which 
efficiency can be significantly higher that of conventional system. 



In conclusion, it should be noted that in the photonic crystals built from metal 
(or dielectric) threads, waves are effectively diffracted even when their lengths 
are much smaller than the distance between the threads. That is why it is 
possible to generate radiation in the terahertz range using diffraction gratings 
with a millimeter period alone. 

Above, we have studied the features of generation for the sources with trans- 
verse dimensions much greater than the radiation wavelength. In the case 
when the finite dimensions of a resonator should be taken into account, it is 
necessary to follow the rules given in [S7|117fl52fl63] . 

24-2 Radiative instability of a relativistic electron beam moving in a finite- 
dimensional photonic (or natural) crystal 

According to the analysis given in [163] , the first and most important step in 
describing the generation process in VFELs (FELs and so on) is the analysis of 
the problem of the electron beam instability in the resonator. The theoretical 
study of the instability of electron beams moving in natural and artificial 
(photonic) crystals was carried out above for the ideal case of an infinite 
medium (see the review in [87] and |49|102|117|152|181] ). The question arising 
in this regard is how the finite dimensions of the photonic crystal placed 
inside the resonator affect the law of electron beam instability. It is known, for 
example, that the discrete structure of the modes in waveguides and resonators 
is crucial for effective generation in the microwave range |83|l84|ll83yi84|ll85] . 
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In present section the radiative instability of a beam moving in a photonic 
crystal is studied. The dispersion equation describing instability in this case 
is obtained. It is shown that the law T ~ p^/('^+'^) is also valid and caused by 
the mixing of the electromagnetic field modes in the finite volume due to the 
periodic disturbance from the photonic crystal. 



The system of equations describing generation of induced radiation in pho- 
tonic (and natural) crystals can be obtained similarly to those in section 124.11 
(equations ffTg^ - dWj) ). 



In the general case, the susceptibility of the photonic crystal reads x (^) = 
S Xceii if — Ti), where Xceii if — ri) is the susceptibility of the crystal unit cell. 

i 

The susceptibihty of an infinite perfect crystal x if) can be expanded into the 
Fourier series as follows: x (^) = Z^Xr^*^^, where f is the reciprocal lattice 

T 

vector of the crystal. 



To be more specific, let us consider in details a practically important case 
when a photonic crystal is placed inside a smooth waveguide of rectangular 
cross-section. 



The eigenfunctions and the eigenvalues of such a waveguide are well-studied 



Suppose the 2;-axis to be directed along the waveguide axis. Make the Fourier 
transform of (11861) over time and longitudinal coordinate z . Expanding thus 
obtained equation for the field E (rl, fc^, w) over a full set of vector eigenfunc- 
tions of a rectangular waveguide Y^^ (rj_, k^) (where m,n = 1, 2, 3...., while A 
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describes the type of the wave [B], one can obtain for the field E the equahty 



E (f^, u) = Y: CL ^) yL {r±, h) . (209) 



As a result, the following equation can be written 

{k1 + ti^nx) ~ ^ ^mn i^z, ^) — 

E / {r±, K) X (r) (r-1, K) e'^^^^-^)^ d^^C^^, (K, u) dk'.dz = (210) 

'n'X' 



m'ra'A' 



where K^nA — ^ImX + ^^nA- 



The beam current and density appearing on the right-hand side of fl21Up are 
complicated functions of the field E. To study the problem of the system 
instability, it is sufficient to consider the system in the approximation linear 
over perturbation, i.e., one can expand the expressions for j and p over the 
field amplitude E and abridge oneself with the linear approximation. 

As a result, a closed system of equations comes out. For further consideration, 
one should obtain the expressions for the corrections 5j and 5p due to beam 
perturbation by the field. Considering the Fourier transforms of the current 
density and the beam charge j {k, u?j and p {k^ u?j , one can obtain from f ll83p 
that 

Sj (k, u;)=eJ2 S^c. (oo - ku^ + ^ , (211) 

Q=i [ w - ku^ J 

where r^g is the original coordinate of the electron, is the unperturbed 
velocity of the electron. 
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For simplicity, let us consider a cold beam, for which Ua ~ m, where u is the 
mean velocity of the beam. The general case of a hot beam is obtained by 
averaging 6j (jt, u?j over the velocity distribution in the beam. 

According to f ll84p . the velocity 6va is determined by the field E (r^, u) taken 
at the electron location point r„. The Fourier transform of the field E {faju) 
has a form 



As a result, the formula for 6j (k,uj) includes the sum J^^ ^)ra ^^^^ ^^^^ 
particle distribution in the beam. Suppose that the electrons in an unperturbed 
beam are uniformly distributed over the area occupied by the beam. Therefore 



^ e-<fc-fc') = (2vr)3 po6{k- k') 



where po is the beam density (the number of electrons per 1 cm^). 

As a result, the following expression for 5j {k,u?j can be obtained [5T|60] : 



6j(k,oj) = - ^ ^ " ^ uEik.uj). (212) 



{uj — ku 



Using the continuity equation, one immediately obtains the expression for 
p (^k, uij . Expression ( I212p . the inverse Fourier transform of E (jc, u?j , and the 
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expansion fl209p enable writing the system of equations fl210p as follows: 



-^l E I Y^n {r±, K) X (r) (r-1, K) e^C^-'^^)^ dh-^C^,^, (K, u) dk'^dz = (213) 



m'ra'A' 



where (^x, A;,) = / e-*^^^^^^f^„ (rl, A;,) d\^. 

Note that within the limit where the transverse dimensions of a photonic 
crystal tend to infinity, the expression between the braces takes the form 
(ew)^, where e is the unit polarization vector of the wave emitted by the beam 

Now let us consider the integrals on the left-hand side of equation f|213p . 
Note that according to [6|16fl88] . the eigenfunctions Y^^{r_!^,kz) of a rect- 
angular waveguide include the combinations of sines and cosines of the form 
sin^x, cos^x {sin^y, cos^y), i.e., in fact, the combinations e*^^,e*^^. 
Hence, the left-hand side of the equation includes the integrals of the type 



i 



Xceii (x -Xi,y- yi, z - Zi) e*"" ""dx. 



The substitution of variables x — Xi = rj gives the sums of the form 

i 

where Xi = d^fi, where d^ is the period of the photonic crystal along the 
X-axis, /i = 1, 2, ...Nx, where A^^; is the number of cells along the x-axis. 
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The above-mentioned sum 

^ ■K{m—m')dxNx 

sm 



^ .Z— / • TT(m—m')d.r. ' ^ ' 



2a 



If m — m' = 0, then 5*2; = A^^.. 



Let us now discuss what this sum is equal to when m — m' = 1. In the 
numerator d^N^ = a, hence, the nominator is equal to 1 {sin^ = 1), while in 
the denominator sin^ ^ tttt- As a result, the relation ^^(™— ™- -i) = 1 ^ g 

With growing difference m — m', the contribution to the sum of the next terms 
diminishes until the following equality is fulfilled 

TT (m — m') dx 



2a 



ttP, (215) 



where P = ±1, ±2... In this case the sum Sx = N^. 

The similar reasoning is valid for summation along the y-axis. 

It follows from the aforesaid that if the equalities like (12 14^ . ( 1215^ are fulfilled, 
that is, the equalities k^m — Km' — '^x are fulfilled (i.e., k'^^, = k^m — 'Tx), where 
= is the x-component of the reciprocal lattice vector of the photonic 
crystal, F = 0, ±1, ±2... and kyn - ky^, = Ty (i.., ky^, = kyn - Ty), where Ty = 
^F' is the ^/-component of the reciprocal lattice vector of the photonic crystal. 



y 



F' = 0, ±1,±2...), then the major contribution to the sums comes from the 
amplitudes C^',„, (A;!,, u) = C^' {k±mn - T_i_, k^ - T^,u?j = C^' (km 



r, u 



In describing the system we shall further consider only those modes that satisfy 
the equalities of the type (12 14p . (12151) . As stated above, the contribution of 
other modes is suppressed. 
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As a result, one can rewrite the system of equations (12130 as 



27r 



(216) 



I.e., 



k^ - ^ 

mn (p- 



1 + X^n(0) 



27r 



A'r ^ 



r, 



where Xmn i^) = tJ ^mn ir±, k,) x (r^, t,)Y^,^, (rl, /c. - r,) dV^, x r,) = 
S / Xcei/ — Xi,y — Hi, e~*'^^^(iC, m' and ra' are found by the conditions hke 



(I215p . is the Langmuir frequency, 



47re po 



This system of equations coincides in form with that describing the instabihty 
of a beam passing through an infinite crystal [5T][60] . The difference between 
them is that the coefficients appearing in these equations are defined differ- 
ently and that in the case of an infinite crystal, the wave vectors k^nn have a 
continuous spectrum of eigenvalues rather than a discrete one. 

These equations enable one to define the dependence k (w), thus defining the 
expressions for the waves propagating in the crystal. By matching the incident 
wave packet and the set of waves propagating inside the photonic crystal using 
the boundary conditions, one can obtain the explicit expression describing the 
solution of the considered equations. 

The result obtained is formally analogous to that given in [55] . 



(217) 



According to (12171) . the expression between the square brackets acts as the 
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dielectric permittivity e of the crystal under the conditions when diffraction 
can be neglected: 



eo = n' = l + xr^ (0) 



n is the refractive index. 



2n 



u 



As is seen, in this case the contribution to the refractive index comes not only 
from the scattering of waves by the unit cell of the crystal lattice, but also 
from the scattering of waves by the beam electrons (the term proportional to 
u1): the photonic crystal penetrated by a beam of electrons is a medium that 
can be described by a ceratin refractive index n (or the dielectric permittivity 

According to fl217p . the beam contribution increases when u — )■ ku. 

Since this system of equations is homogeneous, its solvability condition is the 
vanishing of the system determinant. 

In the beginning, let us assume that the diffraction conditions are not fulfilled. 
Then the amplitudes of diffracted waves are small. In this case the sum over 
r can be dropped, and the conditions for the occurrence of the wave in the 
system is obtained by the requirement that the expression between the square 
brackets equal zero. 

This expression can be written in the form (the velocity mHo^;) 



(u-k.uf(kl -i(e.c^--^)[ 1 



where rio is the refractive index of the photonic crystal in the absence of the 



uYmn [k±,kA d kj 
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beam = = 1 + x^^„ (0), 



I.e., 




'mn 



) 



) 



{uj — kzuf' 




I 



Since the nonlinearity is insignificant, let us consider as the zero approximation 
the spectrum of the waves of equation (12181) with zero right-hand side. 

Let us concern with the case when u — k^u — )■ (i.e., the Cherenkov radiation 
condition can be fulfilled) and (^k'^ — (^^nl — — )■ 0, i.e, the electromag- 
netic wave can propagate in a photonic crystal without the beam. With zero 
right-hand side the equation reads 



Since fcoz = - > in view of the Cherenkov condition, we are concerned with 
the propagation of the wave with kiz > in the photonic crystal. In this case 
in the equation for k^, one can take [k^ — ki^) {kz + kiz) ~ 2kiz {kz — kiz) and 
rewrite equation (I218p as follows: 




(219) 



As a consequence, in this case the roots of the equation are 




(220) 



{kz - kiz) (kz - k2zf 



ulu^ (ng - 1) [ 1 
2kizu'^-fc^ \ 27T 



I 



uY^n{k±,kz)d'k^ (221) 



I.e., 



{kz - kiz) (kz - k2zf = -A 



(222) 
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where A is real and A > (as for the occurrence of the Cherenkov effect, it 
is necessary that nl > 1). We have obtained the cubic equation for k^. Let us 
consider the case when the roots ki^ and k2z coincide kiz = ^22- It is possible 
when the particle velocity satisfies the condition 




Introduction of ^ = k — k^ gives for kiz = k2z 

= -A. (224) 



The solution of equation (12241) gives three roots ,^1 = —\/~A, ^2,3 = | (l i ^Voj v^- 

As a consequence, the state corresponding to the root ^2 = | (l + iV^ ) y/A 
grows with growing z, which indicates the presence of instability in a beam 
[189] . In this case Im kz = Im C,2 ~ i^P- 

Note here that the photonic crystal built from metallic threads has the re- 
fractive index uq < 1 for a wave with the electric polarizability parallel to the 
threads, i.e., in this case the Cherenkov instability of the beam does not exist 
[87] (but if the electric vector of the wave is orthogonal to the metallic threads, 
the refractive index is no > 1 , so for such a wave the Cherenkov instability 
exists |165] ). 

It should be pointed out, however, that, unlike an infinite photonic crystal, the 
field in the crystal placed into the waveguide has a mode character, and so the 
presence of in the denominator of equation (I223p results in reduction of 
the radicand in (12231) to the magnitude smaller than unity even when tIq > 1. 
Hence, u > c, which is impossible. Consequently, the radiative instability of 
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2 2 

the above type in the waveguide can arise under the condition — '^'^j^ > 1 
rather than Uq > 1. 



Suppose now that in the photonic crystal the conditions can be reahzed under 
which the wave amphtude Cmn (kmn + t^) is comparable with the amplitude 
Cmn (kmn^- By analogy with the standard diffraction theory for an infinite 
crystal |165p5] . in the case under consideration, when x << 1, it is sufficient 
that only the equations for these amplitudes remain in ( ]217p . 

To be specific, let us further consider a photonic crystal formed by parallel 
threads. Also assume that they are parallel to the waveguide boundary {y, z). 

Analysis of diffraction of a A-type wave with the electric vector in the plane 
(y, z) (a TM-wave) gives 

,2 



A-^ - —F 



{kmn, ou) - — x™„ (-r) [k^n + T,u)= (225) 

{kmn + r) -Eq (krnn + r,u'^ ^Xmn (^) (kmn, = 0. 



Since the term containing (^co — (ji + uj is small when {u — kuj vanishes, 
in the second equation it is dropped. 

The dispersion equation defining the relation between kz and u is obtained by 
equating to zero the determinant of the system ( 1225^ and has a form: 



{uj — k^uf 



(j^mn ~ ^^o) ({kmn + ^) — ) ^ ^XtX- 

^ ^ ' (226) 



!d. ) A 

7c4 1 27r 



Because the right-hand side of the equation is small, one can again seek the 
solution near the points where the right-hand side is zero that corresponds 
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the condition of occurrence of the Cherenkov radiation and excitation of the 
wave which can propagate in the waveguide: 



(227) 



The roots of equations are sought near the conditions A;^„ ~ {kmn + 't) , 



(A;, + r,)' = [(A;o^ + r,) + = (/co. + + 2 (/cq^ + t,) ^ + f 
Hence, 

(A;o^ + T^)' + («mn + Tx)' + 2 (/Cq. + T,) + 2 (/Cq, + T,) ^ + = 

(229) 

{kmn + rf + 2 {ko, + T,) e + = ^Ln + 2^0mnT + + 2 (/Cq, + T,) ^ + C'- 

And one can get 

2A;o,e (2 (/^oz + T,) e + {'^komnT + —XrX-r = 



4A;oz (A;o. + r,) f + 2A;o. (24^„f + r") ^ - ^XrX-r = (230) 



+ ^ - z = 

(^0^ + T^) (/COz + T;^) 



_ (2feV+r^) 
^1.2 - -T7I — ± 



2 

2koT + r'^ \ _^ XrX-r 



4 {koz + T,) ^ {kQ, + T^) y c4 4A;oz (A;oz + r^,) 
If {koz + Tz) = — l^oz + TzI, the root can cross the zero point. At the same 
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time, the second equation should hold 
ui — kzU — uj — kozU — — 0. 



Consequently, 



^ = = koz^ V^o 5— • 



If £0 < 1, then e = ^ (1 - /3 - > 0, e = ^ - fc, 



Let the roots i^i and ^2 coincide (i^i = ^2)- This is possible at point 

^{koz + Tz) c2 ^4koz\koz + rz\' 
here koz + Tz < 0. 

The roots coincide when the following equality is fulfilled 

^ , _ VXtX-t 

f^oz — =F" 

u 



\l^kQz \kQz+Tz\ 



I.e., 



-- = k,,^^^^^ and k,.^^-L-^^ 



sj^k^z \kQz + Tz 

Let £0 < 1, then ^ > k^z (since w < c), the situation for the solution ^ = 



fcoz — ^ / ^^^^ ^ gets complicated and the Vavilov-Cherenkov condition is 
not fulfilled. 



Now let us consider the solution = /cq^ + ^ Vx^x r — At < the 
difference /cqz + can be reduced so that the sum on the right would appear 
to become equal to -, and so one could obtain 4 coinciding roots. 
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Interestingly enough, for backward diffraction, which is a typical case of fre- 
quently used one-dimensional generators with a corrugated metal waveguide 
(the traveling- wave tube, the backward- wave tube), such a coincidence of roots 
is impossible. 

Indeed, let the roots .^i and .^2 coincide. In this case for the backward Bragg 
diffraction \tz\ ~ S/cqz, < 0. Then by substituting the expressions for 
fcoz = ^y^gQ — '^"'jf and Eq = = 1 + Xmn (0) ^^"^ retaining the first-order 
infinitesimal terms, the relation - k^z + can be reduced to the form 

M|H _ + i^lM^ < the equality does not hold and 

the four-fold degeneracy is impossible. Only ordinary three-fold degeneration 
is possible. 

However, ii Eq > 1 and is appreciably large, then in a one- dimensional case, 
the four-fold degeneracy of roots is also possible in a finite photonic crystal^. 

Thus, the left-hand side of equation f l226p has four roots ^1, ^2, and a double 
degenerated root ^3. Hence, equation ( 12261) can be written as follows: 



If the roots coincide (^1 = ^2 = ^3), one obtains — ^i)*^ = B,, i.e., ^ — ^1 = 




The fourth root of B has imaginary solutions depending on the beam density 
as Imkz ~ pI^^ (the parameter B ~ col, i.e., B ~ pq, see the right-hand side 

^ The authors are grateful to K. Batrakov, who drew our attention to the fact 
that for an infinite crystal with Eq > 1, the intersection of roots is possible in a 
one-dimensional case. 
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of fl226p ). This increment is larger than the one we obtained for the case of 
the three-fold degeneracy. 

The analysis shows that with increasing number of diffracted waves, the law 
established in [25|49fl86] is valid: the instability increment appears to be pro- 
portional to p~, where s is the number of waves emerging through diffraction. 
As a result, the abrupt decrease in the threshold generation current also re- 
mains in this case (the threshold generation current jth ~ furV-i/l — where 
L is the length of the interaction area). 

It is interesting that according to |117j . for a photonic crystal made from 
metallic threads, the coefficients x (''")) defining the threshold current and the 
growth of the beam instability, are practically independent on r up to the 
terahertz range of frequencies because the diameter of the thread can easily 
be made much smaller than the wavelength. That is why photonic crystals 
with the period of about 1 mm can be used for lasing in terahertz range at 
high harmonics (for example, photonic crystal with 3 mm period provides the 
frequency of the tenth harmonic of about 1 terahertz (A=300 micron). 

The analysis of laser generation in VFEL with a photonic crystal when the 
beam moves in an undulator (electromagnetic wave) located in a finite crystal, 
made similarly to the above analysis, shows that in this case the dispersion 
equation and the law of instability also have the same form as in the case 
of an infinite crystal. The procedure for going from the dispersion equations 
describing instability in the infinite case f l217p is similar to that discussed 
earlier in this paper. It consists in replacing the continuous k by the quantified 
values of kmn and redefining the coefficients appearing in equations like (I217p . 

It is important to emphasize the general character of the rules found in this 



173 



paper for obtaining dispersion equations that describe the radiative instabihty 
of the electron beam in a finite photonic crystal. In particular, they are valid 
for describing the processes of instability of an electromagnetic wave in finite 
nonlinear photonic crystals. 

25 Hybrid systems with virtual cathode for high power microwaves 
generation ( |166| ) 

The interest for high power microwave (HPM) sources has emerged in re- 
cent years due to revealing new applications and offering novel approaches to 
existing ones. 

Vacuum electronic sources, which convert the kinetic energy from an electron 
beam into the electromagnetic field energy, are a natural choice for generating 



The high current density electron beam, once generated, propagates through 
an interaction region, which converts the beam's kinetic energy to HPM. It 
is the particular nature of the interaction that distinguishes various classes of 
sources. 

High power HPM sources generating high electromagnetic power density re- 
quire the high power densities in the electron beams, where space-charge effects 
are essential. 

When the magnitude of the current of an electron beam injected into a 
drift tube exceeds the space-charge-limiting current Iiimit, an oscillating virtual 
cathode (VC) is formed |167yi68j . 




HPM. 
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According to |167j . the following formula gives a good approximation for the 
space-charge- limiting current: 



'■limit 



{kA) 



mc 



^(7^/3 - 1) 



3/2 



(231) 



where m and e are the mass and the charge of an electron, c is the speed 
of light, 7 is the Lorentz factor of the electron beam and G depends on the 
geometry [167fl69j . For example, for an annular electron beam in a cylindrical 
drift tube G reads as follows |169] : 

1 



G 



rt 



+ 21n^)(l-sec/if^ 



(232) 



where rj, and rl^ are the outer and inner radius of the electron beam, R and L 
are the radius and length of the cylindrical drift tube and Hi is the first root 
of the Bessel function Jo{fi) = 0. 

When the oscillating virtual cathode is formed, two types of electrons exist: 
those oscillating in the vircator potential well and passing through the vircator 



area (see Figure [261 
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Figure 26. An oscillator with virtual cathode 

For electrons oscillating in the area " cathode-anode- virtual cathode" , two ra- 
diation mechanisms provide a radio- frequency signal |172] : 
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1. one radiation mechanism originates from the oscillations of the reflected 
electrons about the anode foil (electron oscillations in the potential well " cathode- 
anode- virtual cathode"). A microwave signal is generated at a frequency of 
roughly where dd is the anode-to-cathode spacing. 

2. The other radiation mechanism is the oscillation of the virtual cathode at 
a frequency near the plasma frequency ujp of the space charge density that is 
formed. That is 



where Ue is the number density of the electrons in the space charge configu- 
ration (in the plane of the anode grid) |172] . 

The essence of the above radiation mechanisms is bremsstrahlung radiation 
ensuing from electron deceleration. 

It should be emphasized that bremsstrahlung radiation from electrons oscillat- 
ing in an electron beam with a virtual cathode is accompanied by transition 
radiation, which is originated by electron velocity rather than acceleration. 
Use of a photonic crystal enables one to construct several types of hybrid 
systems with a virtual cathode, which could radiate due to different radia- 
tion mechanisms (bremsstrahlung and diffraction (transition) radiation) with 
different frequencies. 

In vircator systems, a grid cathode and anode (or anodes) are commonly used 
|169fl70fl71fl73] . The electron beam oscillates, making electrons periodically 




cross the grid anodes and cathode (see, for example |171j ). It is transition ra- 
diation that occurs when electrons pass through a border between two media 




(233) 
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with different indices of refraction. It is worth noting that periodical excite- 
ment of transition radiation from electrons oscillating in vircator is similar to 
diffracted radiation from a charged particle in a periodic structure. As a result, 
in a system with oscillating virtual cathode, the vircator radiation, which is 
actually the electron beam bremsstrahlung, is accompanied by radiation ex- 
cited by the additional mechanism due to transition (diffraction) radiation 
from oscillating current passing through the grid anodes (cathode). 



Let us turn to that part of the beam, which passes through the virtual cathode 
area. Recall that the oscillation of the virtual cathode can produce a highly 
modulated electron beam, and, as a result, the energy from the bunched trans- 
mitted beam can be recovered using slow- wave structures |174j . 




slow-wave 
structure 



Figure 27. Hybrid system "vircator + travelling wave tube" |174j 

Interaction of the electron beam with the slow-wave structure in, for instance, 
a conventional TWT poses special challenges: interaction is sufficiently effec- 
tive only for electrons moving at the distance 6 from the slow-wave structure 
surface 



6 < 



A/37 
47r ' 



(234) 



6 is the so-called beam impact parameter, A is the radiation wavelength, (3 
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v/c , V is the electron beam velocity, 7 is the electron Lorentz-factor. For 
example, for electrons with the energy of 250 keV (/3 = 0.74 and 7 = 1.49) 
and the radiation wavelength A = 10mm (frequency 30 GHz), the impact 
parameter 6 ~ 0.9mm. It means that for efficient radiation generation by the 
annular electron beam with the thickness A < 6 only a part of the beam would 
contribute to radiation. 

As we showed in [87|89fl03fll2fll7fl75fl76fl77] . this challenge can be over- 
come by applying a photonic crystal, formed by metal grids (grid traveling 
wave tube (grid TWT), grid volume free electron laser (grid VFEL)). 

What is more, in accordance with |178|179|180] . the application of metal in- 
serts (meshes, grids and so on) inside a resonator enables increasing the elec- 
tron beam limit current. 

Therefore, in the grid TWT (grid VFEL), the presence of the metal grid 
(photonic crystal) serves both for forming the resonator, where interaction 
of the beam and radiation occurs, and for potential balancing that makes it 
possible to increase the beam vacuum limit current. 

And for the grid TWT (grid VFEL) with the supercritical current, the elec- 
tron beam executes compound motion exciting two radiation mechanisms con- 
tributing to radiation: bremsstrahlung of oscillating electrons and diffraction 
(transition) radiation from downstream electrons interacting with the periodic 
grid structure(photonic crystal). 

This means that the hybrid system "vircator -|- grid TWT (grid VFEL)" arise 
by analogy with the system described in |T74j , where several vircators could 
appear due to the presence of several anode grids (see also |169j ). But in 
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contrast to the system |169j . the hybrid system "vircator + grid TWT (grid 
VFEL)" uses periodically placed grids with either constant |112|ll75|ll76j or 
variable period |177] . 

The frequency of diffracted radiation excited by an electron beam in a periodic 
structure with the period d is determined by the condition 

UJ — kn{k)v = TV , (235) 



where v is the electron beam velocity, r is the reciprocal lattice vector (|t| = 
n(k) is the refraction index of a periodic structure, p is an integer number 
(p = l,2,3,...). 

When the electron beam velocity v is parallel to the reciprocal lattice vector 
f, ([235]) reads [TT7] 

^ = In — n — M — m 

(2(1 — pn[oo, k) cosu) 

Essentially, a photonic crystal in the grid TWT (grid VFEL) is transparent 
to radiation as well as to an electron beam (see Figure |28|) . Moreover, several 
diffracted waves could exist in a photonic crystal (see Figure I25]) d. which 
makes it possible to introduce a feedback in such a system at the frequency 
of diffracted radiation and, hence, to couple several hybrid "vircator + grid 
TWT (VFEL)" generators making a phase-locked source, in which diffracted 
waves from one photonic crystal (grid resonator) excites oscillations in the 
neighbor resonators. 

The proposed grid systems drastically differ from the system described in 
|169] . where several grids serve only for forming several vircators (see Figure 
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Figure 28. Grid TWT (grid VFEL) and photonic crystal arrangement 




drift tubes 
Figure 29. Several vircators [169j 

Of course, radiation from a hybrid generator "vircator + grid TWT (VFEL)" 
can be excited by several electron beams similar to a phase-locked array. 

The bunched electron beam, which has passed through the virtual cathode 
area, can also be used for excitation of free electron laser (ubitron) (see Figure 
1301) oscillation contributing to the radiation power. 

Thus, the use of a photonic crystal enables one to construct several types 



180 




Figure 30. Hybrid system "vircator + undulator FEL" 
of hybrid systems with virtual cathode, which could radiate due to different 

radiation mechanisms (bremsstrahlung and diffraction (transition) radiation) 

with different frequencies. The photonic crystal makes it possible to create a 

phase-locked array of generators. 

26 Volume Free Electron Laser (VFEL) as a new trend in devel- 
opment of high-power tunable radiation sources: review of the- 
oretical background and experiments in millimeter range 

To conclude, let us mark the most important results obtained in investiga- 
tion of the VFEL. Use of non one-dimensional distributed feedback in vacuum 
electronic devices enables frequency tuning in a wide range and removes limits 
for available output power. It also solves a problem of development of highly 
stable mode-selective overmoded resonators, which open up new avenues for 
extending the operating frequencies of all classes of microwave vacuum elec- 
tronic devices. 

New advances in different areas require the development of tunable, wide- 
band, high-power sources of coherent electromagnetic radiation in gigahertz, 
terahertz and higher frequency ranges. Conventional electron vacuum devices 
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have restricted the possibihty of frequency tuning (usually it does not exceed 
5-10%) for the certain carrier frequency at certain e-beam energy. They also 
have power limits because high levels of power density inside the system cause 
sparkovers and damage of mirrors. 

Volume free electron laser (VFEL) [25)18711181] was proposed as a new type of 
a tunable high-power source of electromagnetic radiation. 

The most essential feature of the FEL and other types of generators is a feed- 
back, which is formed by a system of mirrors, or distributed feedback based 
on diffraction in a spatially periodic medium, when wave vectors of transmit- 
ted and reflected waves are colinear. The distinction of volume FEL (VFEL) 
is non-one- dimensional multi-wave volume distributed feedback (VDFB). Use 
of a non one-dimensional distributed feedback in vacuum electronic devices 
gives the possibility of frequency tuning in a wide range and removes limits 
for available output power. It also solves a problem of development of highly 
stable mode-selective overmoded resonators, which open up new avenues for 
extending the operating frequencies of all classes of microwave vacuum elec- 
tronic devices. 

It is well known that each radiative system is defined by its eigenmodes and 
by the so-called dispersion equation, which in the case of small perturbations 
(linear regime) describes possible types of waves in a system and relation 
between frequency and wave number of the system eigenmodes. Thorough 
analysis of FEL dispersion equation |108j shows that: 

1. dispersion equation for FEL in collective regime coincides with that for 
conventional traveling wave tube amplifier (TWTA) [109J; 

2. FEL gain (increment of electron beam instability) is proportional to p^^^. 
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where p is the electron beam density. 

But the law of instability of the electron beam can essentially change at passing 
through a spatially periodic medium. This fact was first indicated in [49j. Dis- 
persion equations were obtained and investigated for conditions of multi-wave 
diffraction. It was shown that there is a new law of electron beam instability at 
the points of diffraction roots degeneration. The amplification and generation 
gain of the electromagnetic wave are sharply changed at these points. 

This result is also valid for an electron beam, which moves in vacuum close 
to the surface of a spatially periodic medium [22] (or in a vacuum slot made 
inside a periodic medium). 

First lasing of the VFEL was reported at PEL 2001 [HI]. 

26. 1 Volume FEL distinctive features 

The advantages of VFEL are exhibited in a wide spectral range from mi- 
crowaves to X-rays [5^87|110] . Frequency tuning, possibility to use wide elec- 
tron beams (several e-beams) and reduction of threshold current density nec- 
essary for the start of generation provided by VFEL, make it a basis for 
development of more compact, high-power and tunable radiation sources than 
conventional electron vacuum devices could let. 

Benefits given by VFEL: 

1. Volume FEL provides frequency tuning by rotation of diffraction grating; 

2. Use of multi-wave diffraction reduces generation threshold and the size of 
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the generation zone. The starting current j depends on interaction length L as 
[55] : jstart ~ ^/ {{kLY^kXrL)'^'^}, s is the number of surplus waves appearing 
due to diffraction (for example, in the case of two-wave Bragg diffraction s = 1, 
for three- wave diffraction s = 2, and so on). 

3. Wide electron beams and diffraction gratings of large volumes can be used in 
VFEL. Two or three-dimensional diffraction gratings allow one to distribute 
the interaction over a large volume and to overcome power restrictions in 
a resonator. Volume distributed feedback provides mode discrimination in a 
VFEL resonator. 

4. VFEL can simultaneously generate radiation at several frequencies. 
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